Hi, my R friends,
I am going to plot the surface of a bivariate normal distribution and its
contours. I have written the below codes:
library(MASS)
set.seed(69)
n <- 5000
x <- rnorm(n, 0, 15)
y <- 0.5 * x + rnorm(n, 0, 10)
z <- kde2d(x, y, n = 50)
persp(z, theta = 55, phi = 35,
shade = 0.75, col = "gold", expand = 0.5, r = 2,
ltheta = 25, ticktype = "detailed")
contour(z)
It seems that I could use mnormt package for bivariate normal distribution.
How can I use suitable functions of mnormt package in my above codes?
Indeed I want to change my above codes somewhat and use the mnormt package.
Thank you in advance
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
