R语言实现ochiia系数

#读取数据
mydata <- read.csv("C:\\Users\\NAU\\Desktop\\app18.csv")


#查看数据类型
mode(mydata)


#数据格式转换 列表（list）转换矩阵（matrix）
mydata=as.matrix(mydata)


#计算公式
> m <- matrix(mydata, nrow = 72, ncol = 72, byrow = FALSE)
> n <- matrix(1:5184, nrow = 72, ncol = 72, byrow = FALSE)
> for(i in 1:72){
+  for(j in 1:72){
+   n[i,j]=m[i,j]/(m[i,i]*m[j,j])^(1/2)
+   }
+ }


#最大输出行
options(max.print=<数值，比如100000000>)?


#输出
> n


#手动输入数据的要点是（不要打后面的括号）
a<-c(341,45,36,66,108,116,45,242,92,7,0,11,36,92,240,15,33,5,66,7,15,219,22,10,108,0,33,22,199,0,116,11,5,10,0,163)
m<-matrix(a,nrow=6,ncol=6,byrow=FALSE)
p<-matrix(1:36,nrow=6,ncol=6,byrow=FALSE)
for(i in 1:6) {
  for(j in 1:6){
    p[i,j]=m[i,j]/(m[i,i]*m[j,j])^(1/2)
    }
  }