[R] integrate fails with errors

herodote@oreka.com herodote at oreka.com
Tue Jul 19 18:11:37 CEST 2005


Hi all,
i'm new to R,
I need to modelize in R a statistic algorithm,
This algo use Weibull, normal law, linear regression, normalisation, root mean square, to find eta and beta fitting the weibull model (to analyse few results) and further when we will get more information apply bayes model .

the problem is when When i try to integrate it fails with errors.

by the way i like to integrate something like that :

beta0,eta0,n are initialized as single integer, temp is a 1 dimension array containing 9 integer.

integrate(function(beta) ((beta/(eta0)^beta)^n)*prod(temp^(1-beta)*exp(-sum(temp^beta)/(eta^beta)))*(1/(sqrt(2*pi))*exp(((beta-beta0)^2,0,Inf)

R says : The longuest object isn't a multiple of the shortest object in temp^beta , the same for temp^(1-beta).

I don't understand why R fails to take calculate each temp^beta then sum it over again for each beta values.

don't know if my post is explicit my head is burned out with these things today

thks.

////////////////////////////////////////////////////////////
// Webmail Oreka : http://www.oreka.com
////////////////////////////////////////////////////////////




More information about the R-help mailing list