fmg_hmp<-read.table("G:/科研/methy-evo/20140107/fmg_hmp.txt")
colnames(fmg_hmp)<-c("GeneID","ASW","CEU","CHB","CHD","GIH","JPT","LWK","MEX","MKK","TSI","YRI")
smg_hmp<-read.table("G:/科研/methy-evo/20140107/smg_hmp.txt")
colnames(smg_hmp)<-c("GeneID","ASW","CEU","CHB","CHD","GIH","JPT","LWK","MEX","MKK","TSI","YRI")
quan_smg_hmp<-matrix(nrow=11,ncol=3)
quan_fmg_hmp<-matrix(nrow=11,ncol=3)
for (i in 2:12){
  quan_smg_hmp[i-1,1]<-quantile(smg_hmp[,i],na.rm=TRUE)[3]
  quan_smg_hmp[i-1,2]<-quantile(smg_hmp[,i],na.rm=TRUE)[4]
  quan_smg_hmp[i-1,3]<-quantile(smg_hmp[,i],na.rm=TRUE)[2]
  quan_fmg_hmp[i-1,1]<-quantile(fmg_hmp[,i],na.rm=TRUE)[3]
  quan_fmg_hmp[i-1,2]<-quantile(fmg_hmp[,i],na.rm=TRUE)[4]
  quan_fmg_hmp[i-1,3]<-quantile(fmg_hmp[,i],na.rm=TRUE)[2]
}

fmg_genome<-read.table("G:/科研/methy-evo/20140107/fmg_r2_genome.txt")
colnames(fmg_genome)<-c("GeneID","ASW","CEU","CHB","CHS","CLM","FIN","GBR","JPT","LWK","MXL","PUR","TSI","YRI")
smg_genome<-read.table("G:/科研/methy-evo/20140107/smg_r2_genome.txt")
colnames(smg_genome)<-c("GeneID","ASW","CEU","CHB","CHS","CLM","FIN","GBR","JPT","LWK","MXL","PUR","TSI","YRI")
quan_smg_genome<-matrix(nrow=13,ncol=3)
quan_fmg_genome<-matrix(nrow=13,ncol=3)
for (i in 2:14){
  quan_smg_genome[i-1,1]<-quantile(smg_genome[,i],na.rm=TRUE)[3]
  quan_smg_genome[i-1,2]<-quantile(smg_genome[,i],na.rm=TRUE)[4]
  quan_smg_genome[i-1,3]<-quantile(smg_genome[,i],na.rm=TRUE)[2]
  quan_fmg_genome[i-1,1]<-quantile(fmg_genome[,i],na.rm=TRUE)[3]
  quan_fmg_genome[i-1,2]<-quantile(fmg_genome[,i],na.rm=TRUE)[4]
  quan_fmg_genome[i-1,3]<-quantile(fmg_genome[,i],na.rm=TRUE)[2]
}