如何使用loess拟合模型

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。这两个概念从根本上是不相容的。