search for: affectionstatus

Displaying 1 result from an estimated 1 matches for "affectionstatus".

2011 Jul 14
2
R package: pbatR
...ibrary(tcltk) pbat.set("C:/pbat") x <- data.frame(pid = c(1,1,1,2,2,2,2,3,3,3), # three families id = c(1,2,3,1,2,3,4,1,2,3), idfath = c(0,0,1,0,0,1,1,0,0,1), idmoth = c(0,0,2,0,0,2,2,0,0,2), sex = c(1,2,1,1,2,2,2,1,2,1), AffectionStatus = c(0,0,1,0,0,1,0,0,0,1), # 1 for case, 0 for control m1.1 = c(1,1,2,2,1,1,2,2,2,1), # two SNPs with two columns for each SNP m1.2 = c(1,2,1,2,1,2,1,1,2,2), m2.1 = c(1,1,2,2,2,1,1,1,2,1), m2.2 = c(2,1,2,1,2,2,2,1,1,1)) x1 <- as.ped(x)...