在r中运行for循环时设置矩阵单元格的最小值和最大值

Question

我在一个for循环中运行代码，对矩阵中的单个单元格进行一些简单的算术。矩阵的所有单元格都应该有0到1之间的值。矩阵中的值永远不应该变成负值或超过1。我无法使用传统的“最小”和“最大”限制设置来解决这个问题。

示例数据：

x<-1:5
y<-runif(length(x)*length(x),0,1)
M<-matrix(y, length(x), length(x))

例如，从每个单元格中减去0.01（在我的实际情况中，在循环的每次迭代中减去的0.01来自特定的单元格，而不是这里所写的每个单元格，但我认为这不应该影响如何设置最小/最大值？

for (i in 1:50){
M[,]<-M[,]-0.01
}

e、 g.或者如果向矩阵的每个单元格添加0.01：

for (i in 1:50){
M[,]<-M[,]+0.01
}

如果在循环之后返回矩阵，矩阵中的值将不会像我希望的那样包含在0和1范围内。我对如何做到这一点感到困惑。

如果我在做某件特定于矩阵中某个“单元”的事情，例如。

for (i in 1:100){
M[3,4]<-M[3,4]-0.01
}

这会改变设置最小/最大值的方式吗？

Answer 1

剪辑功能怎么样

clip<-function(x, min=0, max=1) {
    x[x<min]<-min; 
    x[x>max]<-max; 
    x
}

这将设置所有值

 clip(M-.9)

添加也一样

 clip(M+.9)

还是你试图减去最小值的问题？当你做这样的事情时

for (i in 1:50){
    M[,]<-M[,]-0.01
}

循环中的语句每次都对向量的每个元素进行操作。所以你真的减去了0.5。这可能太多了。

Answer 2

我认为有两种方法可以做到这一点：

运行循环并获得结果（超出范围），然后使用以下命令强制数据在最小/最大值范围内：

M[which(M<min)] = min
M[which(M>max)] = max

你可以把赋值函数修改为

for (i in 1:50){
    M[which(M<max)] <- M[which(M<max)]+0.01
}

对于添加案例和

for (i in 1:50){
    M[which(M>min)] <- M[which(M>min)]-0.01
}

对于减法的情况。