如何使用黄土gamm拟合模型

Question

我想用黄土伽玛给出方差的想法。但是用我的数据拟合模型会给我一个错误。

我有一个昆虫数量和处理的响应变量，以四个级别作为解释变量。该计数是随时间收集的。

loess_model <- loess(ground_beetle ~ Treatment, span=0.4, data=m_data)
gg<-gamm(ground_beetle ~ lo(Treatment) + s(Sampling_date, k = 5, fx = 
T), family = poisson, data = m_data)

Error in loess(ground_beetle ~ Treatment, span = 0.4, data = m_data) : 
predictors must all be numeric
Error in names(dat) <- object$term :'names' attribute [1] must be the 
same length as the vector [0]

Answer 1

您不能使用gamm()，也不能使用 mgcv 软件包中的其他任何东西。 lo()来自 gam 软件包，这些平滑不能与来自 mgcv 的受罚样条线思想结合使用。

更常见的问题和loess()错误的直接原因是您可以平滑因子变量，例如Treatment。这两个概念从根本上是不兼容的。