R 中的 for 循环之内或之外 - 计算矩阵的对角积

发布于 2024-09-26 17:00:42 字数 402 浏览 2 评论 0原文

我试图找到 20x20 矩阵中 2 位数字的最大对角积。

这给出了一条错误消息：

i <- 17:1
z <- for (j in 1:(18-i))
        {b <- max ((x[i,j]*x[i+1,j+1]*x[i+2,j+2]*x[i+3,j+3]))}}

但这没有：

z <- for (i <- 17:1)
{for (j in 1:(18-i))
        {b <- max ((x[i,j]*x[i+1,j+1]*x[i+2,j+2]*x[i+3,j+3]))}}

但第二个版本给我的数字太小了。为什么第一个不起作用，我认为它会产生正确的答案，但我不明白错误消息。

原文

I'm trying to find the maximum diagonal product of 2 digit numbers in a 20x20 matrix.

This gives an error message :

i <- 17:1
z <- for (j in 1:(18-i))
        {b <- max ((x[i,j]*x[i+1,j+1]*x[i+2,j+2]*x[i+3,j+3]))}}

But this doesn't:

z <- for (i <- 17:1)
{for (j in 1:(18-i))
        {b <- max ((x[i,j]*x[i+1,j+1]*x[i+2,j+2]*x[i+3,j+3]))}}

but the second version gives me a number too small . Why does the first one not work, i think it would yield the correct answer, but i don't understand the error message.

分享到QQ

分享到微博

如果你对这篇内容有疑问，欢迎到本站社区发帖提问参与讨论，获取更多帮助，或者扫码二维码加入 Web 技术交流群。

发布评论

需要登录才能够评论，你可以免费注册一个本站的账号。

颜漓半夏 2024-10-03 17:00:42

这看起来是错误的。

您不能将 for 循环的结果分配给变量。 max() 是在一个标量变量上，这是无意义的。最后，未指定矩阵 x 。我会用更小的东西重试，甚至可能将一些临时结果打印到屏幕上。

先走再跑仍然是个好建议。稍后您仍然可以矢量化以获得冲刺解决方案。

回复收藏 0 原文

单挑你×的.吻 2024-10-03 17:00:42

实际上，与 Dirk 相反，我认为您应该尽快了解 R 中的矢量化。您尝试实现的循环结构远非最佳，而且实际上是多余的。仅应在非常特殊的情况下使用 for 循环。查看此问题中的讨论。查看一下便利函数的帮助文件，例如 diag()、combn()、prod() 和 apply()。

将它们组合起来做你想做的事情很容易：

x <-matrix(1:400,ncol=20)

Diag <- diag(x)

Id <- combn(1:length(Diag),2)

Diag.prod <- apply(matrix(Diag[Id],ncol=2),1,prod)

Max.Diag.prod <- max(Diag.prod)

编辑：你使用数据框，但你可以使用 as.matrix(x) 将其轻松转换为矩阵。

Actually, contrary to Dirk I believe that you should be aware of vectorization in R as soon as possible. The loop structure you try to implement is far from optimal, and actually redundant. Using a for-loop should be done only in very specific cases. Check the discusison in this question. Take a look at the help files of convenience functions like diag(), combn(), prod() and apply().

It's easy to combine them to do what you want :

x <-matrix(1:400,ncol=20)

Diag <- diag(x)

Id <- combn(1:length(Diag),2)

Diag.prod <- apply(matrix(Diag[Id],ncol=2),1,prod)

Max.Diag.prod <- max(Diag.prod)

Edit : You use a data frame, but you can use as.matrix(x) to convert this easily to a matrix.

回复收藏 0 原文

~没有更多了~