如果我使用以下代码,随机效应应该是什么样的变量:
lme(minpH ~ `rumination (min/d)`*`VET EIWIT ratio`,
data = phase_data_compleet, random.effects(experiment))
如果随机效应变量是数字,字符的两倍,它会给出一个错误,也许我做错了什么。我的目的是运行一个以minpH为因变量,反刍、VET EIWIT RATIO为独立变量的混合效应模型,实验是随机效应,因为数据是在过去的几个实验中收集的。
如果我使用没有随机效应的代码,它会起作用:
lm(minpH ~ `rumination (min/d)`*`VET EIWIT ratio`, data = phase_data_compleet)
我想听听你对我做错了什么的建议。
你是说吗
lme(minpH ~ `rumination (min/d)`*`VET EIWIT ratio`,
random = ~ 1|experiment,
data = phase_data_compleet)
理想情况下,experiment
应该是一个因子,但可以想象lme
会自动为您转换它。
random.effects()
完全是转移注意力:它是ranef()
的同义词,在拟合模型后使用,以从总体估计中提取有关组水平偏差的信息。