Message-ID: <48E28BE6.1070004@tiffy.it>
Date: 2008-09-30T20:28:22Z
From: Susanne Pfeifer
Subject: R's integrate function
In-Reply-To: <48E2892E.9090205@tiffy.it>
Hi,
sorry, I was a bit confused and my function looked slightly different
from the one I posted before:
Susanne Pfeifer schrieb:
>>> integrate(function(y,z){
>> + sapply(y, function(y,z){
>> + integrate(function(x,z)
>> +
>> 1/x*dbeta(0.01,x/(0.005/1.005),(1-x)/(0.005/1.005))*dbeta(y,x/(0.005/1.005),(1-x)/(0.005/1.005))*(1-y)^z,0,1
>>
>> + ,423)$value
>> + })
>> + },0,1,423)
integrate(function(y,z){
sapply(y, function(y,z){
integrate(function(x,z)
1/x*dbeta(0.01,x/(0.01/1.01),(1-x)/(0.01/1.01))*dbeta(y,x/(0.01/1.01),(1-x)/(0.01/1.01))*(1-y)^z,0,1,12000)$value
})
},0,1,12000)
Sorry for the confusion...