我正在尝试运行flux.package来计算一些土壤气体通量,并且基本上遵循示例脚本和格式。
但是,这样做之后:
tt.neetest = read.csv("tt.neetest.csv", header = TRUE)
## prepare flux estimation
# make parts with chop
tt.parts <- chop(tt.neetest, factors=c("session", "spot","plot"),
nmes=c("spot", "date", "session","plot"), min.cm=10)
# prepare var.par list (like with flux)
vp <- list(CO2 = "NEE", time = "datetime", area = "area",
volume = "volume", t.air = "t.cham", p.air = 101325)
一切正常,我尝试使用以下方法计算通量:
# do the flux estimation
tt.flux <- fluxx(tt.parts, subset=c(1:30), vp, pdk=0.5,
range.lim=3, out.unit="mg")
这给了我以下错误消息:
Error in runsd(x$ghg, ceiling(n * wndw)) : 'k' must be larger than 2
那是什么?以前从未见过这样的错误。什么是“ k”?我的数据中没有这样的数据。
干杯,克劳迪娅