[BioC] beadarraySNP and GenomeStudio report file
Christian Ruckert
cruckert at uni-muenster.de
Wed Apr 22 12:34:32 CEST 2009
I have followed the instructions in this thread
https://stat.ethz.ch/pipermail/bioconductor/2008-November/025282.html
apart from using GenomeStudio to read in Illumina Human Cyto SNP-12 Bead
Chip data.
But when trying to read in the files I get the following error:
mydata<-read.SnpSetIllumina(Sample_Map2Samplesheet("Sample_Map.txt"),reportfile="FinalReport.txt")
Error in read.SnpSetIllumina(Sample_Map2Samplesheet("Sample_Map.txt"), :
subscript out of bounds
The error occurs when running line:
G <- G[, selected]
content of variable "selected":
[1] "sentrix_id_4554318341_R01C01" "sentrix_id_4554318341_R01C02"
[3] "sentrix_id_4554318341_R02C01" "sentrix_id_4554318341_R02C02"
[5] "sentrix_id_4554318341_R03C01" "sentrix_id_4554318341_R03C02"
[7] "sentrix_id_4554318341_R04C01" "sentrix_id_4554318341_R04C02"
[9] "sentrix_id_4554318341_R05C01" "sentrix_id_4554318341_R05C02"
[11] "sentrix_id_4554318341_R06C01" "sentrix_id_4554318341_R06C02"
G[1:3,]
4554318341_R01C01 4554318341_R01C02 4554318341_R02C01
cnvi0111185 1416 917 1418
cnvi0111186 7080 6714 7680
cnvi0111187 982 773 1216
4554318341_R02C02 4554318341_R03C01 4554318341_R03C02
cnvi0111185 1166 1343 1200
cnvi0111186 7465 7158 8150
cnvi0111187 1016 1075 997
4554318341_R04C01 4554318341_R04C02 4554318341_R05C01
cnvi0111185 1445 1169 1485
cnvi0111186 8428 8732 8995
cnvi0111187 983 1005 1246
4554318341_R05C02 4554318341_R06C01 4554318341_R06C02
cnvi0111185 1819 1269 1449
cnvi0111186 8702 5399 4641
cnvi0111187 1317 1314 1183
I think I miss something in the input files:
Sample_map.txt (Group, Parent1, Parent2 and Replicate are empty):
Index Name ID Gender Plate Well Group Parent1 Parent2 Replicate
SentrixPosition
1 07-019117 4554318341_R01C01 Unknown WG9999912-DNA A01
4554318341_R01C01
2 08-001712 4554318341_R01C02 Unknown WG9999912-DNA D01
4554318341_R01C02
3 07-007836 4554318341_R02C01 Unknown WG9999912-DNA G01
4554318341_R02C01
4 06-013382 4554318341_R02C02 Unknown WG9999912-DNA B02
4554318341_R02C02
5 06-008925 4554318341_R03C01 Unknown WG9999912-DNA B01
4554318341_R03C01
6 06-006744 4554318341_R03C02 Unknown WG9999912-DNA E01
4554318341_R03C02
7 07-020190 4554318341_R04C01 Unknown WG9999912-DNA H01
4554318341_R04C01
8 NA10859_1 4554318341_R04C02 Unknown WG9999912-DNA C02
4554318341_R04C02
9 08-000591 4554318341_R05C01 Unknown WG9999912-DNA C01
4554318341_R05C01
10 07-015573 4554318341_R05C02 Unknown WG9999912-DNA F01
4554318341_R05C02
11 08-004776 4554318341_R06C01 Unknown WG9999912-DNA A02
4554318341_R06C01
12 NA10859_2 4554318341_R06C02 Unknown WG9999912-DNA D02
4554318341_R06C02
The output created by Sample_Map2Samplesheet:
[Header]
Investigator Name,
Project Name,
Experiment Name,
Date,
[Data]
Index,Sample_Name,ID,Gender,Sample_Plate,Well,Group,Parent1,Parent2,Replicate,Sentrix_Position,Pool_ID,Sentrix_ID
1,07-019117,4554318341_R01C01,Unknown,WG9999912-DNA,A01,,,,,4554318341_R01C01,pool_id,sentrix_id
2,08-001712,4554318341_R01C02,Unknown,WG9999912-DNA,D01,,,,,4554318341_R01C02,pool_id,sentrix_id
3,07-007836,4554318341_R02C01,Unknown,WG9999912-DNA,G01,,,,,4554318341_R02C01,pool_id,sentrix_id
4,06-013382,4554318341_R02C02,Unknown,WG9999912-DNA,B02,,,,,4554318341_R02C02,pool_id,sentrix_id
5,06-008925,4554318341_R03C01,Unknown,WG9999912-DNA,B01,,,,,4554318341_R03C01,pool_id,sentrix_id
6,06-006744,4554318341_R03C02,Unknown,WG9999912-DNA,E01,,,,,4554318341_R03C02,pool_id,sentrix_id
7,07-020190,4554318341_R04C01,Unknown,WG9999912-DNA,H01,,,,,4554318341_R04C01,pool_id,sentrix_id
8,NA10859_1,4554318341_R04C02,Unknown,WG9999912-DNA,C02,,,,,4554318341_R04C02,pool_id,sentrix_id
9,08-000591,4554318341_R05C01,Unknown,WG9999912-DNA,C01,,,,,4554318341_R05C01,pool_id,sentrix_id
10,07-015573,4554318341_R05C02,Unknown,WG9999912-DNA,F01,,,,,4554318341_R05C02,pool_id,sentrix_id
11,08-004776,4554318341_R06C01,Unknown,WG9999912-DNA,A02,,,,,4554318341_R06C01,pool_id,sentrix_id
12,NA10859_2,4554318341_R06C02,Unknown,WG9999912-DNA,D02,,,,,4554318341_R06C02,pool_id,sentrix_id
In advance, many thanks.
Christian
More information about the Bioconductor
mailing list