similar to: t.test question

Hi, I am trying to install the BioC package from bioconductor onto a Windows Server 2003 machine. I can connect to bioconductor and when I run the getBioC("affy","release") function it starts to download but then it stops with the following error: Error: unable to create temporary directory 'C:\Program Files\R\R-2.2.1\release\file1eb26e9' The download just ends

Using R-1.8.0 (d/l and compiled on 2003-10-01) on WinXP, I seem to be unable to determine the maximum memory allocated to R. The help still says to use memory.limit(size=NA), but this returns the value NA. In addition, I have set --max-mem-size=2G but I run out of memory somewhere around 500Mb (which is why I am trying to find out how much memory is allocated). I don't have any other programs

I have a very basic question about limma. Assume I have experiments from 3 or more RNA sources in a reference design. It is easy to define individual contrasts but I want to specify a contrast matrix that tests for significant differences among ALL the different RNA sources (i.e. the analogous thing to a simple One-Way ANOVA). How can I do that? Thanks! Max --

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

Hello, I am using the following code to plot an MVA plot. library(affy) library(Biobase) library(limma) library(gcrma) pd<-read.phenoData("Clk.targets.2.txt",header=TRUE, row.names=1,as.is=TRUE,sep="\t") Data <- ReadAffy(filenames=pData(pd)$FileName,phenoData=pd) Print(Data) eset <- gcrma(Data) write.exprs(eset,

Hi together, we received a question about the bioconductor package, maybe anybody could help the guy and I will deliver the message over our forum to him. If this way is allowed ... Regards Knut Here the question: Hello, I have installed R on a linux machine. I have then installed the biocLite.R package from bioconductor. I wanted to install a few other packages - hgu133plus2cdf, rma . I

Dear R-Help, I am using the function "topTable" from the LIMMA package. To estimate adjusted P-values there are several options (adjust="fdr" , adjust="BH") as shown below: topTable(fit, number = 10, adjust = "BH", fit$Name) I guess any of these options (fdr, BH, etc.) is using a default of FDR=0.05 which is quite conservative (i.e., very

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)"

How can I get a single backslash in a character string? My goal is to escape dots in a string that will be used as a regular expression. I thought I could do it this way: gsub(".", "\\.", x) Unfortunately, "\\" does not represent a literal backslash as I expected, but rather a pair of backslashes: > "\\." [1] "\\." > "\\" [1]

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 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

> 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, 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

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)) >

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

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