Entering edit mode
Hi everyone,
I am using the genotype.Illumina function from crlmm. The following is my code and the error that I am getting:
> head(samplesheet)
SampleID SentrixPosition SentrixBarcode Array
TARGET-30-PAAPFA-10A 1906277330_A 1906277330 Illumina.com:HumanHap550v3
TARGET-30-PACRYY-10A 1906277899_A 1906277899 Illumina.com:HumanHap550v3
TARGET-30-PACUGP-10A 1870075248_A 1870075248 Illumina.com:HumanHap550v3
TARGET-30-PADENF-10A 1906277855_A 1906277855 Illumina.com:HumanHap550v3
TARGET-30-PADIEY-10A 1870147515_A 1870147515 Illumina.com:HumanHap550v3
TARGET-30-PADIRB-10A 1906277047_A 1906277047 Illumina.com:HumanHap550v3
> head(arrayNames)
[1] "1906277330_A"
[2] "1906277899_A"
[3] "1870075248_A"
[4] "1906277855_A"
[5] "1870147515_A"
[6] "1906277047_A"
> head(arrayInfo)
$barcode
[1] "SentrixBarcode"
$position
[1] "SentrixPosition"
> head(batch)
[1] "1" "1" "1" "1" "1" "1"
# run genotype.Illumina function
cnSet <- genotype.Illumina(sampleSheet=samplesheet,
arrayNames=arrayNames,
arrayInfoColNames=arrayInfo,
cdfName="human550v3b",
batch=batch)
Instantiate CNSet container.
path arg not set. Assuming files are in local directory, or that complete path is provided
Initializing container for genotyping and copy number estimation
Processing sample stratum 1 of 1
path arg not set. Assuming files are in local directory, or that complete path is provided
Loading chip annotation information.
Quantile normalizing 92 arrays by 10 strips.
|======================================================================| 100%
Loading snp annotation and mixture model parameters.
Calibrating 92 arrays.
|===================================================== | 76%Error in chol.default(crossprod(sweep(matS, 1, z[, 1], FUN = "*"), matS)) :
the leading minor of order 1 is not positive definite
In addition: Warning messages:
1: In readIdatFiles(sampleSheet = sampleSheet[sel, ], arrayNames = arrayNames[sel], :
Chips are not of the same type. Skipping 6110681250_R01C01_Grn.idat and 6110681250_R01C01_Red.idat
2: In readIdatFiles(sampleSheet = sampleSheet[sel, ], arrayNames = arrayNames[sel], :
Chips are not of the same type. Skipping 6184541022_R01C01_Grn.idat and 6184541022_R01C01_Red.idat
3: In readIdatFiles(sampleSheet = sampleSheet[sel, ], arrayNames = arrayNames[sel], :
Chips are not of the same type. Skipping 6110681250_R02C01_Grn.idat and 6110681250_R02C01_Red.idat
4: In readIdatFiles(sampleSheet = sampleSheet[sel, ], arrayNames = arrayNames[sel], :
Chips are not of the same type. Skipping 6110681250_R01C02_Grn.idat and 6110681250_R01C02_Red.idat
5: In readIdatFiles(sampleSheet = sampleSheet[sel, ], arrayNames = arrayNames[sel], :
Chips are not of the same type. Skipping 6110681250_R02C02_Grn.idat and 6110681250_R02C02_Red.idat
So it does say that some chiops where skipped because they are not of the same type but what does the following error mean?
Error in chol.default(crossprod(sweep(matS, 1, z[, 1], FUN = "*"), matS)) : the leading minor of order 1 is not positive definite
Any help would be much appreciated. Thanks!
