[R] 3D Gaussian with different colors

Roberto rmoscetti at unitus.it
Wed Dec 26 08:58:48 CET 2012


I wrote, finding pieces of code in the web, a simple script to draw a 3D
gaussian plot.
For the next step, I need that the 2 gaussian distributions have different
colors.

Can someone help me to do this?

fn <- function(x, y, scale, scale2)
dnorm(x,mean=1,sd=scale)*dnorm(y,mean=-1,sd=scale) +
dnorm(x,mean=2,sd=scale2)*dnorm(y,mean=1,sd=scale2)

x <- seq(-4,4,len=40)
y <- seq(-4,4,len=40)
z <- outer(x, y, fn, scale = 1, scale2 = 0.5)
nrz <- nrow(z)
ncz <- ncol(z)
jet.colors <- colorRampPalette( c("lightblue", "green", "darkgreen") )
nbcol <- 100
color <- jet.colors(nbcol)
zfacet <- z[-1, -1] + z[-1, -ncz] + z[-nrz, -1] + z[-nrz, -ncz]
# Recode facet z-values into color indices
facetcol <- cut(zfacet, nbcol)
persp(x,y,z, theta = 30, phi = 30, expand = 0.2, col = color[facetcol],
      ltheta = 120, shade = 0.75, ticktype = "detailed") 




--
View this message in context: http://r.789695.n4.nabble.com/3D-Gaussian-with-different-colors-tp4653938.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list