###################################################
### chunk number 1: 
###################################################
library(multiscan)


###################################################
### chunk number 2:  eval=FALSE
###################################################
## vignette("multiscan")


###################################################
### chunk number 3: 
###################################################
data(murine)


###################################################
### chunk number 4: 
###################################################
murine[1:10,]


###################################################
### chunk number 5:  eval=FALSE
###################################################
## help(murine)


###################################################
### chunk number 6: 
###################################################
data(murine)
fit<-multiscan(murine)
fit


###################################################
### chunk number 7: 
###################################################
str(fit)


###################################################
### chunk number 8: 
###################################################
gene.exprs<-fit$mu


###################################################
### chunk number 9:  eval=FALSE
###################################################
## plot(fit)


###################################################
### chunk number 10: fitted
###################################################
plot(fit)


###################################################
### chunk number 11:  eval=FALSE
###################################################
## help(multiscan)


###################################################
### chunk number 12:  eval=FALSE
###################################################
## op<-par(mfrow=c(2,2))
## plot(fit,residual=TRUE)
## par(op)


###################################################
### chunk number 13: sdres
###################################################
op<-par(mfrow=c(2,2))
plot(fit,residual=TRUE)
par(op)