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...
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
