partage/GeneralizedAsymmetricLaplace
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Correction vectorisation du calcul

Browse source
This commit is contained in:
François Pelletier 2014-03-08 22:37:11 -05:00
parent a2ff837ec7
commit 2a94c8656e
1 changed files with 4 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

13
R/psaddleapproxGAL.R
View file

@ -22,13 +22,8 @@ psaddleapproxGAL <- function(x,param,eval.time=1,type="mu",log=FALSE)
u <- s * sqrt(diffcgfGAL(s,2,param,eval.time,type,log))
w <- sign(s)*sqrt(2*(s*x-cgfGAL(s,param,type,log)))
if(x==mGAL(1,param,type,log))
{
return(1/2 + diffcgfGAL(0,3,param,eval.time,type,log)/
(6*sqrt(2*pi)*diffcgfGAL(0,2,param,eval.time,type,log)^(3/2)))
}
else
{
return(pnorm(w)+dnorm(w)*(1/w-1/u))
}
(x==mGAL(param,1,type,log)) * (1/2 + diffcgfGAL(0,3,param,eval.time,type,log)/
(6*sqrt(2*pi)*diffcgfGAL(0,2,param,eval.time,type,log)^(3/2))) +
(x!=mGAL(param,1,type,log)) * pnorm(w)+dnorm(w)*(1/w-1/u)
}