[BioC] ExpressionSet, fitPLM confusion

James W. MacDonald jmacdon at med.umich.edu
Thu Dec 14 15:42:00 CET 2006


Hi Georg,

In addition to what Wolfgang told you, the affyPLM package in the devel 
repository won't work with R-2.5.0 yet, because some of the changes that 
have been made to the ExpressionSet object have not been addressed in 
the affyPLM package.

 > dat <- ReadAffy()
 > pset <- fitPLM(dat)
Error in checkSlotAssignment(object, name, value) :
	assignment of an object of class "AnnotatedDataFrame" is not valid for 
slot "phenoData" in an object of class "PLMset"; is(value, "phenoData") 
is not TRUE
In addition: Warning message:
The exprSet class is deprecated, use ExpressionSet instead

Your best bet at this time is to go back to using release versions.

Best,

Jim


Wolfgang Huber wrote:
> 
> Dear Georg,
> 
> to use the devel versions of Bioconductor (ie what will be Bioconductor
> 2.0) you will need to use a more recent version of R (what will be R
> 2.5): the new code in Biobase uses an argument "msg" of the .Deprecated
> function that was not yet there in R 2.4.
> 
> See http://www.bioconductor.org/download/BioC2.0/
> "BioC 2.0 is currently in development. It is being designed to work with
> R 2.5.z and is expected to be released in April of 2007."
> 
>  Best wishes
>  Wolfgang
> 
>  Otto wrote:
> 
>>Dear Bioconductors,
>>
>>Once more, I am a bit confused about an error linked with the new
>>ExpressionSet class. What I am trying to do is to apply fitPLM to an
>>AffyBatch object. First, I generate the AffyBatch object like this:
>>
>>targets<-new("AnnotatedDataFrame")
>>pData(targets)<-read.table("Targets.txt", header=TRUE, row.names=2)
>>
>>
>>Data<-ReadAffy(celfile.path =cel.path, 
>>               filenames = sampleNames(targets),
>>               verbose=TRUE)
>>
>>
>>Then I use fitPLM and get an error:
>>
>>Pset<-fitPLM(Data)
>>
>>Error in .Deprecated(msg = EXPRSET_DEPR_MSG) : 
>>	unused argument(s) (msg = "The exprSet class is deprecated, use ExpressionSet instead")
>>
>>However ?fitPLM tells me:
>>
>>This function converts an 'AffyBatch' into an 'PLMset' by fitting a
>>specified robust linear model to the probe level data
>>
>>And class() confirms that my object is indeed an AffyBatch:
>>
>>class(Data)
>>[1] "AffyBatch"
>>attr(,"package")
>>[1] "affy"
>>
>>
>>I think I am working with the last versions of the needed libraries
>>(see below). They were updated from the devel repository from
>>source. Any hints?
>>
>>Best,
>>
>>Georg
>>
>>
>>sessionInfo()
>>
>>R version 2.4.0 (2006-10-03) 
>>x86_64-redhat-linux-gnu 
>>
>>locale:
>>LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=en_US.UTF-8;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C
>>
>>attached base packages:
>>[1] "splines"   "tools"     "stats"     "graphics"  "grDevices" "utils"    
>>[7] "datasets"  "methods"   "base"     
>>
>>other attached packages:
>>zebrafishcdf      affyPLM        gcrma  matchprobes     affydata         affy 
>>    "1.15.0"     "1.11.8"      "2.7.1"      "1.7.1"     "1.11.1"    "1.13.12" 
>>      affyio      Biobase 
>>     "1.3.1"    "1.13.28"
>>
>>_______________________________________________
>>Bioconductor mailing list
>>Bioconductor at stat.math.ethz.ch
>>https://stat.ethz.ch/mailman/listinfo/bioconductor
>>Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
> 
> 
> 


-- 
James W. MacDonald, M.S.
Biostatistician
Affymetrix and cDNA Microarray Core
University of Michigan Cancer Center
1500 E. Medical Center Drive
7410 CCGC
Ann Arbor MI 48109
734-647-5623


**********************************************************
Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues.



More information about the Bioconductor mailing list