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)

和类似的尝试是我认为的问题。

来自?lme：

na.action：表明数据发生时应该发生什么的功能包含“ NA”

您提供的是数据，而不是函数，因为na.omit(dataset)返回一个data.frame，其中NA包含已删除的行，而不是可应用于指定的data=的行。只是：

na.action=na.omit

或类似的na.*函数就足够了。

确定此类问题的一种方法是肯定使用?debug-debug(lme)，然后逐行浏览该函数，以准确了解错误是如何响应的。

'lme'错误R“尝试应用非功能

1 个答案: