similar to: R-1.8.0 memory.limit()

I've tried looking at ?colors and ?palette and if I'm understanding it correctly, I'm supposed to type in (for example) palette(rainbow(13)) before I type in my plot (of 13 lines) if I want 13 different colors. But this does not work. Other things that i have tried besides rainbow give me errors. Am I just doing something completely wrong? Anna

I get a segfault when using savePlot("foo","pdf"). FWIW, I don't have a problem using the canonical pdf("foo.pdf") plot(something) dev.off() Best, Jim > R.Version() $platform [1] "i386-pc-mingw32" $arch [1] "i386" $os [1] "mingw32" $system [1] "i386, mingw32" $status [1] "Under development (unstable)"

> R.version.string [1] "R version 2.1.0, 2005-03-17" I don't see anything in either https://svn.r-project.org/R/trunk/NEWS or in the Changes file for R-2.1.0 about changes in package.skeleton() (nor in the help page), but when I run this function, all the .Rd files produced are of the data format even if all I have in my .GlobalEnv are functions. A trivial example is to run the

Hi all, In the CHANGES file for R-2.3.0alpha, there is the following statement: winMenuAdd() now has no limits on the number of menus or items, and names are now limited to 500 (not 50) bytes. However, I can reproducibly get a segfault using this (admittedly silly) example: for( i in 1:5) winMenuAdd(paste("Test", letters[i], sep="")) for(i in 1:5) for(j in 1:24)

The function knn.cv in the class package doesn't have error checking to ensure that the length of the classlabel argument is equal to the number of rows in the test set. If the classlabel is short, the result is often a segfault. > library(class) > dat <- matrix(rnorm(1000), nrow=10) > cl <- c(rep(1,5), rep(2,5)) > cl2 <- c(rep(1,5), rep(2,4)) > knn.cv(dat, cl) [1] 2

Hi, I am unable to build a package I maintain using a relatively current build of R-2.4.0 alpha, whereas the package builds just fine on R-2.3.1. Both versions of R were built from source. I'm hoping a guRu might be able to give some help. Some snippets from the build process: R-2.3.1 making DLL ... gcc -Ic:/R-2.3.1/src/extra/zlib -DHAVE_ZLIB -Ic:/R-2.3.1/include -Wall -O2 -c

If one wants to unload compiled code for a package containing a namespace, my understanding is that .onUnload() should be used, with a call to library.dynam.unload(). This is used in e.g., the stats and methods packages, but it appears to me that the compiled code is not being unloaded when the package is detached(). Am I misunderstanding something? Best, Jim > search() [1]

I use dChip and Affymetrix gtype to apply some genotyping functions to some microarrays data. Do you know if similar R functions exist? Thank you! Francesco ----------------------------------------------------------- Francesco Falciano, Ph.D. CINECA (High Performance Systems) via Magnanelli, 6/3 40033 Casalecchio di Reno (BO)-ITALY tel: +39-051-6171724 fax: +39-051-6132198 e-mail:

I am running R 2.0.0, GOstats 1.1.1 and GO 1.7.0, and when I use the function GOHyperG, I have the following error: w1<-as.list(hgu95av2LOCUSID) w2<-unique(unlist(w1)) set.seed(123) myLL<-sample(w2,100) xx <- GOHyperG(myLL) Error in mget(x, env = GOTERM, ifnotfound = NA) : recursive default argument reference In fact first I tried this function with my locusId ' list (with

I am trying to work with 75 affymetrix U133plus2 chips and am running into memory allocation errors when trying to merge or convert probe level data to expression values. I keep getting - Error: cannot allocate vector of size 561011 Kb and that is simply with a data subset. Is there a way around this limitation? -- Rodney A. Staggs Cancer Center Informatics Shared Resource 425 Delaware St S.E.

Hi, I have a package that contains two vignettes that both use saved objects in the examples directory of the package. With previous versions of R I could have a code chunk in the vignette like this: <<echo=false>>= load("../examples/somedata.Rdata") @ followed by a code chunk like <<eval=false>>= foo <- bar("data") @ that simulated the actual

Hi everyone, My purpose is to read a .CEL file into R. The .CEL file was created from a .CAB by using DTT software found on Affymetrix website I read the .CEL file in R using ReadAffy as follows: > d2=ReadAffy(widget=T) and I complete the fields as required. It does not complain. For example I could find the description: > description(d2) Experimenter name: BB Laboratory: FFL Contact

Is there any way to control the format of the axis label numbers on a plot? More specifically, I have some plots that get axes with label numbers in exponential format, and I'd like to change that to non-exponential. Thanks!! -- M. Edward (Ed) Borasky, MS, MNLP, NST, FBG, PGS & PTA znmeb at borasky-research.net http://www.borasky-research.net

Hello, I have a problem, I can't install the package 'mgu74av2cdf'. I downloaded the zip file, yet when asked the R console to install it from a zip file, I got the answer: "Error in file(file, "r") : unable to open connection In addition: Warning messages: 1: error -1 in extracting from zip file 2: cannot open file `mgu74av2cdf/DESCRIPTION'

Hi, I have a data matrix of gene expression data from two groups that I would like to compare using the t-test. The data has been processed using RMA and transformed using log2. I would like to compare the two groups for each gene (N=10,000 genes) and have a result that lists the p-value for each test. Can anyone help me with the t.test setup for this analysis? Thanks, Rich Rich Roth, PhD

Hi everyone, I'm trying to generate an RNA degradation plot of the Estrogen example data plot, but seem to get an error. I've tried defining an ylim value, ylim=c(0,30) , but it doesn't seem to work either. My code is as follows: > RNAdeg<-AffyRNAdeg(Data) > png(DegLoc, width=720, height=720) > par(ann=FALSE) > par(mar=c(3,3,0.1,0.1)) >

Hello does anyone know what algorithm is used to produce the hierarchical clustering in the gplots package using the function heatmap.2? I think it may be the complete linkage clustering algorithm, but I can't find a source that seems reliable. Thank you and sorry if I posted this in the wrong place. If I have, please let me know and I will move it to the appropriate list. -- View this

Can someone throw light on the following problem I am having with RODBC? There's an Excel file I am trying to read from, it has one sheet named 'nameclass'. Thanks in anticipation. Vishal Belsare > library(RODBC) > con = odbcConnectExcel(file.choose()) > tbls <- sqlTables(con) > tbls TABLE_CAT TABLE_SCHEM TABLE_NAME TABLE_TYPE REMARKS 1

All, Anybody can tell how to export a R plot onto Word (or Power Point)? Many thanks in advance. -MY

Hello, I've a question about the fdr method in p.adjust: What is the threshold of the FDR, and is it possible to change this threshold? As I understand the FDR (please correct) it adjusts the p-values so that for less than N% (say the cutoff is 25%) of the alternative hypothesis the Null is in fact true. thanks a lot for help, +regards, Arne