[BioC] limma question

James W. MacDonald jmacdon at med.umich.edu
Thu Feb 23 21:14:28 CET 2006


Hi Sun,

Sun, Yezhou wrote:
> Hi, there,
> 
>  
> 
> I am using limma for some single channel microarray data. Date format
> follows that suggested by user guide: probes in row and arrays in
> column. Data set has no missing value and was read into R by read.table
> function. Design matrix is as following:
> 
>  
> 
>             P          R
> 
> P1        1          0
> 
> P2        1          0
> 
> P3        1          0
> 
> P4        1          0
> 
> P5        1          0
> 
> P6        1          0
> 
> P7        1          0
> 
> R1        0          1
> 
> R2        0          1
> 
> R3        0          1
> 
> R4        0          1
> 
> R5        0          1
> 
> R6        0          1
> 
>  
> 
> When I ran lmFit(data, design), I got this error: 
> 
>  
> 
> 
>>fit<-lmFit(data, design)

It might help if you also showed what is in the data.frame (note that 
'data' isn't a good name - you are masking an existing function). Show 
us what you get from

data[1,]

There may be extra things in there that you don't intend.

Best,

Jim


> 
> 
> Error in qr(x) : NA/NaN/Inf in foreign function call (arg 1)
> 
> In addition: Warning message:
> 
> NAs introduced by coercion
> 
>  
> 
> Traceback gave this: 
> 
> 
>>traceback()
> 
> 
>  
> 
> 3: qr(x)
> 
> 2: nonEstimable(design)
> 
> 1: lmFit(data, design)
> 
>  
> 
> I have converted data and design to numeric matrices but they never
> worked. Could you please what's the problem? I really don't know what's
> going on with it. Thanks a lot!
> 
>  
> 
> Yezhou Sun
> 
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor


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



More information about the Bioconductor mailing list