###load data
data(matchBoxExpression)
###the column name for the identifiers and the ranking statistics
idCol <- "SYMBOL"
byCol <- "t"
####filter the redundant features using SYMBOL and t-statistics
matchBoxExpression <- lapply(matchBoxExpression, filterRedundant, idCol=idCol, byCol=byCol)
###select and merge into a matrix
mat <- mergeData(matchBoxExpression, idCol=idCol, byCol=byCol)
###COMPUTE CAT
cpH2L <- computeCat(mat, idCol=1, size=round(nrow(mat)/1),
decreasing=TRUE, method="equalRank")
###CATplot without probability intervals
par(mar=c(3,3,2,1))
plotCat(cpH2L, main="CAT-plot, decreasing t-statistics",
cex=1, lwd=2, cexPts=1.5, spacePts=15,
legend=TRUE, where="center",
legCex=1, ncol=1)
###compute probability intervals
confInt <- calcHypPI(data=mat)
###CATplot with probability intervals
par(mar=c(3,3,2,1))
plotCat(cpH2L, main="CAT-plot, decreasing t-statistics, probability intevals",
cex=1, lwd=2, cexPts=1.5, spacePts=15,
legend=TRUE, where="center",
legCex=1, ncol=1)
Run the code above in your browser using DataLab