if (length(dir(pattern=".cel$")) != 0) {
## Read in the *.CEL files
abatch <- ReadAffy()
## default calling method
SScores <- SScore(abatch)
## specifying SF and SDT (gives same results as above)
SfSdt <- computeSFandSDT(abatch)
SScores <- SScore(abatch,SF=SfSdt$SF,SDT=SfSdt$SDT)
## specifying outlier and masked values should be included in calculations
SScores <- SScore(abatch,rm.outliers=FALSE,rm.mask=FALSE)
## round results to 3 significant digits
SScores <- SScore(abatch,digits=3)
## show verbose output
SScores <- SScore(abatch,verbose=TRUE)
}
Run the code above in your browser using DataLab