我正在使用素食包中的varpart函数。这是我的脚本：
CR_all.fun_var_part <- varpart(all_fungi, ~ age, ~ lat +long, data = Cr_var)
CR_all.fun_var_part
plot(CR_all.fun_var_part, digits=2)
我的反应变量，所有的真菌，都是我用hellinger方法转化的丰度矩阵。 年龄（X1）代表我的森林的不同年龄，我用尺度函数来缩放。 如果我想用我的GPS坐标作为X2我应该如何对待它们？我应该改变它们吗？ varpart函数的脚本应该如何？ 我希望我的问题不是太模糊，我感谢大家的帮助。
答案 0 :(得分：0)
~ lat + long与拟合线性趋势曲面相同。它只考虑线性地理趋势。如果您只对线性趋势感兴趣，那就没关系。如果您认为需要更复杂的趋势，那么您需要更复杂的模型：坐标上的多项式，坐标上的样条或（与许多人一样）PCNM变量。坐标转换会在对原始坐标绘制时改变响应的形状，但很难说如何转换。
可以使用base函数poly(lon, 2)找到多项式（对于二次多项式）。 PCNM（邻域矩阵的主坐标）可以在许多包中找到，包括素食主义者（pcnm）。
这个答案只关注如何来做到这一点。我不想赞同这些做法。