[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