Walmes,
Encontrei comandos que poderiam ser úteis para a CDF da skew t distribution,
para os meus dados. O problema é que não estou conseguindo operacionalizar
direito (eu acho).
Como "data" tenho 1824 valores de erros padronizados advindos de uma skew t
distribuition, nu e xi são os graus de liberdade e assimetria, respectivamente.
Uso o comando:
cdf = psstd(data, mean = 0, sd = 1, nu = 7.07895, xi = -0.15122)
o erro que aparece é:
Erro em pt(q, df, lower.tail, log.p) :
Argumento não-numérico para função matemática
Tentei ordenar então, da seguinte forma:
cdf = psstd(data[1:1824, ], mean = 0, sd = 1, nu = 7.07895, xi = -0.15122)
o que fez o comando funcionar, mas o resultado ficou bem estranho, com diversos
números muito próximos....devo não estar entendendo direito a rotina.
daí testei outra rotina
#### cdf= st.cumulants(data, scale = 1, shape = -0.15122, df = 7.07895, n=1)
Mas os resultados tbm formam "estranhos" , acredito que não devo estar usando
corretamente as funções....alguém saberia me dizer????
Atc
Julio
________________________________
De: Walmes Zeviani <[email protected]>
Para: [email protected]; julio cesar araujo <[email protected]>
Enviadas: Sábado, 5 de Novembro de 2011 22:21
Assunto: Re: [R-br] cdf da skewed-t - tem alguma rotina pronta??
Rode esse comando no R e navegue nos resultados
RSiteSearch("skew t distribution")
À disposição.
Walmes.
==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: [email protected]
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================_______________________________________________
R-br mailing list
[email protected]
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.
