27

ループ内のエラー (存在する場合) をスキップして、次の反復を続行したいと考えています。{0, 1, 2} からランダムにサンプリングされた要素を持つ 2 行 2 列の行列の 100 個の逆行列を計算したいと考えています。特異行列を持つことが可能です (たとえば、

1 0
2 0

これが私のコードです

set.seed(1)
count <- 1
inverses <- vector(mode = "list", 100)
repeat {
    x <- matrix(sample(0:2, 4, replace = T), 2, 2)
    inverses[[count]] <- solve(x)
    count <- count + 1
    if (count > 100) break
}

3 回目の反復で、行列は特異であり、コードは実行を停止し、エラー メッセージが表示されます。実際には、このエラーを回避して次のループに進みたいと考えています。tryまたは関数を使用する必要があることtryCatchはわかっていますが、それらの使用方法がわかりません。同様の質問がここで尋ねられましたが、それらはすべて非常に複雑であり、答えは私の理解をはるかに超えています. 誰かがこの質問専用の完全なコードを教えてくれたら、本当に感謝しています。

4

4 に答える 4

9

を使用する代わりにtryCatch、関数を使用して行列式を単純に計算できますdet。行列式がゼロの場合に限り、行列は特異です。

したがって、行列式がゼロと異なるかどうかをテストし、テストが正の場合にのみ逆行列を計算できます。

set.seed(1)
count <- 1
inverses <- vector(mode = "list", 100)
repeat {
  x <- matrix(sample(0:2, 4, replace = T), 2, 2)
  # if (det(x)) inverses[[count]] <- solve(x)
  # a more robust replacement for the above line (see comment):
  if (is.finite(determinant(x)$modulus)) inverses[[count]] <- solve(x)
  count <- count + 1
  if (count > 100) break
}

更新：

ただし、特異行列の生成を回避することは可能です。2行2列の行列式matは次のように決定されmat[1] * mat[4] - mat[3] * mat[2]ます。この知識を使用して、乱数をサンプリングできます。特異行列を生成する数をサンプリングしないでください。もちろん、これは以前にサンプリングされた数値によって異なります。

set.seed(1)
count <- 1
inverses <- vector(mode = "list", 100)

set <- 0:2 # the set of numbers to sample from

repeat {

  # sample the first value
  x <- sample(set, 1)
  # if the first value is zero, the second and third one are not allowed to be zero.
  new_set <- ifelse(x == 0, setdiff(set, 0), set)
  # sample the second and third value
  x <- c(x, sample(new_set, 2, replace = T))
  # calculate which 4th number would result in a singular matrix
  not_allowed <- abs(-x[3] * x[2] / x[1])
  # remove this number from the set
  new_set <- setdiff(0:2, not_allowed)
  # sample the fourth value and build the matrix
  x <- matrix(c(x, sample(new_set, 1)), 2, 2)

  inverses[[count]] <- solve(x)
  count <- count + 1
  if (count > 100) break
}

この手順は、生成されたすべての行列が逆行列になることを保証します。

于 2012-12-27T19:31:54.143 に答える
8

try「次の括弧内でエラーをコミットした場合はR、それをスキップして先に進みます。」

したがってx <- matrix(sample(0:2, 4, replace = T), 2, 2)、エラーが発生するのではないかと心配している場合は、次のことを行うだけです。

try(x <- matrix(sample(0:2, 4, replace = T), 2, 2))

ただし、xこれを行うと未定義になり、答えを計算できなくなることに注意してください。これは、到達したときに問題を引き起こす可能性があるため、事前にsolve(x)定義するか、すべてを「試す」ことができます。xtry

try(
      {
      x <- matrix(sample(0:2, 4, replace = T), 2, 2)
      inverses[[count]] <- solve(x)
      }
    )
于 2012-12-27T19:11:34.370 に答える
5

try のドキュメントは、問題をかなりよく説明しています。完全に通過することをお勧めします。

Edit:ドキュメントの例は非常に簡単で、op の質問に非常によく似ていました。しかし、提案をありがとう。ドキュメントページの例に続く答えは次のとおりです。

# `idx` is used as a dummy variable here just to illustrate that
# all 100 entries are indeed calculated. You can remove it.
set.seed(1)
mat_inv <- function(idx) {
    print(idx)
    x <- matrix(sample(0:2, 4, replace = T), nrow = 2)
    solve(x)
}
inverses <- lapply(1:100, function(idx) try(mat_inv(idx), TRUE))
于 2012-12-27T18:56:42.627 に答える