Displaying 1 result from an estimated 1 matches for "geddenstandard".
2007 Feb 27
1
Additional args to fun in integrate() not found?
...e
cdf<-matrix(0, length(zvec),1)
for(i in 1:length(zvec))
{
z<-zvec[i]
if(numint==0)
{
if(z<=0)
{
t<-0.5*(1-1/gamma(1/p)*Igamma((1/p),(-z)^p,lower=TRUE))
}
else
{
t<-1-(0.5*(1-1/gamma(1/p)*Igamma((1/p),(z)^p,lower=TRUE )))
}
}
else
{
t<-integrate(geddenstandard, -35, z, subdivisions=1000,
rel.tol=100*.Machine$double.eps, abs.tol=rel.tol,
stop.on.error=TRUE, keep.xy=FALSE, aux=NULL,p)
}
cdf[i]<-t
}
cdf
}
#-------------------------------------------------------------------------
geddenstandard = function(z,p)
{
f<-p/(2*gamma(1/p))*exp(-abs(z)...