library(data.table)
  snpdat =(fread("data2/1:145538184-145786438.gen.SNPDAT", header=T))
  ld = as.matrix(setDF(fread("data2/1:145538184-145786438.gen.LD", header=F)))

  g = prep_dataset_coco(snpdat,ld,var_y = 1.6421,exact = F,hwe_variance = T)
  step_coc = stepwise_coco(g,joint = F)
  joint_from_ids(g,"rs1471633",return_only_these = T)


theboocock/coco documentation built on May 31, 2019, 9:10 a.m.