library(CustomCDF)
require(RPA)
#if (!requireNamespace("BiocManager", quietly=TRUE))
#install.packages("BiocManager")
#BiocManager::install("hgu133ahsensgcdf")
#BiocManager::install("hgu133ahsenstcdf")
set.seed(11122)
#cels <- list.celfiles("/share/mi/data/GSE3526/CEL/", full.names = T)[1:3]
cels <- list.celfiles("/share/mi/data/GSE3307/", full.names = T)[1:3]
abatch <- ReadAffy(filenames = cels)
#abatch@cdfName <- "HGU133A_Hs_ENSG"
#abatch@cdfName <- "HGU133A_Hs_ENST"
#HGU133A_Hs_ENSG
np1 <- length(unlist(pmindex(abatch)))
print(paste("Total probes:", np1))
xy <- getsnpprobe(chip = "HG-U133A")
report <- removeprobe(abatch, xyMatrix = xy)
np2 <- length(unlist(pmindex(abatch)))
print(paste("Total probes:", np2))
print(paste("Difference: ", np1 - np2))
print( c(np1, np2, np1-np2))
# With original affy sets (cdfName not set)
#> 247965 - 245657
#[1] 2308
# length(featureNames(abatch))
# ~ 22 283
# ENSG
#> c(np1, np2, np1-np2)
#[1] 171666 171651 15
# length(featureNames(abatch))
#[1] 11869
# ENST
#> print( c(np1, np2, np1-np2))
#[1] 701051 696521 4530
#length(featureNames(abatch))
# > 51956
#> dim(xy)
#[1] 2681 2
In total 2681 SNPs
Probesets SNP probes
AFFY 22,283 2308
ENSG 11,869 15
ENST 51,956 4530
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.