Re: [R] integration of two normal density

2010-06-28 Thread Bert Gunter
: Re: [R] integration of two normal density Your intuition is wrong and R is right. Why should the product of two probability density functions be a normalized pdf also? -- as is trivially seen with two uniforms on [0,2], with pdf= 1/2, product = 1/4 on [0,2] . -- Bert -Original Message

Re: [R] integration of two normal density

2010-06-28 Thread Matt Shotwell
Biostatistics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of bill.venab...@csiro.au Sent: Friday, June 25, 2010 10:53 PM To: carrieands...@gmail.com; R-help@r-project.org Subject: Re: [R] integration of two normal density Your

Re: [R] integration of two normal density

2010-06-26 Thread Ravi Varadhan
Ph. (410) 502-2619 email: rvarad...@jhmi.edu - Original Message - From: Carrie Li carrieands...@gmail.com Date: Friday, June 25, 2010 11:29 pm Subject: [R] integration of two normal density To: r-help R-help@r-project.org Hello everyone, I have a question about integration of two

Re: [R] integration of two normal density

2010-06-26 Thread Matt Shotwell
On Fri, 2010-06-25 at 23:28 -0400, Carrie Li wrote: Hello everyone, I have a question about integration of two density function Intuitively, I think the value after integration should be 1, but they are not. Am I missing something here ? t=function(y){dnorm(y, mean=3)*dnorm(y/2,

[R] integration of two normal density

2010-06-25 Thread Carrie Li
Hello everyone, I have a question about integration of two density function Intuitively, I think the value after integration should be 1, but they are not. Am I missing something here ? t=function(y){dnorm(y, mean=3)*dnorm(y/2, mean=1.5)} integrate(t, -Inf, Inf) 0.3568248 with absolute error

Re: [R] integration of two normal density

2010-06-25 Thread Bill.Venables
-help Subject: [R] integration of two normal density Hello everyone, I have a question about integration of two density function Intuitively, I think the value after integration should be 1, but they are not. Am I missing something here ? t - function(y){dnorm(y, mean=3)*dnorm(y/2, mean=1.5