similar to: gridBase and heatmap

Dear all, I have a simple plot using "gridBase" like this. The problem occurs whenever I execute this code there is always a blank page created before the actual plot. How can we disable that blank page? I am using: R version 2.7.2 (2008-08-25) and gridBase version: 0.4-3 __ BEGIN__ library(grid) library(gridBase) opar <- par(no.readonly=TRUE) par(opar) grid.newpage()

I have written a wrapper for heatmap.2 called heatmap.w.row.and.col.clust which auto-generates breaks using breaks<-round((c(seq(from=(-20 * stddev), to=(20 * stddev))))/20, digits = 2) #(stddev in this case = 2.5) This has always worked well in the past but now I am getting an error that non-finite breaks are being generated. Drilling down, it seems that my wrapper is generating finite

Hi all, I am learning R by "doing". And this is my first post. I want to use R: 1- to fetch a DNA sequence from a databank (see bellow) and 2- store it as FASTA file. The problem: neither an error is prompted nor the fasta file is created. Testing the code (see bellow), I notice that everything works until the *"write.dna" *command - which is not creating the fasta file.

I'm suppressing the labeling of my rows and columns in heatmap.2 using the commands: labRow = " ", labCol = " " But I'd like to annotate them again using the axis command: axis(1, at=seq(500, 1000, 500)) mtext("Group 1", "Group 2") For some reason however it appears that the axis command is having no effect. My complete

I suspect you meant WD <- "~/Documents/Scripting/R_Studio/Sequences/" but I am entirely unfamiliar with the packages you are using, and know nothing about what is on your hard drive. For future reference: A) Read the Posting Guide. This is a plain text email list, and your html formatting gets removed leaving a mess that is not always readable. B) Most frequent users of R

The function heatmap() allows to specify row/column labels via the options labRow/labCol. From the code of heatmap(), I understand that when no labels are specified, the row/column labels (or indices) of the input matrix are taken as labels and re-ordered together with the rows and columns of the matrix before plotting, whereas labels supplied via labRow/labCol are plotted in the original order.

I am using attributes colsep and sepwidth from heatmap.2 function, however when using larger values for sepwidth the following columns are not moved to the right but the seperator is covering the column(s) with a white bar and the column labels and ColSideColors remain at their positions, is this intended? I expected the following columns to be shifted to the right by sepwidth. Is there a way

Dear list, I have produced a fairly intricate plot arrangement for use in a publication using layout() and gridBase, and out of curiosity I'd like to learn whether a more elegant and robust solution could be obtained with grid to avoid the layout() function. > library(gridBase) > x <- seq(-pi, pi, length=10) > y <- cos(x) > > par(mar = c(0, 0, 0, 0), pty =

All, I have a simple plot(x,y) and I would like to then insert rectangles of some length (in native coordinates) and height fixed to 0.5 in native coordinates. I can't quite get the code right to do this. Can anyone give me a quick example of how to do this? I looked the gridBase index and the tutorial (from R-news?) but just haven't gotten it down yet. > plot(1:10,1:10)

Hello, I'm trying to follow the documentation of how to use gridBase, and I've reached the minimal code example below as my best effort. Can someone explain how to keep the column of boxplots on the same page as the rectangles (even though I've tried new = TRUE) ? Also, would it be hard / possible to match up the middle of each boxplot to the middle of each rectangle ?

I am running R 2.2.1 on a University-supported linux installation based on Redhat EL3. I am sorry that it did not occur to me to mention this before; I updated R very recently, with the most recent version available for EL3 at http://cran.cnr.berkeley.edu/bin/linux/redhat/el3/. Looking at the gridBase documentation, I find that 2.2.1 is not in fact the most recent version. I have now spent

Hello all, I'm trying to write a function using the gridBase package. I'd like to push several base subplots to a larger plot constructed with grid. However, I'm having trouble getting consistent results when running the function when the plotting window (quartz) is closed, when it is left open and the plot function is repeated to the same window, and when the output is saved to a

Dear All, There might be a bug in heatmap(): x <- matrix(rnorm(10*5), 10, 5) heatmap(x, labCol=1:5) Error in axis(1, 1L:nc, labels = labCol, las = 2, line = -0.5, tick = 0, : 'at' and 'labels' lengths differ, 5 != 10 Works fine under 4.0.1. Looking at the code from 4.0.1 vs. Rdevel, esp. this part sticks out: labRow <- labRow[rowInd] %||% rownames(x) %||%

I am setting up base plots -- one in viewport A and and one in B. This part works fine. But if I go back to A after having done B and add horizontal lines it seems to not use the correct coordinates. How do I tell it to resume using A's coordinates? I am already using par(fig = gridFIG()) but it seems that that's not enough to reestablish them. What happens is that when I go back to

I am setting up base plots -- one in viewport A and and one in B. This part works fine. But if I go back to A after having done B and add horizontal lines it seems to not use the correct coordinates. How do I tell it to resume using A's coordinates? I am already using par(fig = gridFIG()) but it seems that that's not enough to reestablish them. What happens is that when I go back to

Hi, I have searched for a solution but I failed to find an answer. I am hoping you may be able to help me. I have a data set where I have observations for a number of units (n =~40) over a period of time (t =~100) and I have a variable (Z) that codes a categorical variable for each observation. I want to produce a 2D plot where time is on the x-axis and units are on the y-axis. Then each block

Dear All, I ma having a trouble in generating a figure containing 3 insets with the gridBase package. I always get an error message of the kind: Error in gridPLT() : Figure region too small and/or viewport too large No matter which parameters I choose. The plots works nicely with two insets only, but when I try adding the third one, my troubles begin. I am probably doing something wrong in the

Dear All: I tried to use heatmap.2 to generate hierarchical clustering using the following command: heatmap.2(datamatrix, scale="row", trace="none", col=greenred(256), labRow=genelist[,1], margins=c(10,10), Rowv=TRUE, Colv=TRUE) datamatrix is subset of a RMA normalized data subset by a genelist. The problem is a lot of times, the z-score in key are from, like -5 to 15 or

Dear R users, I have a question regarding how to make row labels readable in a heat map. I have successfully made a heat map using function "heatmap.2" in the package "gplots". However, as there are many rows in the heat map, I have difficulties labeling them (heatmap.2 provides a parameter "labRow" to label the row names, but as I have too many rows, I can not make

I'm trying to display the probability space of a function and wanted to see specifically where the maximum posterior probability is. data=read.table("PosteriorData",header=F) mat=as.matrix(data) heatmap.2(mat/max(mat,na.rm=T),dendrogram=c("none"),trace=c("none"), Rowv=F,Colv=F,labRow=3*c(10:-10),labCol=3*c(-10:10),symm=T,col=rainbow(100)) As the