# # # R代码从装饰图案的data-input-vignette来源。Rnw“# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块1号:自由# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #要求(snpStats) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块2号:sysfile # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # pedfile < -执行(“extdata / sample.ped。广州”,包= " snpStats”) pedfile infofile < - - - - - -执行(“extdata / sample.info”、包=“snpStats”) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块3号:redpedfile # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #示例< -阅读。pedfile (pedfile snp = infofile) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块数量4:sample1 # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #样本基因型col.summary美元美元(样本基因型美元)加头(col.summary(样本基因型)美元)# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块5号:sample2 # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #头(样本fam)美元# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块6号:sample3 # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #头(示例图美元)# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块7号:叮铃声# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # fam < -执行(“extdata /样品。fam公司”,包= " snpStats”) bim < -执行(“extdata /样品。荡妇”,包= " snpStats床”)< -执行(“extdata /样品。床”,包= " snpStats”)示例< -阅读。叮铃声(床,荡妇,fam) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块8号:plinkout # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #美元样本基因型col.summary(样本基因型美元)加器头部(样本fam)美元美元(美元样本地图)# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块9号:plinkselect # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # < -读子集。叮铃声(床,荡妇,家人,select.snps = 6:10)基因型col.summary美元子集(子集基因型美元)加器子集映射美元# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块10号:longfile # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # longfile < -执行(“extdata / sample-long。广州”,包= " snpStats”) longfile # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块11号:入围# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #猫(readline (longfile 5), 9 =“\ n”) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块12号:readlong # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # gdata < -阅读。长(longfile字段= c (snp = 1,样品= 2,基因型= 3,信心= 4),gcodes = c(“1”、“2”、“3”),阈值= 0.95)gdata总结(gdata) # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #代码块13号:readlongallele # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # allelesfile < -执行(“extdata / sample-long-alleles。广州”,包= " snpStats”)猫(readline (allelesfile 5), 9 =“\ n”) gdata < -阅读。长(allelesfile字段= c (snp = 1,样品= 2,等位基因。= 3,等位基因。B = 4,信心= 5)阈值= 0.95)gdata gdata基因型gdata等位基因美元