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


###################################################
### chunk number 2: 
###################################################
x<-read.table("Xba.rlmm")
x[1:2,1:5]


###################################################
### chunk number 3: 
###################################################
plot_theta(genotypefile="Xba.rlmm",thetafile="Xba.theta",plotfile="",snpsfile="snps.lst")