### R code from vignette source 'multiscan.Rnw' ################################################### ### code chunk number 1: multiscan.Rnw:70-71 ################################################### library(multiscan) ################################################### ### code chunk number 2: multiscan.Rnw:85-86 (eval = FALSE) ################################################### ## vignette("multiscan") ################################################### ### code chunk number 3: multiscan.Rnw:107-108 ################################################### data(murine) ################################################### ### code chunk number 4: multiscan.Rnw:113-114 ################################################### murine[1:10,] ################################################### ### code chunk number 5: multiscan.Rnw:119-120 (eval = FALSE) ################################################### ## help(murine) ################################################### ### code chunk number 6: multiscan.Rnw:128-131 ################################################### data(murine) fit<-multiscan(murine) fit ################################################### ### code chunk number 7: multiscan.Rnw:135-136 ################################################### str(fit) ################################################### ### code chunk number 8: multiscan.Rnw:142-143 ################################################### gene.exprs<-fit$mu ################################################### ### code chunk number 9: multiscan.Rnw:149-150 (eval = FALSE) ################################################### ## plot(fit) ################################################### ### code chunk number 10: fitted ################################################### plot(fit) ################################################### ### code chunk number 11: multiscan.Rnw:168-169 (eval = FALSE) ################################################### ## help(multiscan) ################################################### ### code chunk number 12: multiscan.Rnw:178-181 (eval = FALSE) ################################################### ## op<-par(mfrow=c(2,2)) ## plot(fit,residual=TRUE) ## par(op) ################################################### ### code chunk number 13: sdres ################################################### op<-par(mfrow=c(2,2)) plot(fit,residual=TRUE) par(op)