seu codigo geral valores negativos parauma variavel weibull e coloca
isto em um vetor (e nao uma lista) no caso de eixtirem NA's
portanto nao tem como funcionar memso

sugiro revisar seus objetos examinando os conteúdos entes de rodar a
fitdistr para ver se de fato posseui dados compativeis com isto no
objeto passado a fitdistr



 Dom, 2012-02-19 às 20:38 +0000, Vanúcia Schumacher escreveu:
> Pessoal,
> 
> 
> estou tentando gerar os parametros de weibull no R, mas não estou
> conseguindo
> Conforme orientaçao do Bruno obtive um erro, e não estou entendendo.
> alguem pode ajudar?
> 
> 
> library(MASS)
> > dados<-read.table("teste.txt",header=FALSE)
> > is.na(dados) <- dados == -999.9
> > vento50<-fitdistr(dados[[1]],densfun="weibull")
> Mensagens de aviso perdidas:
> In dweibull(x, shape, scale, log) : NaNs produzidos
> > is.na(dados) <- dados == -999
> > vento50<-fitdistr(dados[[1]],densfun="weibull")
> Mensagens de aviso perdidas:
> In dweibull(x, shape, scale, log) : NaNs produzidos
> > str(vento50)
> List of 5
>  $ estimate: Named num [1:2] 1.39 2.11
>   ..- attr(*, "names")= chr [1:2] "shape" "scale"
>  $ sd      : Named num [1:2] 0.191 0.292
>   ..- attr(*, "names")= chr [1:2] "shape" "scale"
>  $ vcov    : num [1:2, 1:2] 0.0366 0.0184 0.0184 0.0854
>   ..- attr(*, "dimnames")=List of 2
>   .. ..$ : chr [1:2] "shape" "scale"
>   .. ..$ : chr [1:2] "shape" "scale"
>  $ loglik  : num -46.9
>  $ n       : int 30
>  - attr(*, "class")= chr "fitdistr"
> > 
> 
> 
> 
> 
> 
> _______________________________________________
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.