符合条件

Question

我想用一个包含函数与矩阵mapply使用mapply使用的，如果是这样的：

dim_rf <- 10 

rf <- matrix(runif(dim_rf*dim_rf),ncol=dim_rf) 


fire_rf <- function(i,j) { 
    if(rf[i,j+1] + rf[i+1,j] + rf[i,j-1] + rf[i-1,j]>1) 
    rf[i,j] <-2 
} 

mapply(fire_rf, row(rf), col(rf) ) 

，但它给出了一个错误

Error en if (rf[i, j + 1] + rf[i + 1, j] + rf[i, j - 1] + rf[i - 1, j] > : 

参数长度为零

我怎么能修改功能，使其工作？

我用这个做火渗透的一个简单的模型，因此需要快速，因为100或1000

Answer 1

这将帮助您调试功能：

fire_rf <- function(i,j) { 
    print(paste(i, j)) 
    print(paste(rf[i,j+1], rf[i+1,j], rf[i,j-1], rf[i-1,j])) 
} 

R> mapply(fire_rf, row(rf), col(rf)) 
[1] "1 1" 
[1] "0.431833460461348 0.723794676829129 " 
[1] "2 1" 
[1] "0.61886831978336 0.764425948029384 0.486255120486021" 
[1] "3 1" 
[1] "0.98786562983878 0.416759668383747 0.723794676829129" 
[1] "4 1" 
[1] "0.516932599013671 0.66259386530146 0.764425948029384" 
[1] "5 1" 
[1] "0.770237588090822 0.844016372924671 0.416759668383747" 
[1] "6 1" 
[1] "0.519682829733938 0.273759116884321 0.66259386530146" 
[1] "7 1" 
[1] "0.0957026474643499 0.32785635185428 0.844016372924671" 
[1] "8 1" 
[1] "0.380909610772505 0.164877543691546 0.273759116884321" 
[1] "9 1" 
[1] "0.509409713326022 0.614429801469669 0.32785635185428" 
[1] "10 1" 
Error in rf[i + 1, j] : subscript out of bounds 
Calls: mapply -> <Anonymous> -> print -> paste 

你看，而你认为​​你传递的参数，有时你是不是有时候你逝去的不止是存在的。你需要考虑那些 '我 - 1'， '的J - 1' 和 'I + 1' 和 'J + 1'。如果你真的想这样做，我建议你写一个简单的循环。

Answer 2

我猜，这就是你想要的矩阵的大小应大。你的循环条件是错误的，使算法超出界限。

dim_rf <- 10 
rf <- matrix(runif(dim_rf*dim_rf),ncol=dim_rf) 
fire_rf <- function(i,j) { 
    if((rf[i,j+1] + rf[i+1,j] + rf[i,j-1] + rf[i-1,j])>1) 
     rf[i,j] <<-2 
} 
mapply(fire_rf, 2:(nrow(rf)-1), 2:(ncol(rf)-1)) 

希望这会有所帮助。 :)

PS。另请注意，我将本地分配更改为全局分配。如果这不是你想要的，你可以把它改回来。