Displaying 1 result from an estimated 1 matches for "ex_sigma".
Did you mean:
estsigma
2013 Jan 28
1
Adding 95% contours around scatterplot points with ggplot2
...s.google.com/forum/?fromgroups=#!topic/ggplot2/85q4SQ9q3V8
using the mvtnorm package and the dmvnorm function, but haven't been able
to get it work for my data example. The calculated densities are always
zero (at this step of Hadley's code: dgrid$dens <-
dmvnorm(as.matrix(dgrid), ex_mu, ex_sigma) )
I appreciate any assistance.
Thanks,
Nate
x<-c(seq(0.15,0.4,length.out=30),seq(0.2,0.6,length.out=30),
seq(0.4,0.6,length.out=30))
y<-c(0.55,x[1:29]+0.2*rnorm(29,0.4,0.3),x[31:60]*rnorm(30,0.3,0.1),x[61:90]*rnorm(30,0.4,0.25))
data<-data.frame(level=c(rep(1, 30),rep(2,30), rep(3,30...