Showing posts with label 3 vartiable plots. Show all posts
Showing posts with label 3 vartiable plots. Show all posts

Wednesday, May 1, 2013

RG#99: cloud 3D bars with heatmap


require(lattice)
require(latticeExtra)

data(VADeaths)

cloud(VADeaths, panel.3d.cloud = panel.3dbars,
      xbase = 0.4, ybase = 0.4, zlim = c(0, max(VADeaths)),
      scales = list(arrows = FALSE, just = "right"), xlab = NULL, ylab = NULL,
      col.facet = level.colors(VADeaths, at = do.breaks(range(VADeaths), 20),
                               col.regions = cm.colors,
                               colors = TRUE),
      colorkey = list(col = cm.colors, at = do.breaks(range(VADeaths), 20)),
      screen = list(z = 40, x = -30))


Sunday, April 7, 2013

RG#30: Level plot

require (lattice)
# data = volcano
levelplot(volcano)

#data
 xvr <- seq(-5, 5, len = 50)
yvr <- seq(-5, 5, len = 50)
gdr <- expand.grid(xvar = xvr, yvar = yvr)
gdr$zvar <- rnorm (nrow(gdr), 4, 1)


# color brewing 
library("RColorBrewer"); library("lattice");
brewer.div <- colorRampPalette(brewer.pal(11, "Spectral"), interpolate = "spline")

#plot 
levelplot (zvar ~ xvar * yvar, data = gdr, cuts = 199, col.regions = brewer.div(200), aspect = "iso")



# ggplot2 version 
require(ggplot2)
pl<- ggplot(gdr, aes(xvar, yvar, z= zvar))
pl + geom_tile(aes(fill = zvar)) + theme_bw()



RG#28: Contour plot


#data
 xvr <- seq(-5, 5, len = 50)
yvr <- seq(-5, 5, len = 50)
gdr <- expand.grid(xvar = xvr, yvar = yvr)
gdr$zvar <- rnorm (nrow(gdr), 4, 1)


#plot
require(lattice)
contourplot(zvar ~ xvar * yvar, data = gdr, cuts = 10)




#ggplot2 version
require(ggplot2)
plt <- ggplot(gdr, aes(xvar, yvar, z= zvar))
plt + stat_contour() + theme_bw()

# filled with zvar
 plt + geom_tile(aes(fill = zvar)) + stat_contour() + theme_bw()
plt + stat_contour(geom="polygon", aes(fill=..level..)) + theme_bw()






RG#22: heatmap plot using ggplot2

#data
set.seed(1234)
xm1 <- matrix(rnorm(100*10, rnorm(100, 0.5, 0.1)), nrow=100, ncol=10, byrow=FALSE)
xm2 <- matrix(rnorm(100*10, rnorm(100, 0.5, 0.1)), nrow=100, ncol=10, byrow=FALSE)
xm3 <- matrix(rnorm(100*10, rnorm(100, 0.5, 0.1)), nrow=100, ncol=10, byrow=FALSE)
dd <- cbind(xm1, xm2, xm3)


#heatmap plot
require(ggplot2)

# first need to reshape data to long form
require(reshape)
dd.melt <- data.frame(melt(dd) )
ggplot(dd.melt , aes(x=X1,y=X2, z= value)) + geom_tile(aes(fill= value)) + scale_fill_gradient(low="red", high="green") + theme_bw()



Friday, April 5, 2013

RG#14: 3D plot - wireframe plot

 # data
XV <- seq(-4, 4, len = 40)
 YV <- seq(-4, 4, len = 40)
datf <- expand.grid(xvar = XV, yvar = YV)
 datf$ZV<-  datf$xvar^2 + datf$yvar^2 - 2*datf$xvar* datf$yvar


# plot
require(lattice)
wireframe(ZV ~ xvar* yvar, data = datf)


wireframe(ZV ~ xvar* yvar, data = datf, drape = TRUE, perspective = FALSE, aspect = c(3,2), colorkey = TRUE)