On Tue, 30 Sep 2008, Susanne Pfeifer wrote:
Hi,
sorry, I was a bit confused and my function looked slightly different
from the one I posted before:
In any case, you will not get far trying to integrate this one for values
of z like 12000 without doing some homework. Simply plugging the func
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,
Hi Charles,
hi Victor
Charles C. Berry schrieb:
> What verson of R?
>
> Works for me:
>
>> 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,
I got the same error.
> sessionInfo()
R version 2.6.0 (2007-10-03)
i386-pc-mingw32
locale:
LC_COLLATE=Spanish_Colombia.1252;LC_CTYPE=Spanish_Colombia.1252;LC_MONETARY=Spanish_Colombia.1252;LC_NUMERIC=C;LC_TIME=Spanish_Colombia.1252
attached base packages:
[1] stats graphics grDevices utils
What verson of R?
Works for me:
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)
18.9513 with abs
Hello,
I am trying to use R's integrate function to calculate the following
integral for z=423:
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,
Hello,
I am trying to use R's integrate function to calculate the following
integral for z=423:
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,
7 matches
Mail list logo
