当前,我正尝试在代码中存储每次循环的纬度矩阵。
当前,我正在尝试按以下方式存储变量:
latitude[p] <- lat
其中纬度是我要存储的变量,p是循环号,而lat是随机生成的纬度矩阵。
以下是我的可复制代码:
p=0
latitude = 0
while(p<5)
{
p = p + 1
lat <- round(rnorm(n=50, m = 55, sd = 8),0)
lon <- round(rnorm(n=50, m = 76, sd = 20),0)
no2 <- round(rnorm(n=100, m = .3, sd = .1),2)
partial_data <- data.frame(lat,lon,no2)
coordinates(partial_data) = ~lat+lon
spplot(partial_data,'no2')
latitude[p] <- lat
}
当前,我收到错误消息:
In latitude[p] <- lat :
number of items to replace is not a multiple of replacement length
,当我调用“ latitude [1]”时,得到一个数字作为输出。我想要为循环生成的完整随机矩阵作为输出。
谢谢!