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)...