[BioC] Limma - Result output - Error - Help

Gordon K Smyth smyth at wehi.EDU.AU
Sat Apr 6 01:08:51 CEST 2013


Dear Atul,

You have uses ebayes() when you intended eBayes().

Best wishes
Gordon

> Date: Thu, 04 Apr 2013 19:01:24 -0400
> From: Atul Kakrana <atulkakrana at gmail.com>
> To: "bioconductor at r-project.org" <bioconductor at r-project.org>
> Subject: [BioC] Limma - Result output - Error - Help
> Message-ID: <515E0644.9050304 at gmail.com>
> Content-Type: text/plain; charset=ISO-8859-1
>
> Hello All,
>
> I wrote a script to analyse affymetrix data but its giving an error
> while using the topTable from limma:
>
> Error in array(x, c(length(x), 1L), if (!is.null(names(x)))
> list(names(x),  :
>  'data' must be of a vector type
>
> I checked all the steps that might be causing the error but couldn't
> find out. So, I trimmed my script to minimal basic and it still gives me
> same error.
>
>> mydata <- ReadAffy()
>> esetRMA <- rma(mydata)
> Background correcting
> Normalizing
> Calculating Expression
>> esetRMA
> ExpressionSet (storageMode: lockedEnvironment)
> assayData: 45101 features, 12 samples
>  element names: exprs
> protocolData
>  sampleNames: 00.0A.CEL 00.0B.CEL ... 12.5C.CEL (12 total)
>  varLabels: ScanDate
>  varMetadata: labelDescription
> phenoData
>  sampleNames: 00.0A.CEL 00.0B.CEL ... 12.5C.CEL (12 total)
>  varLabels: sample
>  varMetadata: labelDescription
> featureData: none
> experimentData: use 'experimentData(object)'
> Annotation: mouse4302
>> sampletype <- c('1','1','1','2','2','2','3','3','3','4','4','4')
>> group <-factor(sampletype)
>> design<- model.matrix(~0+group)
>> colnames(design) <- c('WB','Ten','Eleven','Twelve')
>> fit <- lmFit(esetRMA, design)
>> contrast.matrix<-makeContrasts(Ten-WB, levels=design)
>> fit2<- contrasts.fit(fit,contrast.matrix)
>> fit2 <- ebayes(fit2)
>
>> topTable(fit2, adjust = "fdr", number = 50)
> Error in 1:ncol(fit) : argument of length 0
>
>> topTable(fit2,coef = 1, adjust = "fdr", number = 50)
> Error in array(x, c(length(x), 1L), if (!is.null(names(x)))
> list(names(x),  :
>  'data' must be of a vector type
>
>> results <- decideTests(fit2)
> Error in decideTests(fit2) : Need MArrayLM object
>
>
> Could anybody please explain me or point out what I am missing here? Its
> strange that I wrote a similar script for Illumina which works just fine.
>
> AK
>
> -- 
> Atul Kakrana
> Delaware Technology Park
>

______________________________________________________________________
The information in this email is confidential and intend...{{dropped:4}}



More information about the Bioconductor mailing list