[BioC] how to get oligoSnpSet object?

Yupu Liang liang at cbio.mskcc.org
Fri Jul 18 23:18:41 CEST 2008


Hello,

  I am new to snp array and I am trying to work on affymetrix's  
SNP6.0 data through oligo package. But I couldn't find a way to  
either get the copy number or get the data into oligoSnpSet object:

library(oligo)
fullFilenames <- list.celfiles(path=getwd(),full.name=TRUE)

#normalization and summarization
data <- just SNPRMA(fullFilenames)

#genotype calling
outputDir <- file.path(getwd(),'crlmmTest')
crlmm(fullFilenames,outputDir)
crlmmOut <-getCrlmmSummaries(outputDir)



 >data
SnpCnvQSet (storageMode: lockedEnvironment)
assayData: 906600 features, 7 samples
   element names: thetaA, thetaB
phenoData
   sampleNames: CA_SNP_Array_6_Gist194.CEL,  
CA_SNP_Array_6_Gist237_11.CEL, ..., C
   A_SNP_Array_6_Gist363_2.CEL  (7 total)
   varLabels and varMetadata description: none
featureData
   featureNames: SNP_A-1780270, SNP_A-1780271, ..., SNP_A-8717615   
(906600 total)
   fvarLabels and fvarMetadata description: none
experimentData: use 'experimentData(object)'
Annotation: pd.genomewidesnp.6

 > crlmmOut
SnpCnvCallSetPlus (storageMode: lockedEnvironment)
assayData: 906600 features, 7 samples
   element names: calls, callsConfidence, thetaA, thetaB
phenoData
   sampleNames: CA_SNP_Array_6_Gist194.CEL,  
CA_SNP_Array_6_Gist237_11.CEL, ..., C
   A_SNP_Array_6_Gist363_2.CEL  (7 total)
   varLabels and varMetadata description: none
featureData
   featureNames: SNP_A-4270094, SNP_A-8282305, ..., SNP_A-8433021   
(906600 total)
   fvarLabels and fvarMetadata description: none
experimentData: use 'experimentData(object)'
Annotation: pd.genomewidesnp.6


Am I doing something Wrong? Have I missed anything?

Thanks,
Yupu



More information about the Bioconductor mailing list