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. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.