test/test_DS.R

In tpbilton/GUSrelate: Genotyping Uncertainty with Sequencing data and RELATEdness

library(GUSrelate)

ra = readRA("~/Academia/Otago_University/GUSrelate/test/simDS.vcf.ra.tab")
grm = makeGRM(ra, samfile = "~/Academia/Otago_University/GUSrelate/test/test_sim_ploid.csv")

grm$addSampleInfo("~/Academia/Otago_University/GUSrelate/test/test_sim_addSamInfo.csv")

grm$HWEtest(nThreads = 3)
grm$computeGRM(name = "VR", ep=0)

grm$computeGRM(name = "WG", method="WG", ep=0)