[R] something weird in integration (pracma library)
casperyc
casperyc at hotmail.co.uk
Sat May 5 18:42:53 CEST 2012
Hi,
################################################################################
library(pracma)
k=20
mu=4.5
casigma=17000
myint=function(j) {
quadinf(function(x)
(1/(1+exp(-x)))^j*(1-1/(1+exp(-x)))^(k-j)*dnorm(x,mu,casigma),-Inf,Inf)
}
sapply(0:k,myint)
################################################################################
works fine
################################################################################
casigma=50500
sapply(0:k,myint)
################################################################################
casigma too large!
so try change of variable, y= (x-mu)/sigma
################################################################################
myint3=function(j) {
quadinf(function(y)
(1/(1+exp(-y*casigma-mu)))^j*(1-1/(1+exp(-y*casigma-mu)))^(k-j)*dnorm(y),-Inf,Inf)
}
sapply(0:k,myint3)
################################################################################
works again, but maybe precision is reduced??
HOWEVER, the problem now is
################################################################################
casigma=101
sapply(0:k,myint3)
casigma=100
sapply(0:k,myint3)
casigma=99
sapply(0:k,myint3)
casigma=98
sapply(0:k,myint3)
casigma=97
sapply(0:k,myint3)
################################################################################
does NOT work when casigma is 99 or 100. (when casigma is 'small')
I wonder if there are 'many' other small values of casigma that have the
same problem???
and why???
Casper
-----
######################
PhD candidate in Statistics
Big R Fan
Big LEGO Fan
Big sTaTs Fan
######################
--
View this message in context: http://r.789695.n4.nabble.com/something-weird-in-integration-pracma-library-tp4611381.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list