Question: How to get SNP 6.0 array data in R?
gravatar for agicict
4.1 years ago by
Korea, Republic Of
agicict150 wrote:



I currently downloaded a number of SNP datasets from GEO (CCLE)


Also I downloaded Genomewide affymetrix 6.0 SNP array annotation data from affymetrix site.


I imported the SNP 6.0 CEL files using affy package's ReadAffy.

Strangely, every CEL files has approximately 7,000,000 rows but affymetrix annotation data has around 900,000 rows.

What is right way to match probes?

annotation snp R • 3.6k views
ADD COMMENTlink modified 4.1 years ago by Henrik Bengtsson60 • written 4.1 years ago by agicict150
gravatar for Henrik Bengtsson
4.1 years ago by
United States
Henrik Bengtsson60 wrote:

The GenomeWideSNP_6 chip type has 6,892,960 probes (never changes). These probes are arranged in probe sets ("units") corresponding to 934,946 bi-allelic SNPs and 946,371 non-polymorphic single-probe CN loci. The exact number differ slightly between genome builds.

Check out the aroma.affymetrix R package, cf. It has several ready pipelines for GenomeWideSNP_6, especially for copy-number analysis.

(author of aroma.affymetrix)

ADD COMMENTlink modified 7 weeks ago by RamRS25k • written 4.1 years ago by Henrik Bengtsson60

Thank you for your kind explanation.

ADD REPLYlink written 4.1 years ago by agicict150
gravatar for Jan Oosting
4.1 years ago by
Jan Oosting910
Leiden, NL
Jan Oosting910 wrote:

You could use the crlmm Bioconductor package

ADD COMMENTlink written 4.1 years ago by Jan Oosting910
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1869 users visited in the last hour