################################################### ### chunk number 1: readfile1 eval=FALSE ################################################### ## ## cset <- readCodelinkSet("targets.txt") ## ################################################### ### chunk number 2: annotations ################################################### library(codelink) data(codelink.exprset) annotation(codelink.exprset) annotation(codelink.exprset) <- "rwgcod" annotation(codelink.exprset) ################################################### ### chunk number 3: preprocess ################################################### cset <- codPreprocess(codelink.exprset) ################################################### ### chunk number 4: plot ################################################### codPlot(cset, what = "density") codPlot(cset, what = "ma") codPlot(cset, what = "image") ################################################### ### chunk number 5: fit eval=FALSE ################################################### ## ## design <- model.matrix(~ -1 + factor(cset$Treatment)) ## fit <- lmFit(cset, design) ## ################################################### ### chunk number 6: fit eval=FALSE ################################################### ## ## writeCodelink(codelink.exprset, file = "intensities.txt") ## ################################################### ### chunk number 7: fit eval=FALSE ################################################### ## ## foo.avg <- averageProbes(codelink.exprset) ## ################################################### ### chunk number 8: sessioninfo ################################################### sessionInfo()