Re: [R] Bayesian Integration (stat question)

From: francogrex <>
Date: Sun, 30 Dec 2007 11:16:14 -0800 (PST)

Yacas and Ryacas does exactly what I want. Actually I had asked previously a question about algebra (not integration precisely) and guys here have pointed out package Ryacas. Sorry to disturb.

francogrex wrote:
> This may be far-fetched:
> In Bayesian analysis to find the marginal posterior distribution of a
> parameter it requires integration out of the so-called nuisance
> parameters. Is there in R (like an equivalent to the deriv function) but
> to find the expression of the integration (or an anti-derivative)? Like
> for example I want to integrate out mu in the normal distribution, I would
> introduce (note:punctuation marks are made up):
> integral(prod[i to n] 1/sqrt(2*pi*s^2)*exp(-(xi-mu)/2*s^2), dmu)
> and it would give me:
> 1/(sqrt(n)*(2*pi*s^2)^((n-1)/2))*exp(-sum[i to n]((xi-xbar)^2)/2*s^2)
> ?

View this message in context:
Sent from the R help mailing list archive at

______________________________________________ mailing list
PLEASE do read the posting guide
and provide commented, minimal, self-contained, reproducible code.
Received on Sun 30 Dec 2007 - 19:21:14 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sun 30 Dec 2007 - 19:31:24 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.