multiscan.R
In multiscan: R package for combining multiple scans

### 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)