require(beadarrayExampleData)
data(exampleSummaryData)
exampleSummaryDatalog2 <- channel(exampleSummaryData, "G")
exampleSummaryDataNorm <- normaliseIllumina(exampleSummaryDatalog2,
method = "quantile", transform = "none")
aaa<-na.omit(data.frame(I(rownames(exprs(exampleSummaryDataNorm))),exprs(exampleSummaryDataNorm)))
ProbeID <- aaa[,1][1:100]
eSet <- na.omit(exprs(exampleSummaryDataNorm))[1:100,]
stddev <- na.omit(se.exprs(exampleSummaryDataNorm))[1:100,]
nSet <- na.omit(attributes(exampleSummaryDataNorm)$assayData$nObservations)[1:100,]
seSet<-stddev/sqrt(nSet)
iccResults<-iccFun(eSet,seSet,nSet,ProbeID =ProbeID ,iccQuant=c(0,0.25,0.5,1),
diffIcc=TRUE,keepData=TRUE)
emCDF(iccResults, iccQuant= c(0,0.25,0.5,1))
Run the code above in your browser using DataLab