#计算相关系数
setwd('C:/Users/Tao/Desktop')
expr=read.delim('TPM1.txt',header=T,row.names = 1,sep = "\t",na.strings = '-') #读取表达矩阵
r=cor(expr,method='pearson',use='pairwise')
# method可选 pearson 或者spearman
# use="pairwise" 表示剔除缺失
write.table(r,file='pearson.r.txt',col.names = NA,quote=F,sep="\t")
网友评论