r-RothC模型描述和soiR实现之间的参数计算差异

Question

我正在尝试使用SoilR包运行RothC模型。其中一个步骤涉及估计水分对分解速率的影响。

根据模型描述(第14页(，参数b计算为：

如果符合TSMD<0.444 max.TSMD，
b = 1.0
否则，
b = 0.2 + (0.8) * (max.TSMD - acc.TSMD)/(max.TSMD - 0.444 max.TSMD)

但是soilR实现，fW.RothC.R，是：

b=ifelse(Acc.TSMD > 0.444*Max.TSMD, 1, (0.2+0.8*((Max.TSMD-Acc.TSMD)/(Max.TSMD-0.444*Max.TSMD))))

正如您所看到的，该条件具有>而不是<

如何解释？

我没有足够的信誉来创建标签，但添加标签soilr和rothc会很有用。

Answer 1

R函数是可以的。RothC的手册有一个错误(或者TSMD句子后面可能有一个不可见的字符减号(。检查RothC的另一个实现(solassestination(，有相同的公式。https://search.r-project.org/CRAN/refmans/soilassessment/html/carbonTurnover.html