if (FALSE) {
ped <- read.ped("mydata")
phe <- read.phe("myphe")
cur <- pbat.work( ped )
res <- pbat.m(trait ~ NONE, ped=ped, phe=phe, fbat="gee")
...
## whatever you want to do with the results
pbat.unwork( cur ) ## sends you back to where you were
}
Run the code above in your browser using DataCamp Workspace