Question

我想使用Adist函数计算矩阵的所有行对之间的aitchison距离。我正在使用iris数据进行测试，我使用for脚本进行了测试，但是当我运行代码时，我收到了此消息：

adist [i]＆lt; - aDist（X [i - 1，]，X [i，]）出错：类型对象＆＃39;闭合＆＃39;不是子集表格

我感谢任何帮助。任何人都可以帮助我，提前谢谢这是代码：

X <- iris [ ,1:3]
for(i in 1: nrow (X)) {
    adist[i] <- aDist( X[i-1, ], X[i, ])
}

Answer 1

使用for循环：

library('robCompositions') 
X <- iris [ ,1:3]

a <- vector('numeric', length=nrow(X) )
for(i in 1: nrow (X)) {
  a[i] <- aDist( X[i-1, ], X[i, ])
}

更像R的方法

library('robCompositions') 
X <- iris [ ,1:3]

aDistFun <- function(x, r){
  return(aDist( x[r-1, ], x[r, ])) 
}

a <- sapply(seq(1, nrow(X)), FUN=aDistFun, x=X )

计算aitchison距离的代码出错

1 个答案: