On Mon, May 9, 2011 at 8:23 PM, Dat Mai <dat.d....@gmail.com> wrote: > Thanks for the help! Though, when I do this, I get a message saying "Error > using packet 1 value where TRUE/FALSE." Would that have anything to do with > this set of code?
It might. (As we don't know what you "did", we can't really say.) > On Sun, May 8, 2011 at 3:00 PM, Deepayan Sarkar <deepayan.sar...@gmail.com> > wrote: >> >> On Fri, May 6, 2011 at 10:21 PM, Dat Mai <dat.d....@gmail.com> wrote: >> > Hello all, >> > >> > I'm trying to create a heatmap using 2 matrices I have: z and v. Both >> > matrices represent different correlations for the same independent >> > variables. The problem I have is that I wish to have the values from >> > matrix >> > z to be represented by color intensity while having the values from >> > matrix v >> > to be represented by circle size. I currently have the following in >> > front of >> > me and an unsure of what to add or change in order to achieve that goal. >> > >> > panel.corrgram.2 = >> > function(x, y, z, subscripts, at = pretty(z), scale = 0.8, ...) >> > { >> > require("grid", quietly = TRUE) >> > x <- as.numeric(x)[subscripts] >> > y <- as.numeric(y)[subscripts] >> > z <- as.numeric(z)[subscripts] >> > zcol <- level.colors(z, at = at, ...) >> > for (i in seq(along = z)) >> > { >> > lims <- range(0, z[i]) >> > tval <- 2 * base::pi * seq(from = lims[1], to = lims[2], by = 0.01) >> > grid.polygon(x = x[i] + .5 * scale * c(0, sin(tval)), >> > y = y[i] + .5 * scale * c(0, cos(tval)), >> > default.units = "native", >> > gp = gpar(fill = zcol[i])) >> > grid.circle(x = x[i], y = y[i], r = .5 * scale, >> > default.units = "native") >> > } >> > } >> >> Probably this modification >> >> panel.corrgram.2 <- >> function(x, y, z, v, subscripts, at = pretty(v), scale = 0.8, ...) >> { >> require("grid", quietly = TRUE) >> x <- as.numeric(x)[subscripts] >> y <- as.numeric(y)[subscripts] >> z <- as.numeric(z)[subscripts] >> zcol <- level.colors(v, at = at, ...) >> >> <Rest unchanged> >> >> >> followed by >> >> levelplot(signif(z,3), v = v, panel = panel.corrgram.2, ...) >> >> will do it. >> >> -Deepayan >> >> >> > >> > k=levelplot(signif(z,3), scales = list(x = list(rot = 45)), >> > col.regions=col.sch, panel = panel.corrgram.2, >> > label = num.circ, xlab="", ylab="", main=paste(output,"receptor >> > response")) >> > print(k) >> > >> > -- >> > Best, >> > Dat Mai >> > PhD Rotation Student >> > Albert Einstein College of Medicine >> > >> > [[alternative HTML version deleted]] >> > >> > ______________________________________________ >> > R-help@r-project.org 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. ______________________________________________ R-help@r-project.org 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.