data1<-read.table('clipboard',header=T) (其中,将EXCEL表格中的数据先复制,然后进行导入,即可),
第二步:将数据框转化为矩阵(A)
A<-as.matrix(data1)
第三步:将数据框的第一行的标题去掉
colnames(A)<-NULL
第四步:求矩阵A的转置阵,设转置阵为A1
A1<-t(A)
第五部:计算矩阵A与其转置阵A1相乘
H1<-A%*%A1
R语言如何组成矩阵#利用已有数据组建新的矩阵
v1 <- c(1:4) #创建名为v1的向量
v2 <- c(5:8) #创建名为v2的向量
m5 <- rbind(v1,v2) #将向量v1和v2按行合并为一个矩阵
m6 <- cbind(v1,v2) #将向量v1和v2按列合并为一个矩阵
#参考R语言中常见的几种创建矩阵形式总结_R语言_脚本之家 (jb51.net)
修改列名,
如果是矩阵使用colnames(chr1_edite) <- c("chrom","map","Likelihood")
如果是数据框使用names(chr1_edite) <- c("chrom","map","Likelih111")
但是数据框也可以用colnames()进行修改
x%in%y 返回逻辑向量,表示x中是否有元素在y