To Jeff: a tip of the hat. I have another question after Jeff's solution:
On Jan 3, 2007, at 5:09 PM, Jeffrey Horner wrote: > Michael Kubovy <kubovy <at> virginia.edu> writes: >> I tried and it gave a strange result. See >> http://people.virginia.edu/~mk9y/mySite/twoGaussian.R >> and >> http://people.virginia.edu/~mk9y/mySite/twoGaussian.pdf > > > Try the image function. The smoothness of the plot will be > proportional to the length of x and y. For instance 200 isn't bad: > > mu1 <- 0 > mu2 <- 5 > s <- 1 > x <- seq(-2.5, 7.5, length = 200) > y <- seq(-2.5, 2.5, length = 200) > f <- function(x,y){ > term1 <- 1/(2*pi*sqrt(s*s)) > term2 <- -1/2 > term3 <- (x - mu1)^2/s > term4 <- (y - mu1)^2/s > term5 <- (x - mu2)^2/s > term1*(.5 * exp(term2*(term3 + term4)) + .5 * exp(term2*(term5 > + term4))) > } # setting up the function of the multivariate normal density > z <- outer(x, y, f) > # persp(x, y, z) > require(grDevices) > #pdf('twoGaussian.pdf') > #filled.contour(x, y, z, axes = F, frame.plot = F, asp = 1, > # col = gray(seq(0, 0.9, len = 25)), nlevels = 25) > image(x,y,z,col=gray(seq(0,0.9,len=200))) Is there a simpler way to get rid of axes, frame, and axis labels than image(x, y, z, col = gray(seq(0, 0.9, len = 200)), asp = 1, xaxt = 'n', yaxt = 'n', bty = 'n', ann = F) ? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903 Office: B011 +1-434-982-4729 Lab: B019 +1-434-982-4751 Fax: +1-434-982-4766 WWW: http://www.people.virginia.edu/~mk9y/ ______________________________________________ R-help@stat.math.ethz.ch 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.
