回答情况：
提问时间：
问题标签：
<- gl(n = 5, k = 10, length = 180)我的目标是找到正确的lme4::lmer() (和一个nlme+\gamma_{11}(得分)+ u_1$
由于一些未知的原因，我很难将普通方程转换为lme4/nlme
浏览 2提问于2018-01-26得票数 0
回答已采纳
也就是说，我想在oats.lmer中拟合模型(参见下面)，但使用来自nlme包的函数lme。代码是require("lme4")
oats.lmer <- lmer(yield~nitro + (1|Block:Variety), data = Oats)
浏览 4修改于2016-12-13得票数 2
回答已采纳
1回答
我有一个有关联的NLME对象。我想知道如何从模型中提取关联的范围。我正在运行模拟，所以我不能只是阅读摘要，并得到它的手动。我的模型看起来是这样的：library(nlme)
lme(fixed=temp ~ time, random=~1|day,correlation=corExp(form
浏览 1修改于2014-05-31得票数 3
回答已采纳
我已经安装了一个模型使用nlme()从package nlme。为此，我需要提取固定效果的方差矩阵。编辑:下面是一个代码示例a=c(.5,1.5)t=c(pts,pts)d=data.frame(t,y,batch)
nl
浏览 3修改于2014-05-08得票数 6
回答已采纳
我尝试使用stargazer来格式化nlme包中lmList函数的输出。如果我手动索引每个列表元素，就可以让stargazer工作，但如果我只传递完整的列表给stargazer，就不行了。library(nlme)data("iris")
m <- lmList(Sepal.Length ~ Sepal.Width | Species, data
浏览 1修改于2018-11-18得票数 1
在nlme的正畸数据集中，共有27名受试者，每个受试者在4个不同的年龄阶段进行测量。我希望利用这些数据来探讨在什么情况下模型会被高估。以下是模型：library(lme4)
m2Orthodont )
m1nlme
浏览 11提问于2022-08-16得票数 1
回答已采纳
library(nlme) data = Loblolly,1, start = c(Asym = 103, R0 = -8.5, lrc = -3.3))
我正在用R中的nlme
浏览 1提问于2017-01-16得票数 1
回答已采纳
当我在nlme中拟合数据时，我从来不会在第一次尝试中就成功，在nlme(fit.model)之后，我习惯于看到这样的情况：fit.nlme.bgf<-nlme(fit.beta, random=list(w.max + t.e + t.m ~1))
fit.nlme.bgf<-nlme(fit.beta, rando
浏览 2修改于2014-04-16得票数 3
我试图拟合一个交叉的非线性随机效应模型，作为线性随机效应模型，，如提到的，以及在这个中使用nlme软件包。不过，不管我尝试了什么，我都会犯错误。下面是一个例子# simulate datana <- 30sigma_bVariance StdDev#R b1 0.221 0.4701
# now try to do the same with `
浏览 0修改于2018-08-10得票数 7
回答已采纳
1回答
我试着使用r中的nlme包来做一个多层线性模型。
我将产量作为响应变量，降雨量作为预测变量，对6个不同的地点(州)进行了60年的预测。我正在试着看看降雨对所有地区的产量的影响程度是相同的还是不同的。
浏览 2提问于2014-06-05得票数 0
library(nlme) data = Loblolly,LME step: Loglik: -312.2787, nlminb iterations: 23 Seed Error in nlme.formulaage, Asym, R0, lrc), data = Loblolly, :
Singularity in backs
浏览 3修改于2021-05-22得票数 19
回答已采纳
如果使用"nlme“而不是"nlmer”，我如何才能适应以下模型？
数据(在文章的末尾，您可以找到数据来再现这里编写的代码)。
浏览 2修改于2018-10-02得票数 0
1回答
在其他人的帮助下，我使用R中的nlme包创建了以下logistic增长模型，它依赖于四个模型参数。天仅仅是进行测量的日子。MaxPop =每个给定站点的最大人口(覆盖率百分比在0-100%之间)。Global_Model <- nlme(PercentCover ~ I(100 *((inv.logit(exp(MaxPop) *Days^Hill/(exp(Days50)^Hill + Days问题：一个人会怎么做？我非常希望得到一些建议。
干杯!
浏览 2修改于2017-05-23得票数 1
回答已采纳
library(nlme) structure(list(class = c(91L, 91L, 91L, 91L, 91L, 91L, 92L, 92L,) model = nlme(outcome ~ exp(beta1) * s * days / (s * days + exp(1 - beta3/2) :
the leading minor of order 2 is not positive definite In addition:
浏览 0提问于2017-08-09得票数 5
我正在尝试执行一个具有混合模型效果和响应变量变化的循环。我来自和。我应该说，我尝试过创建一个函数，然后应用或应用(成功)。data<- structure(list(paciente = structure(c(6134, 6099, 6457, 6164, 6127, 6507, 6513, 6486, 6080, 6101, 6007, 6023, 6516, 6001
浏览 6修改于2022-03-04得票数 0
我正在尝试运行gls函数，但一直收到以下错误消息：我的数据如下所示：Achyranthes_aspera 37 117 0Alternanthera_adscendens -999 -999 -999
浏览 0提问于2016-10-04得票数 2
我想发布一个由nlme生成的augPred图，所以它需要看起来更好看一些(比如评论员)。我想重新排序面板和更改标题。我试着用我发现的点阵信息，但不知怎么的，它在nlme中不起作用。这是最初的数字：我用这个dataset 尝试的代码，您将需要这个函数：grow<-read.table("cobsgddv8.txt", header=T)mass ~ SSbgf(gdd, w.max, t.e, t.m), data = grow10G)
plot(intervals(fit.
浏览 1修改于2013-09-17得票数 1
回答已采纳
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号