第一步：導(dǎo)入數(shù)據(jù)

?data1<-read.table('clipboard',header=T)? （其中，將EXCEL表格中的數(shù)據(jù)先復(fù)制，然后進(jìn)行導(dǎo)入，即可），

第二步：將數(shù)據(jù)框轉(zhuǎn)化為矩陣（A）

A<-as.matrix(data1)??

第三步：將數(shù)據(jù)框的第一行的標(biāo)題去掉

?colnames(A)<-NULL

第四步：求矩陣A的轉(zhuǎn)置陣，設(shè)轉(zhuǎn)置陣為A1

A1<-t(A)

第五部：計(jì)算矩陣A與其轉(zhuǎn)置陣A1相乘

?H1<-A%*%A1