r在R中合并PLINK表型和协变量

篇首语：本文由编程笔记#小编为大家整理，主要介绍了r 在R中合并PLINK表型和协变量相关的知识，希望对你有一定的参考价值。

famfile <- read.table(&＃39;/home/cnieverg/MEG2/qc/pts_meg2_mix_am.fam&＃39;,header=F,stringsAsFactors=F)
names(famfile) <- c("FID","IID","F","M","Gender","Pheno")
mds <- read.table(&＃39;/home/cnieverg/MEG2/qc/pca/pts_meg2_mix_am-qc-eu1_pca.menv.mds_cov&＃39;,header=T,stringsAsFactors=F)
dat <- merge(famfile,mds,by=c("FID","IID"))
#Number of analyzed subjects of European ancestry
dim(dat)[1]
#Number of analyzed cases/controls of European ancestry, overall
table(dat$Pheno)
#Cross tabs for pool of analyzed subjects, by gender
table(dat$Gender,dat$Pheno)