################################################### ### chunk number 1: setup1 ################################################### #line 79 "vignettes/cellHTS/inst/doc/cellhts.Rnw" library("cellHTS") ################################################### ### chunk number 2: setup2 ################################################### #line 83 "vignettes/cellHTS/inst/doc/cellhts.Rnw" ## for debugging: options(error=recover) ## for software development, when we do not want to install ## the package after each minor change: ## for(f in dir("~/huber/projects/Rpacks/cellHTS/R", full.names=TRUE, pattern=".R$"))source(f) ################################################### ### chunk number 3: dataPath ################################################### #line 126 "vignettes/cellHTS/inst/doc/cellhts.Rnw" experimentName = "KcViab" dataPath=system.file(experimentName, package="cellHTS") ################################################### ### chunk number 4: dirDataPath ################################################### #line 137 "vignettes/cellHTS/inst/doc/cellhts.Rnw" dataPath rev(dir(dataPath))[1:12] ################################################### ### chunk number 5: readPlateData ################################################### #line 144 "vignettes/cellHTS/inst/doc/cellhts.Rnw" x = readPlateData("Platelist.txt", name=experimentName, path=dataPath) ################################################### ### chunk number 6: showX ################################################### #line 148 "vignettes/cellHTS/inst/doc/cellhts.Rnw" x ################################################### ### chunk number 7: plateFileTable ################################################### #line 159 "vignettes/cellHTS/inst/doc/cellhts.Rnw" cellHTS:::tableOutput(file.path(dataPath, "Platelist.txt"), "plate list") cellHTS:::tableOutput(file.path(dataPath, names(x$intensityFiles)[1]), "signal intensity", header=FALSE, dropColumns=1) ################################################### ### chunk number 8: writeReport1Show eval=FALSE ################################################### ## #line 191 "vignettes/cellHTS/inst/doc/cellhts.Rnw" ## out = writeReport(x) ################################################### ### chunk number 9: writeReport1Do eval=FALSE ################################################### ## #line 194 "vignettes/cellHTS/inst/doc/cellhts.Rnw" ## out = writeReport(x, force=TRUE) ################################################### ### chunk number 10: printout eval=FALSE ################################################### ## #line 202 "vignettes/cellHTS/inst/doc/cellhts.Rnw" ## out ################################################### ### chunk number 11: browseReport1 eval=FALSE ################################################### ## #line 208 "vignettes/cellHTS/inst/doc/cellhts.Rnw" ## browseURL(out) ################################################### ### chunk number 12: annotatePlateRes ################################################### #line 236 "vignettes/cellHTS/inst/doc/cellhts.Rnw" x = configure(x, "Plateconf.txt", "Screenlog.txt", "Description.txt", path=dataPath) ################################################### ### chunk number 13: plateConfscreenLogTable ################################################### #line 250 "vignettes/cellHTS/inst/doc/cellhts.Rnw" cellHTS:::tableOutput(file.path(dataPath, "Plateconf.txt"), "plate configuration", selRows=25:28) cellHTS:::tableOutput(file.path(dataPath, "Screenlog.txt"), "screen log", selRows=1:3) ################################################### ### chunk number 14: ################################################### #line 292 "vignettes/cellHTS/inst/doc/cellhts.Rnw" table(x$plateConf$Content) ################################################### ### chunk number 15: normalizePlateMedian ################################################### #line 346 "vignettes/cellHTS/inst/doc/cellhts.Rnw" x = normalizePlates(x, normalizationMethod="median") ################################################### ### chunk number 16: summarizeReplicates ################################################### #line 391 "vignettes/cellHTS/inst/doc/cellhts.Rnw" x = summarizeReplicates(x, zscore="-", summary="mean") ################################################### ### chunk number 17: boxplotzscore ################################################### #line 400 "vignettes/cellHTS/inst/doc/cellhts.Rnw" ylim = quantile(x$score, c(0.001, 0.999), na.rm=TRUE) boxplot(x$score ~ x$wellAnno, col="lightblue", outline=FALSE, ylim=ylim) ################################################### ### chunk number 18: normalizePlateMedianWithZscore ################################################### #line 423 "vignettes/cellHTS/inst/doc/cellhts.Rnw" xalt = normalizePlates(x, normalizationMethod="median", zscore="-") xalt = summarizeReplicates(xalt, summary="mean") ################################################### ### chunk number 19: geneIDs ################################################### #line 451 "vignettes/cellHTS/inst/doc/cellhts.Rnw" x = annotate(x, geneIDFile="GeneIDs_Dm_HFA_1.1.txt", path=dataPath) ################################################### ### chunk number 20: geneIDsTable ################################################### #line 455 "vignettes/cellHTS/inst/doc/cellhts.Rnw" cellHTS:::tableOutput(file.path(dataPath, "GeneIDs_Dm_HFA_1.1.txt"), "gene ID", selRows = 3:6) ################################################### ### chunk number 21: printxagain ################################################### #line 477 "vignettes/cellHTS/inst/doc/cellhts.Rnw" x ################################################### ### chunk number 22: savex ################################################### #line 483 "vignettes/cellHTS/inst/doc/cellhts.Rnw" save(x, file=paste(experimentName, ".rda", sep=""), compress=TRUE) ################################################### ### chunk number 23: writeReport2 eval=FALSE ################################################### ## #line 491 "vignettes/cellHTS/inst/doc/cellhts.Rnw" ## out = writeReport(x, force=TRUE, ## plotPlateArgs = list(xrange=c(0.5, 1.5)), ## imageScreenArgs = list(zrange=c( -2, 6.5), ar=1)) ################################################### ### chunk number 24: browseReport2 eval=FALSE ################################################### ## #line 499 "vignettes/cellHTS/inst/doc/cellhts.Rnw" ## browseURL(out) ################################################### ### chunk number 25: imageScreen eval=FALSE ################################################### ## #line 551 "vignettes/cellHTS/inst/doc/cellhts.Rnw" ## imageScreen(x, ar=1, zrange=c(-3,4)) ################################################### ### chunk number 26: exportData eval=FALSE ################################################### ## #line 616 "vignettes/cellHTS/inst/doc/cellhts.Rnw" ## writeTab(x, file="Data.txt") ################################################### ### chunk number 27: exportOtherData eval=FALSE ################################################### ## #line 628 "vignettes/cellHTS/inst/doc/cellhts.Rnw" ## # determine the ratio between each well and the plate median ## y = array(as.numeric(NA), dim=dim(x$xraw)) ## nrWell = dim(x$xraw)[1] ## for(p in 1:(dim(x$xraw)[2])) { ## samples = (x$wellAnno[(1:nrWell)+nrWell*(p-1)]=="sample") ## y[, p, , ] = apply(x$xraw[, p, , , drop=FALSE], 3:4, ## function(w) w/median(w[samples], na.rm=TRUE)) } ## y=signif(y, 4) ## out = matrix(y, nrow=prod(dim(y)[1:2]), ncol=dim(y)[3:4]) ## out = cbind(x$geneAnno, x$wellAnno, out) ## colnames(out) = c(names(x$geneAnno), "wellAnno", ## sprintf("Well/Median_r%d_ch%d", rep(1:dim(y)[3], dim(y)[4]), ## rep(1:dim(y)[4], each=dim(y)[3]))) ## write.tabdel(out, file="WellMedianRatio.txt") ################################################### ### chunk number 28: transfplots ################################################### #line 686 "vignettes/cellHTS/inst/doc/cellhts.Rnw" library("vsn") par(mfcol=c(3,2)) myPlots=function(z,...) { hist(z[,1], 100, col="lightblue", xlab="",...) meanSdPlot(z, ylim=c(0, quantile(abs(z[,2]-z[,1]), 0.95, na.rm=TRUE)), ...) qqnorm(z[,1], pch='.', ...) qqline(z[,1], col='blue') } dv = matrix(x$xnorm, nrow=prod(dim(x$xnorm)[1:2]), ncol=dim(x$xnorm)[3]) myPlots(dv, main="untransformed") xlog = normalizePlates(x, normalizationMethod="median", transform=log2) dvlog = matrix(xlog$xnorm, nrow=prod(dim(xlog$xnorm)[1:2]), ncol=dim(xlog$xnorm)[3]) myPlots(dvlog, main="log2") ################################################### ### chunk number 29: sessionInfo ################################################### #line 710 "vignettes/cellHTS/inst/doc/cellhts.Rnw" toLatex(sessionInfo())