Valeu Paulo.<br><br><div class="gmail_quote">Em 27 de abril de 2011 00:42, Paulo Justiniano <span dir="ltr"><<a href="mailto:paulojus@leg.ufpr.br">paulojus@leg.ufpr.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Entre parentesis está o erro padrao da estimativa, tipicamente obtido pelo hessiano numérico (exceto pelos casos em que o estimador tem forma fechada)<div><div></div><div class="h5"><br>
<br>
<br>
On Tue, 26 Apr 2011, Cristiano Melo wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Walmes,<br>
<br>
Fiz o seguinte:<br>
library(MASS)<br>
fitdistr(vetor_observado,"weibull") - que acredito ser suficiente.<br>
<br>
O resultado foi:<br>
shape scale<br>
1.377413 210.784742<br>
(0.167581) (23.479204)<br>
<br>
Que novamente acredito ser o que quero. Só não entendi esta informação entre parênteses.<br>
<br>
Em 26 de abril de 2011 09:21, Walmes Zeviani <<a href="mailto:walmeszeviani@gmail.com" target="_blank">walmeszeviani@gmail.com</a>> escreveu:<br>
Cristiano,<br>
<br>
A mensagem sobre os empates vem do fato da sua amostra possuir valores repetidos. Às vezes, a precisão<br>
dessa medidas (casas decimais) é pequeno, imagine medir altura de 100 pessoas, é bem provável ter duas<br>
com 1,78 m, ou outro valor.<br>
<br>
Na ks.test(vetor_observado, distribuição, parametro1, parametro2, demais_opções), você precisa passar o<br>
valor dos parâmetros sob hipótese. Normalmente os valores usando são as estimativas obtidas com os<br>
dados. Então você precisa estimar. Para o caso da normal, mean(x) e sd(x) são os estimadores. Para<br>
outras distribuições você pode usar a função MASS::fitdistr(). Consulte a documentação para instruções<br>
de uso.<br>
<br>
À disposição.<br>
Walmes.<br>
<br>
==========================================================================<br>
Walmes Marques Zeviani<br>
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)<br>
Departamento de Estatística - Universidade Federal do Paraná<br>
fone: (+55) 41 3361 3573<br>
VoIP: (3361 3600) 1053 1173<br>
e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a><br>
twitter: @walmeszeviani<br>
homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a><br>
linux user number: 531218<br>
==========================================================================<br>
<br>
<br>
2011/4/25 Cristiano Melo <<a href="mailto:cristianogmelo@gmail.com" target="_blank">cristianogmelo@gmail.com</a>><br>
É o seguinte: tenho em um arquivo txt um vetor que representa tempos até a falha de equipamentos.<br>
Gostaria de fazer alguns teste de aderência para verificar se estes dados se aproximam de algumas<br>
distribuições de probabilidade. Usei o lillie.test(dados) para verificar se dos dados aderem a uma<br>
distribuição nomal. No entanto, gostaria de verificar se estes mesmos dados (e algumas variações) se<br>
aderem a uma exponencial, gamma e weibull.<br>
<br>
Sei que a função é a ks.test(x, y,..., alternative=c("two.sided" "less" or "greater")) para o teste de<br>
Kolmogorov-Smirnov.<br>
Para testar normalidade com a função ks, fiz o seguinte: ks.test(vetor, "pnorm", sd=sd(vetor),<br>
mean=mean(vetor),alternative=c("two.sided")). Curioso que o resultado foi bem diferente da lillie.test,<br>
a seguinte mensagem foi apresentada:<br>
Warning message:<br>
In ks.test(vetor, "pnorm", sd=sd(vetor), mean=mean(vetor),alternative=c("two.sided")), : não é possível<br>
calcular os níveis descritivos corretos com empates.<br>
O que isso quer dizer???????<br>
<br>
Quando tentei usar outra distribuição, pweibull por exemplo, o p-value foi menor que 2.2e-16, ou seja,<br>
nada a ver, e repetindo a mesma frase anterior. O mesmo resultado foi com as outras. Como não sei a<br>
sintaxe para weibull fiz o seguinte:<br>
ks.test(vetor, "pweibull", 1.129, 2,alternative=c("two.sided"))<br>
<br>
Onde estou errando? Vi que para montar uma pweibull são necessários o vetor de quantis e os parâmetros<br>
shape e scale. É necessário fazer separado e quardar em uma variável e depois jogar na ks? Como consigo<br>
esse vetor de quantis? Achei que seria automático.<br>
Estou correto se fizer assim:<br>
<br>
ks.test(vetor, "pweibull",10,2)<br>
<br>
E tem alguma forma de estimar os parâmetros shape e scale desta função?<br>
<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
<br>
<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
<br>
<br>
<br>
</blockquote>
</div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
<br></blockquote></div><br><div style="visibility: hidden; left: -5000px; position: absolute; z-index: 9999; padding: 0px; margin-left: 0px; margin-top: 0px; overflow: hidden; word-wrap: break-word; color: black; font-size: 10px; text-align: left; line-height: 130%;" id="avg_ls_inline_popup">
</div>