[BioC] dupcor.series

Gordon Smyth smyth at wehi.edu.au
Wed Mar 17 00:52:49 MET 2004


Jason,

I am a bit dubious about the merits of between array normalization when 
there are so many missing values (nearly 6000 on one of the arrays), but it 
shouldn't in principle cause an error to the duplicateCorrelation function. 
So there may be a bug in the function. Would you be prepared to send me the 
data matrix so that I can try it out myself?

Gordon

At 08:38 PM 16/03/2004, Jason Skelton wrote:
>>>*Error in chol(dinfo) : the leading minor of order 2 is not positive
>>>definite.*
>>
>>I've never seen this error either. What version of R/limma are you using? 
>>What do you see if you type
>>
>>summary(nbaResults$M)?
>
>Hi Gordon
>
>Using limma version 1.5.1
>& R version 1.8.0
>
>If I look at the normalized within Arrays data (Which works with 
>dupcor.series)
>
> > summary(nwam0m2MA$M)
>B737.10.03.2004.m0m2.5um.0019m0m2 B737.10.03.2004.m2m0.5um.0046m0m2
>Min.   :-1.417e+01                Min.   :-1.638e+01
>1st Qu.:-1.224e+00                1st Qu.:-6.466e-01
>Median :-2.273e-03                Median : 8.657e-03
>Mean   :-2.393e-01                Mean   :-1.095e-02
>3rd Qu.: 1.105e+00                3rd Qu.: 6.254e-01
>Max.   : 7.943e+00                Max.   : 8.213e+00
>NA's   : 5.811e+03                NA's   : 1.685e+03
>
>The normalised between arrays data looks like this
>I can't see that much difference but the quality of the second array 
>(B737.10.03.2004.m2m0.5um.0046m0m2)
>wasn't that great in the first place.
>
> > summary(nbamom2MA$M)
>B737.10.03.2004.m0m2.5um.0019m0m2 B737.10.03.2004.m2m0.5um.0046m0m2
>Min.   :-1.048e+01                Min.   : -22.15456
>1st Qu.:-9.056e-01                1st Qu.:  -0.87431
>Median :-1.681e-03                Median :   0.01171
>Mean   :-1.770e-01                Mean   :  -0.01481
>3rd Qu.: 8.174e-01                3rd Qu.:   0.84567
>Max.   : 5.874e+00                Max.   :  11.10602
>NA's   : 5.811e+03                NA's   :1685.00000
>
>Cheers
>
>Jason



More information about the Bioconductor mailing list