search for: jjmichael

---------- Forwarded Message ---------- Subject: Re: [R] conditional coloring of image labels Date: Friday 11 November 2005 1:04 pm From: jim holtman <jholtman at gmail.com> To: Jake Michaelson <jjmichael at cc.usu.edu> Use 'mtext': genes=cbind(ABC1=c(3,4,4,5,6,3), ABC2=c(4,3,4,7,7,8), ABC3=c(8,7,8,6,3,2)) ###plot the image image(1:nrow(genes), 1:ncol(genes), genes, axes = FALSE, xlab = "", ylab = "", col=cm.colors(256)) sig=cbind(name=c("ABC1", "ABC2...

I've been using R for a while under Mac OS X, which thanks to the R on OS X developers, is probably the best platform for learning R. I recently built a Linux box with a Pentium D processor, and am running an AMD64 port of Ubuntu with the SMP kernel. After setting up the basics on the box, I thought I'd install R as well. I've had a very hard time finding good documentation

Hi all, In heatmap's documentation, it mentions that the output value is actually an invisible list...how would one access this list? Thanks, Jake

I've been using R for a while under Mac OS X, which thanks to the R on OS X developers, is probably the best platform for learning R. I recently built a Linux box with a Pentium D processor, and am running an AMD64 port of Ubuntu with the SMP kernel. After setting up the basics on the box, I thought I'd install R as well. I've had a very hard time finding good documentation

I'm working on some heatmaps, and the person I'm working with would prefer a red-black-green color palette (red denoting gene induction and green denoting gene repression). Does such a palette exist already? If not, is there an easy way to create one? Thanks, Jake

Hi all, I'm trying to build R 2.1.1 on Ubuntu 5.04 i686-SMP. Configure goes well with: ./configure --with-BLAS --with-readline=no but once I run 'make', I get the following error: In file included from devX11.c:64: devX11.h:57:74: X11/Intrinsic.h: No such file or directory Any ideas? Thanks in advance, Jake

Hi all, I've got a set of gene expression data, and I'm plotting several heatmaps for subsets of the whole set. I'd like the heatmaps to have the same color distribution, so that comparisons may be made (roughly) across heatmaps; this would require that the color distribution and distance functions be based on the entire dataset, rather than on individual subsets. Does

Hi all, Does anyone know of a fairly easy way to "stretch" a heatmap vertically? I've got 42 arrays and would like to be able to see as many significant genes as possible (right now I can only get 50 genes with it still being readable). In some comparisons there are several hundred significant genes. I've fiddled with the "asp" argument, but that doesn't

I want to plot a heatmap without reordering the columns. This works fine in "heatmap": > heatmap(meanX[selected,], col=cm.colors(256), Colv=NA) But in "heatmap.2" I get: > heatmap.2(meanX[selected,], col=cm.colors(256), Colv=NA) Error in if (!is.logical(Colv) || Colv) ddc <- reorder(ddc, Colv) : missing value where TRUE/FALSE needed (Note that instructions for

Hi All, I've worked on some code to take a heatmap with 1000 row entries, and split this up into 20 pages, each with 50 rows from the original heatmap. I want to preserve the row order such that all 20 pages, if put together, would comprise the original heatmap. Here's what I've done: ##make the initial heatmap, with all 1000 rows and write it to an object 'heatAll'

Hi all, I've got some image plots on which I'd like to include some gene information (in the margins using mtext). Unfortunately, the description is rather long and will need to be wrapped to fit on several lines. From what I know about mtext, it's really only meant for single-line labels, not paragraphs. Here's some sample code of the idea I'm trying to accomplish:

Does anyone know how to extract residuals in lmer? Here's the error I get: > crop.lme=lmer(response~variety*irrigation*pesticide+(1|rep)+(1|rep: pesticide)+(1|rep:pesticide:irrigation), crop.data) > qqnorm(crop.lme) Error in qqnorm.default(crop.lme) : y is empty or has only NAs > resid(crop.lme) NULL Thanks! --Jake

Hi all, I am interested in plotting a heatmap of a set of genes. I would like the text labels of these genes to be colored red rather than black if they meet certain statistical criteria (using an if statement). I'm not sure how to change individual color labels without changing them all. Can anyone provide some insight on how to do this? Thanks in advance, Jake

Hi All, I'm taking an Experimental Design course this semester, and have spent many long hours trying to coax the professor's SAS examples into something that will work in R (I'd prefer that the things I learn not be tied to a license). It's been a long semester in that regard. One thing that has really frustrated me is that lme has an extremely counterintuitive way for

I did a very quick search of the archive and couldn't find a readily available answer to this one: I'd like to convert, for example: c("a", "b", "a", "b") to c(1, -1, 1, -1) In the case of the first vector, it may be any length, but will always only have two unique values. It must always be replaced by corresponding values of 1 and -1. Any

Quick question: how can I reverse the order of the rows in a matrix? i.e. make the last row first and the first row last, etc.?