[R] Interpreting R -results for Bivariate Normal

beetle2 samandbrendan at aapt.net.au
Mon Jun 8 03:00:14 CEST 2009


HI Guys,
I know that this forum is not for homework but I am trying to interpret R
output code.
I was just wondering if someone might be able to help.

I have been given the following.

For (X1,X2) distributed bivariate normal with parameters 

mu1 = 5.8  
mu2 = 5.3  
sd1 = sd2 = 0.2

and p = 0.6

The r-code and inpit/output are as follows

input

m <- 5.3 + 0.6*(6.3 - 5.8)
s  <- sqrt(0.2^2*(1-0.6^2))
q  <- seq(5.12,6.08,0.16)

print(rbind(q,pnorm(mean=m,sd=sd=s,q=q)))

output

q 5.1200  5.280   5.44   5.6     5.76     5.92    6.1
   0.0013  0.023   0.16   0.5     0.84     0.98    1

I have been asked to interpret 

E[X2|X1 = 6.3] and varE[X2|X1 = 6.3]


I take it that s<- = 0.16 is the standard variation

So I am assuming that varE[X2|X1 = 6.3] = 0.16^2 = .0256

m <- 5.3 + 0.6*(6.3 - 5.8) = 5.6    this the Expected value of E[X+Y]


I see from the output that this would be correct because the probability of
5.6 = 0.5


to interpret E[X2|X1 = 6.3] I can't see it in the output.  And I'm not sure
how to find the conditional probabilty from the output.

Any help would be greatly appreciated




-- 
View this message in context: http://www.nabble.com/Interpreting-R--results-for-Bivariate-Normal-tp23916967p23916967.html
Sent from the R help mailing list archive at Nabble.com.




More information about the R-help mailing list