Question

这个问题，可能有一个简单的解决方案，但我想不出怎么做......

所以我有一个如下脚本：

# ------------------ MODEL SETUP ----------------------------------------# simulation length
t_max <- 50

# arena
arena_x <- 100
arena_y <- 100

# plant parameters

a <- 0.1
b <- 0.1
g <- 1

# list of plant locations and initial sizes
nplants <-dim(plantLocsX)[1]*dim(plantLocsX)[2]

iterations<-5
totalBiomass<-matrix(0,nrow=iterations,ncol=1)

# starting loop 

sep <- 10


# Original matrix
plantLocsX <- matrix(rep(seq(0,arena_x,sep), arena_y/sep),
                    nrow=1+arena_x/sep,
                     ncol=1+arena_y/sep)


plantLocsY <- t(plantLocsX)

plantSizes <- matrix(1,nrow=nplants,ncol=1)

# Plot the plants
radius <- sqrt( plantSizes/ pi )
symbols(plantLocsX, plantLocsY, radius, xlim = c(0,100), ylim=c(0,100), inches=0.05, fg = "green", 
        xlab = "x domain (m)", ylab = "y domain (m)", main = "Random Plant Locations", col.main = 51)


# Calculate distances between EACH POSSIBLE PAIR of plants
distances <- matrix(0,nrow=nplants,ncol=nplants)

for (i in 1:nplants){
  for (j in 1:nplants){
    distances[i,j] <- sqrt( (plantLocsX[i]-plantLocsX[j])^2 + (plantLocsY[i]-plantLocsY[j])^2 )
    }
}


# ------------------ MODEL RUNNING ---------------------------------------

我需要更改包含植物位置和植物大小的数组，以便在每个时间步骤中删除并添加条目（分别模拟死亡率/繁殖）。每次迭代后，“距离”必须用工厂位置和大小更新......我只能想到复杂的方法：在每个时间步骤破坏和构造新的矩阵以适应新的元素数量但必须有函数使这更简单....任何建议？非常感谢！！

更改数组以在R中的每个时间步骤添加/删除条目

0 个答案: