similar to: problem with creation of eSet

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

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

Full_Name: Colin A. Smith Version: 1.8.0 OS: Mac OS X 10.2.6 Submission from: (NULL) (216.102.90.18) Both Biobase and my package annaffy use S4 classes to define methods for "[". Both packages use the save image method of installation. (See annaffy 1.0.3 in BioC CVS.) Depending on how both packages are loaded, the Biobase definitions seem to be getting masked out: >

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

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

-------------- next part -------------- An embedded message was scrubbed... From: Barbara Cegielska <barceg at ibch.poznan.pl> Subject: loading SPOT file Date: Wed, 06 May 2009 16:48:35 +0200 Size: 1151 URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20090506/e35e3f7e/attachment-0002.eml>

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

Hi i'm a french medical student, i have some data that i import from excel. My colomn of the datafram are the localisations of metastasis. If there is a metatsasis there is the symbol "_". i want to exclude the row without metastasis wich represent the NA data. so, i wrote this mela is the data fram mela1=ifelse(mela[,c(11:12,14:21,23,24)]=="_",1,0) # selection of the

Dear R-helpers & bioconductor Sorry for cross-posting, this concerns R-programming stuff applied on Bioconductor context. Also sorry for this long message, I try to be complete in my request. I am trying to write a subset method for a specific class (ExpressionSet from Bioconductor) allowing selection more flexible than "[" method . The schema I am thinking for is the following:

Hello, I want to have the x-axis title of my plot in 2 lines, centered: experiment 1: log2(Ratio H/L) I know that in principle that works with '\n'. However, I am also using the 'substitute' command for my axis title. However, it does not make a new line. What I have so far: logbase <- 2 test <- "bait" cellline <- "cellline" plot(

Hi, I have a data.frame that looks something like this. feature 5637 ? ? ? ? ? ? ? ?1321N1 feature1 -0.568750616 -0.934748758 feature2 -0.913080902 -0.941455172 feature3 0.442477294 -0.257921866 I want to change it to look like this. feature ? ? ?cell.line ? ? value feature1 5637 -0.568750616 feature2 5637 -0.913080902 feature3 5637 0.442477294 feature1 1321N1

Hello, I'm wondering how I can merge two featuresets into one. My dataset is two sets of microarray data and it looks like followings: > rawData $v1 TilingFeatureSet (storageMode: lockedEnvironment) assayData: 2197815 features, 59 samples element names: channel1, channel2 protocolData rowNames: LT290677RU_D1_2011-02-16 LT286300LU_D1_2010-07-24 ... LT003990RU_D1_2010-11-04 (59

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

Hello, I will appreciate any suggestion for this simple problem. I have a matrix of characters "A", "P", "M". I want to convert those characters to integers 0,1, 1 respectively. I am using the following R statements: exprs.new<-gsub("P",1,exprs) exprs.new1 <- gsub("A",0,exprs.new) exprs.new2 <-

Full_Name: Manisha Desai Version: 1.3.0 OS: UNIX Submission from: (NULL) (156.111.36.213) There may be a bug in the table command. The value of 0 does not appear to show up (at least in this version of R). The following is an example: > blah<-c(rep(1,100),rep(0,100)) > table(blah) Gives the following: blah 1 100 Doing a summary shows that blah does contain values of 0. >

Hi, Thanks in advance for any advice you can give me, I am very stumped on this problem... I use R every day and consider myself a confident user, but this seems to be an elementary problem.. Outline of problem: I am analysing the results of a study on protein expression in cancer tissues. I have raw intensities from 2 different types of cancer and normal tissue, which can be taken from several

I'm seeing weird issues in methods initialization, i.e. loading marrayClasses loads Biobase, and when explicitly done, as in library(Biobase) library(marrayClasses) is fine, but when Biobase is loaded via a require statement in marrayClasses' .First.Lib, I end up with: Warning message: In the method signature for function "coerce", class "exprSet" has no

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