similar to: topTable function from LIMMA

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

Dear R Users, In the expresso function, which combination of these methods for data pre-processing (when using affymetrix oligo arrays) is the best: bgcorrect.metod = rma rma2 mas normalize.method = qspline quantiles loess pmcorrect.method = pmonly subtractmm mas summary.method = liwong avgdiff medianpolish mas There are many options within each method. I would appreciate a hint on the best

Hi Martin, In using the Cluster Package, I have results for PAM and DIANA clustering algorithms (below "part" and "hier" objects): part <- pam(trout, bestk) # PAM results hier <- diana(trout) # DIANA results GeneNames <- show(RG$genes) # Gene Names are in this object But

Dear R users, I am using the beadarray package. I am trying to upload raw bead-level data using these commands: ######################################################## library(beadarray) datadir <- ("C:/Computer_programs/R/beadarray/cecilia") targets = read.table("targets.txt", sep = "\t", header = TRUE, as.is = TRUE) BLData = readIllumina(arrayNames =NULL,

In using the NLME package (R 2.6.1 for Windows), I am having a problem in running an R script that used to run with no problems using a Linux OS in 2004. So I am wondering if during these last ~3 yrs we had major changes in the syntax of the NLME package that I am not aware. This is the R script: library(nlme) treat=as.factor(c(1,2,1,2,1,2,1,2)) mouse=as.factor(c(1,1,2,2,3,3,4,4))

I try to send this message To Gordon Smyth at smyth at vehi,edu.au but it bounced back, so here it is to r-help I am trying to use limma, just downloaded it from CRAN. I use R 2.0.1 on Win XP see the following: > library(RODBC) > chan1 <- odbcConnectExcel("D:/Data/mgc/Chips/Chips4.xls") > dd <- sqlFetch(chan1,"Raw") # all data 12000 > # > nzw <-

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

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

Hallo, > fit12<-lmFit(qrg[,1:2]) > t12<-toptable(fit12,adjust="fdr",number=25,genelist=qrg$genes[,1]) > t12 ID logFC t P.Value adj.P.Val B 522 PLAU_OP -6.836144 -8.420414 5.589416e-05 0.01212520 2.054965 1555 CD44_WIZ -6.569622 -8.227938 6.510169e-05 0.01212520 1.944046 Can anyone tell me what the difference is between P.Value

>>>>> "GS" == Gordon K Smyth <smyth@wehi.edu.au> >>>>> on Sat, 8 Jan 2005 01:11:30 +1100 (EST) writes: <.............> GS> p.adjust() unfortunately gives incorrect results when GS> 'p' includes NAs. The results from topTable are GS> correct. topTable() takes care to remove NAs before GS> passing

Hi all, My name is Amy, I am a masters student in Bioinformatics at North Carolina State University. I am working on a project and I am trying to use the lumi R package for microarray data analysis. I have shown the sample code here and have questions about modifying the sample code for my own data. lumi package in R, example.lumi, the sample data has 8000 features and 4 samples I have

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

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

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

You asked the same question on the Bioconductor mailing list back in August. At that time, you suggested yourself a solution for how the adjusted p-values should be interpreted. I answered your query and told you that your interpretation was correct. So I'm not sure what more can be said, except that you should read the article Wright (1992), which is cited in the help entry for p.adjust(),

Hi, I have a question concerning the analysis of some affymetrix chips. I downloaded some of the data from GEO GSE11324 (see below). In doing so I'm stuck after I identified the probesets with significant changes. I have problems in assigning probeset specific gene names as well as getting the genomic coordinates. Furthermore I have no clue how to deal with the fact, that most genes have

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 - When I export my files using the write.table command and name the file.xls and use the sep="/t" to try and separate the columns in excel, everything just stays in one column. how do i get them all into separate columns (it did this before but i can't figure out what's different between this time exporting and other times when i enter the same thing) example: >

I append below a suggested update for p.adjust(). 1. A new method "yh" for control of FDR is included which is valid for any dependency structure. Reference is Benjamini, Y., and Yekutieli, D. (2001). The control of the false discovery rate in multiple testing under dependency. Annals of Statistics 29, 1165-1188. 2. I've re-named the "fdr" method to "bh" but