[BioC] limma question
James W. MacDonald
jmacdon at med.umich.edu
Fri Feb 24 01:47:46 CET 2006
In my experience lmFit() never randomly adds anything including an
intercept when you have directly specified a design matrix.
Jim
Naomi Altman wrote:
> Isn't the problem here just that lmFit adds an intercept so that the X
> matrix is singular?
>
> --Naomi
>
> At 03:14 PM 2/23/2006, James W. MacDonald wrote:
>
>> 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
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>
>
> Naomi S. Altman 814-865-3791 (voice)
> Associate Professor
> Dept. of Statistics 814-863-7114 (fax)
> Penn State University 814-865-1348 (Statistics)
> University Park, PA 16802-2111
>
--
James W. MacDonald
University of Michigan
Affymetrix and cDNA Microarray Core
1500 E Medical Center Drive
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