nlme

4热度

1回答

我通过RPy与nlme和lme4 R函数接口，并且我想从我的python控制台访问输出摘要。我运行下面的代码： test1=nlme.lme(r.formula('Pupil~CoI*Time'), random=r.formula('~1|ID'),data=dfr) test2=nlme.lme(r.formula('Pupil~CoI*measurement'),random=r.for

0热度

1回答

lme - 错误：类型'closure'的对象不是子集合

我想运行线性混合效果模型（包nlme），但我重复得到错误：类型'closure'的对象不是子集。 > apoeht <- read.csv("apoeht.csv") > library(nlme) > model.a <- lme(Timmrec ~ age, data = apoeht, random = ~ age | pathid, + na.exclude) Error:

3热度

2回答

打印适合lmList的所有模型的R平方

我使用lmList来适应480关系，我希望每个这些的R2。下面是一个例子数据集和模型，它是非常接近其真正的样子，但我有480欧（实验单位）： eu mass day 11 .02 1 11 .03 2 11 .04 3 11 .06 4 12 .01 1 12 .03 2 12 .04 3 12 .05 4 fit<-lmList(mass ~ day | eu, data=

6热度

1回答

nlme fit：vcov与汇总

我已经安装了package nlme的nlme()模型。现在我想模拟一些预测区间，考虑到参数的不确定性。为此，我需要提取固定效应的方差矩阵。据我所知，这样做有两种方式： vcov(fit) 和 summary(fit)$varFix 这两个产生相同的矩阵。不过，如果我检查 diag(vcov(fit))^.5 这是不一样的summary(fit) 报道的标准误我错了期待这两个是

2热度

1回答

如何在lme4或nlme中编码对比度？

我正在寻找使用lme4时在R中运行对比度的最有效方法。我一直在和一位我真正信任的统计顾问合作，她给了我下面的代码。我在6种治疗方法之间进行了对比，我在6个不同年份运行这些对比。所以我最终写出了90个对比。现在我将在模型中包含另一个因素（抽样深度），这将导致我写450个对比度。必须有更好的方法吗？我一直在阅读在R中运行对比度的方法，但与lme4没有太大的关系。 nlme也适用于我，但它也不清楚它

7热度

1回答

如何在lme中分别放置两个随机效果？

我正在做nlme包装中REML适合的线性混合效果模型。而这些都是为我工作代码： # Linear mixed-effects model fit by REML (intercept and not slope) x <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker) summary(x) # Linear

2热度

1回答

如何为视觉表示绘制二元混合效果模型

我想绘制一个二进制混合效果模型的结果在视觉表示中的一篇论文。我使用LME以适应混合模型： M2 <- lme(Pass ~ zone.time + length + Fat, random =~ 1 | Year) 通行证=二进制1/0 zone.time，长度&脂肪=连续得到： Linear mixed-effects model fit by maximum likelih

0热度

1回答

nlme：不正确的因子名称与NA？

基本上，我正在运行nlme::lme，缺少一个值。该模型与na.action=na.omit吻合良好，但拟合/残差/ coef的名称似乎都是连续改变的？ ## Generate data --------------------- X1=gl(2,4) X2=gl(2,2,8) Y=rnorm(8) dat=data.frame(Y=Y,X1=X1,X2=X2) dat ## mis

0热度

1回答

使面板填满格子？

我正在绘制一些来自nlme的分组数据，我有120个面板。默认绘图plot(dataG)将它们放在2行60列中，填满屏幕，但难以阅读。当我指定布局plot(dataG), layout= c(12,10))时，我得到了正确数量的行和列，但列全部都是一起擦亮的。我不确定问题是否因为一切都发生在nlme内部，但我还没有在我的nlme书中找到解决方案。您可以在这里找到 https://www.drop

0热度

1回答

两个随机条款与NLME

我执行与nlme包河混合模式我的情况是：混合模式是： MY = DFC + DFC2, random=~DFC|Animal, data=my_data) 其中Animal是随机效应。然而，如果我写这样的模型，我只能获得随机截取，并且倾斜DFC（Animal），但不是DFC2。我想也有DFC2随机斜率（Animal）！你能帮我吗？非常感谢你，