similar to: standard errors for expression intensities

Dear R-help Hi, I'm Won. I try to do microarray normalization by R. I use justRMA function within affy package, got error about segment fault. I don't know why it happen. I attached error below. Please help me. Thank you. Cheers, Won ======================= OS : Redhat linux Cpu : intel xeon X5570 Memory : 26Gb & OS : Ubuntu Cpu : intel q6600 Memory : 8Gb

I am trying to preprocess a large dataset of affymetrix data. Creating an affybatch is not possible with the computer I am running it on, so I have used the justRMA command to run RMA. I have read the affy document describing the justRMA command and the help documentation but I am unclear as to whether this command uses median polish after normalization. I assume this is the case but would like

Dear Community, I am very new to the world of Linux and R and I have stumbled upon a problem that I cannot seem to resolve on my own. Here is the relevant background: I am working on a 64-bit Linux Fedora Core 6 OS. I using R version 2.5.1. I have 3.8 Gb of RAM and 1.9 Gb of swap. As I see it, there are no restraints on the amount of memory that R can use imposed by this particular OS build.

Hi all , I have one query. i have list of some .cel files. in my program i have to mention the path of these .cel files part of my program is, rna.data<-exprs(justRMA(filenames=file.names, celfile.path=*datadir*, sampleNames=sample.names, phenoData=pheno.data, cdfname=cleancdfname(hg18_Affymetrix U133A))) in the place of "datadir" i have to mention the character string of the

Hello. By way of background, I am running out of memory when attempting to normalize the data from 160 affymetrix microarrays using justRMA (from the affy package). This is despite making 6 gigabytes of swap space available on our sgi irix machine (which has 2 gigabytes of ram). I have seen in various discussions statements such as "you will need at least 6 gigabytes of memory to normalize

I have CEL files from Affymetrix Mouse Array 430_2 and am trying to get the the individual PM intensities (11 per gene) for each sample. I would like to write out this into a tab delimited text file. Where am I stalling? This is what I've done: Change dir(to where CEL files are saved) Data <- ReadAffy() eset <- rma(Data) write.exprs(eset, file="mydata.txt") With this I am

Hi, I am a Ph.D. student from Québec, Canada. I’m a beginner with R and Bioconductor. Until now the only experience I have is in analyzing microarray data using affy and limma packages. Now I am trying to analyze Rat Gene 10 st arrays and I would like to run RMA analysis and Smyth moderated t test on those arrays. Since no cdf official package is available for those arrays, after reading many

The Bioconductor development team announces release 1.2 of the Bioconductor packages for the analysis of genomic data. Bioconductor is an open source bioinformatics software project based on the R language. Version 1.2 features: ===================== * All packages from the 1.1 release are included. All current bug fixes have been applied, and most have been upgraded and provide enhanced

Hello I am fitting a Cox PH model using the function coxph(). Does anyone know how to obtain the estimate of the covariance matrix under the null hypothesis. The function coxph.detail() does not seem to be useful for this purpose. Thanks, KD. [[alternative HTML version deleted]]

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

Dear R-mailing list Hope you can help me. I am using R for windows to analyze my 107 HGU133Plus2.0 chips, however, R chrash when I try to use ReadAffy(). I want to buy a computer that can handle all these arrays, do you know how big a computer I need to buy? Best, Skov, Denmark [[alternative HTML version deleted]]

Hi, I have built R (current development version) and BioConductor 1.7 with portland group compiler on a AMD Opteron. When I ran qc assessment on Affymetrix latin square data set, I got the following output, Loading required package: affy Loading required package: Biobase Loading required package: tools Welcome to Bioconductor Vignettes contain introductory material. To view,

Hey, I have code that can check the quality of a data set we're working with (expression data), and I'm having some trouble writing code that would make the expression data we have tie to other data we want to link it to (called phenotype data). Does anyone have any advice on how I could make a single object that would do this? Other relevant info: I want to use the pdata() function,

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,

I run R 1.9.0 on windows 2000, and have the following libraries installed: affydata_1.3.1 affy_1.4.23 Biobase_1.4.10 DynDoc_1.3.14 gcrma_1.0.6 hgu133acdf_1.4.3 hgu95av2cdf_1.4.3 hgu95av2probe_1.0 matchprobes_1.0.7 moe430acdf_1.4.3 multcomp_0.4-6 mvtnorm_0.6-6 rae230acdf_1.4.3 reposTools_1.3.29 rgu34acdf_1.4.3 tkWidgets_1.5.1 widgetTools_1.2.7 1. The rma function (in affy library) always crashes.

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.

Hello, I am using survreg for a simulation where data is randomly generated from a Weibull distribution and the model parameters estimated. When this is embedded within a loop and the process repeated for different datasets, the program aborts after a few runs in most cases giving the following message: Error in survreg.fit(X, Y, weights, offset, init = init, controlvals = control, :

Dear all, I was trying to normalize a subset of affy data those transcribe are either P or M (called pm_filter). I am able to normalize pm_filter subset by using RMA method, however MAS5 is not working. For RMA method, I used the following commend: est<-rma(affydata, subset=pm_filter). Could any help me; how do I do this by using MAS5 method? Regards, Anup Som -- Dr. Anup Som,

Full_Name: Edward J. Oakeley Version: 1.7.1 OS: Windows XP Submission from: (NULL) (212.47.183.3) This bug occurs when using the (D)COM server to connect to the "expresso" command of the Bioconductor Affy package. It may be a bug of R, (D)COM or Affy ut I will report it here anyway as it feels like an R bug. The Affy package when invoked will read a series of large (10Mb) text files

I am getting the following error in my script. I am very very new to R and have obtained this script from another person. #read file in (dummy data) starburst.plot<-function(affy.fold, affy.FDR)(ifelse( ((affy.fold) >=0), -1*log10(affy.FDR), 1*log10(affy.FDR))) starburst.plot<-function(meth.fold, meth.FDR)(ifelse( ((meth.fold) >=0), -1*log10(meth.FDR), 1*log10(affy.FDR))) At my next