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 DataLab