From: Thomas Lumley <tlumley_at_u.washington.edu>

Date: Wed 20 Jul 2005 - 03:53:46 EST

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jul 20 03:59:53 2005

Date: Wed 20 Jul 2005 - 03:53:46 EST

On Tue, 19 Jul 2005, [iso-8859-1] herodote@oreka.com wrote:

*>
*

> 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.
**>
*

The help page for integrate() says

f: an R function taking a numeric first argument and returning a numeric vector of the same length.

Your f is not of that form.

-thomas

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jul 20 03:59:53 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:33:50 EST
*