setwd("G:/¿ÆÑÐ/methy-evo/20140911/PCA/")
require(graphics)
data<-read.table("gene_mean_var.txt",row.names=1,header=T)
data_new<-matrix(nrow=21231,ncol=2)
colnames(data_new)<-c("Mean","Var")
pc.res<-princomp(data,cor=T,scores=T)
pc.score<-pc.res$scores
write.table(pc.score,"PCA_result.txt",row.names=T,col.names=F,quote=F,sep="\t")