Nothing
##########################################################################
# Allele-specific CRMAv2
##########################################################################
library("aroma.affymetrix");
verbose <- Arguments$getVerbose(-8, timestamp=TRUE);
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# Setup
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
dataSet <- "GSE13372";
chipType <- "GenomeWideSNP_6,Full";
cdf <- AffymetrixCdfFile$byChipType("GenomeWideSNP_6", tags="Full")
print(cdf)
csR <- AffymetrixCelSet$byName(dataSet, chipType=chipType, cdf=cdf);
print(csR);
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# AS-CRMAv2
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
## Extract only one tumor and two matching normals
nms <- c(T="GSM337641", N="GSM337662", R="GSM337663")
idxs <- indexOf(csR, nms)
csE <- extract(csR, idxs)
res <- doASCRMAv2(csE, verbose=verbose)
# nn <- names(csR)
# resR <- doASCRMAv2(csR[which(nn=="GSM337641")], verbose=verbose);
# resN <- doASCRMAv2(csR[which(nn=="GSM337663")], verbose=verbose);
# resRef <- doASCRMAv2(csR[which(nn=="GSM337662")], verbose=verbose);
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.