similar to: How to avoid losing a sample as reference running an individual differential expression analysis (LIMMA)

Dear Sirs, How can I revert to an older limma version? Typing "install.packages("limma")" in R gives a list of mirrors. How can I install the version I want after I obtain and untar the file (e.g, limma_2.9.1.tar.gz)? I am running R 2.5.0 on a Linux machine (CentOS 5). When using limma it will not go past the read.maimages command. I get this error: Error in

Dear All, How to cite the PDF user's guide for the LIMMA package? This is not about how to cite the LIMMA package. Roger Roger L. Vallejo, Ph.D. Computational Biologist & Geneticist U.S. Department of Agriculture, ARS National Center for Cool & Cold Water Aquaculture 11861 Leetown Road Kearneysville, WV 25430 Voice: (304) 724-8340 Ext. 2141 Email: roger.vallejo@ars.usda.gov

Hi all. I've got a really dumb question for anyone. How do I write the output of a limma analysis (basically the topTable) to a text file? I want to output the topTable for the entire microarray (not really a topTable anymore I suppose..). Thanks for any advice! -Josh

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

Dear R-Help, I am using the LIMMA User's Guide 5 January 2007 PDF version. For the example show in Section 7.4 DIRECT TWO-COLOR DESIGNS (pgs. 33-34), I could not grasp the rationale in developing the contrast.matrix with these R statements (">" indicates the R command prompt): > contrast.matrix <-

Hi, I have a question about limma. I have data from spotted arrays. I have class A and Class B on the same slide. In Limma if I put class A as the reference (ref), are the list of genes I get in the output for class A or class B? #design <- modelMatrix(targets, ref="classA") design Josephine Brennan [[alternative HTML version deleted]]

Dear All, I am using the LIMMA package to create 2 contrasts for my data and then calculating the vennCounts of the decideTests from the contrast.fit to be able to create venn Diagrams. The code works fine but the summary(results) shows zeros for all i.e. no gene were up regulated or downregulated. This is not true for my data since toptable output shows Log fold change greater than > 2. I am

Hello Limma users A quick question, I hope: I have dual-channel spotted expression arrays in a simple loop design (no dye swaps), viz: 1 vs reference 2 vs 1 3 vs 2 reference vs 3 There are 4 replicate spots for each probe on each array. It seems as if getting meaningful spot summary results using duplicateCorrelation() is not possible in my case as the spots are not regularly spaced, and I

Hi, I am following the protocol outlined here for analysis of single channel Agilent microarray data: http://matticklab.com/index.php?title=Single_channel_analysis_of_Agilent_microarray_data_with_Limma I keep getting the following error message when using Limma's read.maimages function to load my data into an RGList object: Error in RG[[a]][, i] <- obj[, columns[[a]]] : number of

Hi, following up on my own question, I found smaller example that does not require LIMMA: setClass("FOOCLASS", representation("list") ) ma = new("FOOCLASS", list(M=matrix(rnorm(300), 30,10))) > ma * ma$M Error: C stack usage 7970512 is too close to the limit > library(xlsx) Loading required package: rJava Loading required package: xlsxjars

I am posting this to both R and BioC communities because I believe there is a lot of confusion on this topic in both communities (having searched the mail archives of both) and I am hoping that someone will have information that can be shared with both communities. I have seen countless questions on the BioC list regarding limma (Bioconductor) and its calculation of FDR. Some of them involved

I am posting this to both R and BioC communities because I believe there is a lot of confusion on this topic in both communities (having searched the mail archives of both) and I am hoping that someone will have information that can be shared with both communities. I have seen countless questions on the BioC list regarding limma (Bioconductor) and its calculation of FDR. Some of them involved

Hello, I have just completed an analysis of microarray data using the limma package. The analysis appears to have worked fine. However, when I use the topTable function to get the significant genes, I get the following error: > topTable(fit2,coef=5,adjust="fdr") Error in array(x, c(length(x), 1), if (!is.null(names(x))) list(names(x), : attempt to set an attribute on NULL

Hi, I need some help with the interpretation of B statistics generated by eBayes in the limma package. I want to compare gene expression in three groups of Affy samples. The expression set has been imported from GeneSpring using GSload.expBC and a linear model fitted to the data using a design based on the three groups (6, 5, and 5 samples in each group, respectively). I have then made 3

Hi, I have two data set of normalized Affymetrix CEL files, wild type vs Control type.(each set have further three replicates). > wild.fish AffyBatch object size of arrays=712x712 features (10 kb) cdf=Zebrafish (15617 affyids) number of samples=3 number of genes=15617 annotation=zebrafish notes= > Dicer.fish AffyBatch object size of arrays=712x712 features (10 kb) cdf=Zebrafish (15617

On Wed, December 22, 2004 12:11 am, r.ghezzo at staff.mcgill.ca said: > ----- Forwarded message from r.ghezzo at staff.mcgill.ca ----- > Date: Mon, 20 Dec 2004 15:45:11 -0500 > From: r.ghezzo at staff.mcgill.ca > Reply-To: r.ghezzo at staff.mcgill.ca > Subject: [R] problems with limma > To: r-help at stat.math.ethz.ch > > I try to send this message To Gordon

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

Hi Everyone, I have been trying to read some Arrayvision files( 2 channel cDNA) and am having some problem. My code is : setwd('C:/work/data/limma/ndd1'); files <- c('ndd1_1.txt','ndd1_2.txt','ndd1_3.txt'); RG=read.maimages(files,"arrayvision",sep="\t"); #Normalisation MA=normalizeWithinArrays(RG); #plotPrintTipLoess(MA); #Fit Linear

Hi all, I'm working on microarray and currently analyzing the microrarray data using limmaGUI. Loop design has been applied in this experiment. This is a 2X2 factorial experiment where there are control and treatment at 2 different time points, week 6 and 9. The experimental design is almost the same as the limmaGUI work example: Weaver Data set. I would like to look at the effect of

Dear mailing group, This is my first time here. Glad to have this resource! I am currently trying to load an Excel file into R (limma package loaded) using the source(*name of directory*) command, but it cannot open the file. I renamed the file as .R and .RData, to no avail. The Excel data contains one gene name per row and about 100 data points per gene (columns). I am only used to