''lme'error R“试图应用非功能

Question

我正在使用以下代码对我的数据集进行lme分析

M1 <- lme(VT ~ visit + sx + agevis + c_bmi + gpa + qa + BP + MH + ethn, data = Cleaned_data4t300919, random = ~ 1 + visit |id, corAR1(),method = "ML", na.action = na.omit(Cleaned_data4t300919))

并且我收到以下错误消息：

model.frame.default（formula =〜visit + sx + agevis + c_bmi +：尝试应用非功能

我不确定我做错了什么或如何运行模型。我真的很感谢一个答案。谢谢。

[我正在尝试以VT作为我的因变量，以visit作为我的时间变量来运行线性混合效应模型，该模型具有一阶自回归相关性，即ML估计器，用于缺少观测值的数据。

我曾尝试通过以下方式更改代码，但收到相同的错误消息

library(nlme)
?lme
fm2 <- lme(VT ~ visit + sx + agevis + c_bmi + gpa + qa + BP + MH + ethn, data = Cleaned_data4t300919, random = ~ 1|id, corAR1(),method = "ML", na.action = na.pass(Cleaned_data4t300919))

fm2 <- lme(VT ~ visit + sx + agevis + c_bmi + gpa + qa + BP + sfnMH + ethn, data = Cleaned_data4t300919, random = ~ 1 + visit |cenid, corAR1(),method = "ML", na.action = na.omit(Cleaned_data4t300919))

fm2 <- lme(VT ~ visit + sx + agevis , data = Cleaned_data4t300919, random = ~ 1 + visit |id, corAR1(),method = "ML", na.action = na.omit(Cleaned_data4t300919))

fm2 <- lme(VT~visit + sx + agevis + c_bmi + gpa + qa + BP + MH + ethn, data = Cleaned_data4t300919, na.action = na.exclude(Cleaned_data4t300919))

fm2 <- lme(formula= sfnVT ~ visit + sx + agevis , data = Cleaned_data4t300919, random = ~ 1 + visit |cenid, corAR1(),method = "ML", na.action = na.omit(Cleaned_data4t300919))

我想获得代码的估算值并使用ggplot绘制估算值。

Answer 1

na.action = na.omit(Cleaned_data4t300919)