比较 R 中的两个值

Question

我正在检查矩阵的输入在 R 中是倒数的，即一侧的值是 = 到 1/值..

到目前为止我有..

    AHP <- function(pairwisematrix){

    ## check the input data

    if (!((is.matrix(pairwisematrix) || (is.data.frame(pairwisematrix))))) 
            stop("pairwise matrix must be a matrix or a data frame")

    if (!(nrow(pairwisematrix) == ncol(pairwisematrix))) 
      stop("pairwise matrix must be a square matrix or a data frame")

    for (i in 1:nrow(pairwisematrix))   {
    for (j in 1:ncol(pairwisematrix))       {
        if (i == j) {   if (pairwisematrix[i,j] != 1)                                                                       {   stop("the central values in the reciprocal matrix must be 1.")                              }}
        else { if ((pairwisematrix[i,j] == 1 / pairwisematrix[j,i]) || (pairwisematrix[j,i] == 1 / pairwisematrix[i,j]))))  {   stop("the pairwise matrix must be reciprocal (i.e. value on one side must = 1/value)")     }}
                                            }
                                        }


  out <- "all worked"

    return(out)
}

但是当我测试时：

check1 <- matrix(c(1,1/3,5,3,1,0.2,0.2,5,1),3,3,byrow=T)
test <- AHP(check1)

我得到错误：

成对矩阵必须是倒数（即一侧的值必须 = 1/值)0.333333 & 0.333333

我尝试将值转换为字符串、部分字符串并尝试相同（a，b，）但没有成功。

有人有什么想法吗？

Answer 1

在双 for loops 中有许多 if s，如果它按预期工作，我会感到惊讶。 R 是为处理矩阵而设计的，所以你可以写类似

AHP <- function(pairwisematrix){
  if(!all(pairwisematrix == t(1/pairwisematrix)))
    stop("the pairwise matrix must be reciprocal (i.e. value on one side must = 1/value)")
  else
    return("all worked")
}

AHP(check1)
#[1] "all worked"