我有一个数据dat，该数据具有15000个响应变量，并且在多个组和年份中收集了10个预测变量。
我使用lmer在R中拟合了一个模型，现在我想使用挖泥函数创建多个模型，并使用MuMIn软件包进行模型平均。这是我的方法。
global.model <- lmer(y ~ x1 + I(x1^2) + x2 + I(x2^2) + x3 + I(x3^2) + x4 + I(x4^2) + x5 + I(x5^2) + x6 + x7 + x8 + x9 + x10 + (time|group), REML = F, data = dat)
为了避免计算时间，我想设置一些条件，以使模型中的所有平方项仅在存在线性项时才会出现（例如I(x2^2)仅在{{1 }} 存在）。我还制定了一些规则，下面的评论对此进行了解释
x2我想问一下上面的语法是否正确，以指定我想要的不同规则。