Hi everybody!
I´m with a problem that probably is easy for you, but I really don´t know
how to solve.
On the following script:
for(j in 1:length(limiares))
{
excessos<-limiares[j]-estacao[estacao<limiares[j]]
par.ests<-gpd(-(estacao),threshold=-limiares[j],method=c("pwm"))$par.est
GOF.test<-ks.test(excessos,"pgpd",xi=par.ests[1],beta=par.ests[2])$p.value
tabs[j,]<-c(par.ests,gpd(-(estacao),threshold=-limiares[j],method=c("pwm"))$par.ses,GOF.test,length(excessos))
}
I´ve found the error for some values of i: Erro em ks.test(excessos, "pgpd",
xi = par.ests[1], beta = par.ests[2]) : NA/NaN/Inf em chamada de função
externa (argumento 1)
*My question is: This warning stop the "for" and I don´t want it, is there
some way to continue the "for", and for the cases where the function cannot
calculate the ks.test for the i just leave a NA as answer???*
Thank you *very much*!!!
Alexandra Almeida
--
Alexandra R M de Almeida
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
