# Download data files
load_datafiles()
# mQTL mapping results
load(results)
## Plot 3D profile
dev.new(width=5,height=5,pointsize=5)
probs=c(0.95,0.99,0.999,0.9999) ## probabilities
pplot(results$res,"Full 2D Profile", results$ppm, results$best, quantile(results$res,probs=probs))
Run the code above in your browser using DataLab