# NOT RUN {
data(res)
data(design)
data(covset)
d <- data.frame(cbind(y = res$ori.res$v$E[1,] ,Subject = covset$ear,
Time = covset$time, w = 1/res$ori.res$v$weights[1,], design))
glsout <- rmRNAseq:::glsCAR1(d)
glsout
# }
Run the code above in your browser using DataLab