[R-br] cdf da skewed-t - tem alguma rotina pronta??

julio cesar araujo julio_economia em yahoo.com.br
Terça Novembro 8 02:00:58 BRST 2011


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 <walmeszeviani em gmail.com>
Para: r-br em listas.c3sl.ufpr.br; julio cesar araujo <julio_economia em yahoo.com.br>
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: walmes em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111107/54261098/attachment.html>


Mais detalhes sobre a lista de discussão R-br