similar to: custom subset method / handling columns selection as logic in '...' parameter

Hello Everyone, I am writing programs in R from 7 months and I am able to solve most of the errors/issues except for this current post. My Task is to read a Microsoft Excel file(textE_to_affy.csv) which contains the Microarray Expression Values collected from the Illumina Microarray experiment. These collected intensity values need to be normalized(Rank Invariant Normalization) by using the R

Hi, I am having troubles with creating an eSet and would appreciate any help on the following problem. I am trying to create an eSet using the following code pd <- read.table(file="pdata.txt",header =TRUE,row.names=1); colnames(pd) <- c("type","tumor","time","id"); pdN <- list(type =

Hi, Can anyone help me of the affylmGUI package, I can't get it work and searched for google but can't find any proper solutions. I get the error information each time when I load my cells files, which are shown in the following links. ------------- http://clarezoe.googlepages.com/1.png http://clarezoe.googlepages.com/2.png http://clarezoe.googlepages.com/3.png ------------- Errors also

Hi all, I have this microarray large microarray data set (ALL) from which I would like to subset or extract a set of data based on a factor ($mol.biol). I looked up some example of subsetting in, picked up two commands and tried both but I got error messages as follows > testset <- subset(ALL, ALL$mol.biol %in% c("BCR/ABL","ALL1/AF4")) >> Error in

Hi all, What are current methods people use in R to identify mis-spelled column names when selecting columns from a data frame? Alice Johnson recently tackled this issue (see [BioC] posting below). Due to a mis-spelled column name ("FileName" instead of "Filename") which produced no warning, Alice spent a fair amount of time tracking down this bug. With my fumbling fingers

Full_Name: Guy W. Tillinghast Version: 2.8.0 OS: Windows XP professional Submission from: (NULL) (24.248.24.3) I am encountering difficulty asking R to manipulate the correct columns in Expression Set class (object 4). I download the ALL data with: library(golubEsets) data(Golub_Merge) Note, the data has the samples not in order. This is not R's fault (at least not that I can tell): >

Hi, Reading help("Documentation"), I'm led to believe that a help call like: ?myFun(x, sqrt(wt)) Will search for help on the appropriate method in the case that myFun is generic. This isn't working for me. Here is an example using the Biobase package: ## If Biobase is not installed source("http://bioconductor.org/biocLite.R") biocLite("Biobase")

On 10/11/2006 2:48 PM, Seth Falcon wrote: > Hi, > > Reading help("Documentation"), I'm led to believe that a help call > like: > > ?myFun(x, sqrt(wt)) > > Will search for help on the appropriate method in the case that myFun > is generic. This isn't working for me. Here is an example using the > Biobase package: > > ## If Biobase is

Hi, I want to change .RDA file to a text file. So I did as follows. >load("my.rda") >ls() ---> then it showed [1] exprs >write.table(exprs,"C:\\my.txt",sep="\t") I was successful with the first .RDA file. Then I used the same commands with another .RDA file (172 MB)which is 4 times bigger than the first file (41.2 MB). When I put the last command

Hi, I'm very confused about R structures and the methods to go with them. I'm using R for microarray analysis with Bioconductors. Suppose without reading the documentations, what's the best way to explore a data structure when you know nothing about it? I am currently using is() / class() to see what the object is. str() / attributes() to probe inside the object, and

Hello, I am to run this R script but i keep getting this error. > expr<-exprs(golubMerge) Warning message: The exprSet class is deprecated, use ExpressionSet instead I tried to find information on the website but no luck. (exprSet...etc) thank you. -- View this message in context: http://www.nabble.com/expression-matrix-tp15912874p15912874.html Sent from the R help mailing list archive

Running R2.4.0 on Apple Mac OS X 10.4.8, in Emacs ESS mode, and also R.app. In an attempt to learn a bit more about a particular method (geneNames in package affy) I invoked getMethods("geneNames") which produced geneNames methods, but not the one in affy (output below). I had to know the signature (AffyBatch) in order to find the method > getMethod("geneNames",

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,

Hi, I started this post on bioc-devel but this seems to be more general: https://stat.ethz.ch/pipermail/bioc-devel/2013-May/004311.html See reproducible example from Martin below. Thank you. Renaud ---------- Forwarded message ---------- From: Martin Morgan <mtmorgan at fhcrc.org> Date: 7 May 2013 19:55 Subject: Re: [Bioc-devel] ExpressionSet and LumiBatch: inheritance problem in S4

Hi, I m not able to plot normalized data(normalization by rma) using boxplot. I don't know why? basically, object(formed of normalized data) belong to ExpressionSet class. It is showing error Error in x[!xna] : object of type 'S4' is not subsettable In addition: Warning messages: 1: In is.na(x) : is.na() applied to non-(list or vector) of type 'S4' 2: In is.na(x) : is.na()

Hi all,   I have loaded the LIMMA and Biobase package and tried these commands:   library(limma) library("Biobase") data <- read.table("c:/temp/data.txt",header=T,row.names=1) ExpressionData <- as.matrix(data[,c(2,3,4,6,7,8)]) eset <- new("ExpressionSet", exprs = ExpressionData) design <- cbind(WT=1,P=c(0,1,1,0,1,1),G=c(0,1,0,0,1,0)) fit <-

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

dear all, I have a problem with a masked object in a package we created here. we make a package for a workflow of internal analysis of microarray data. to create the package we used: > install.packages(pkgs="affyAnalysis", repos=NULL) > R CMD INSTALL affyAnalysis Erzeuge Verzeichnisse ... Erzeuge DESCRIPTION ... Erzeuge NAMESPACE ... Erzeuge Read-and-delete-me ... Kopiere

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 have the following code, from which I get the following error message: Error in dim(data) <- dim : attempt to set an attribute on NULL I think the error is coming from the part of my code in BOLD RED. The script works fine until then. #Load libraries source("http://bioconductor.org/biocLite.R") biocLite() library(limma) library(Biobase) #change directory to folder where