mixed-models

    5热度

    2回答

    转换混合模型公式从SAS至R

    重复测量有几个问题和有关混合模式为更复杂的实验设计的职位，所以我想这更简单的模型将帮助其他新手则在这个过程中，以及一 所以，我的问题是我想从SAS制定中的R重复测量的协方差分析PROC MIXED程序： proc mixed data=df1; FitStatistics=akaike class GROUP person day; model Y = GROUP X1/solution a
    r sas mixed-models 2012-08-05

    22热度

    5回答

    如何获得混合效应模型中的系数及其置信区间？

    在lm和glm模型，我使用的功能coef和confint要实现的目标： m = lm(resp ~ 0 + var1 + var1:var2) # var1 categorical, var2 continuous coef(m) confint(m) 现在我添加到模型中随机效应 - 利用lme4包lmer功能使用混合效应模型。但是，功能coef和confint不再适用于我！ > mix1
    r lme4 random-effects mixed-models 2012-06-17

    3热度

    2回答

    入门广义最小二乘均值为NLME或lme4

    最小二乘与他们的标准误差固定效应的aov对象可以用model.tables函数来获得： npk.aov <- aov(yield ~ block + N*P*K, npk) model.tables(npk.aov, "means", se = TRUE) 我不知道如何获得广义最小二乘与手段从nlme或lme4对象的标准误差： library(nlme) data(Machines) f
    r lme4 mixed-models nlme lsmeans 2011-12-02

    2热度

    1回答

    如何制作lmer模型obejcts的列表以用于R中的for循环？

    我试图在R（我的第一个！）中编写一个for循环，以生成并保存使用函数lmer在包lme4中拟合的几个混合效果模型的诊断图。这是我做了什么至今所sleepstudy数据为例： require(lme4) mod1<-lmer(Reaction ~ Days + (1|Subject),sleepstudy) mod2<-lmer(Reaction ~ 1 + (1|Subject),sleep
    r list for-loop lme4 mixed-models 2013-02-22

    4热度

    3回答

    R：与GLMM的连续和分类变量的交互作用图（lme4）

    我想制作交互图以直观显示分类变量（4个级别）和标准化连续变量来自回归模型的结果。 with(GLMModel, interaction.plot(continuous.var, categorical.var, response.var)) 是不是我正在寻找。它会生成一个曲线，其中斜率随连续变量的每个值而变化。我期待，使具有恒定的斜坡曲线如下面的情节： 任何想法？ 我适合的形式fit<-glmer
    r plot interaction lme4 mixed-models 2012-05-05

    2热度

    1回答

    重复测量的混合效果模型

    我正在尝试在重复测量的数据集上开发混合效果模型。 Met上提交给3个处理24个样品的一系列随机选择的天测量（Treat，具有水平c，uc和ga）。 Met由于日间天气条件的差异而导致的水平变化（Date）。日期因此成为模型的第二个随机效应（以及取样项目（ID））。 我的主要兴趣是看看Treat在几天内对Met是否有显着影响。 一些样本数据： # create example data frame
    r mixed-models nlme 2013-01-09

    1热度

    1回答

    从lme摘录N摘要

    我正在寻找一种快速方法来提取在lme测试摘要中生成的观察值数量（N），而不是复制和粘贴？我发现this在线，但该命令似乎不是内置的，并且texreg程序包未上传到此站点。 n <- nobs(model) #extract number of observations 非常感谢。
    r summary mixed-models 2012-10-04

    21热度

    1回答

    提取预测波段适合

    我有以下型号 x <- rep(seq(0, 100, by=1), 10) y <- 15 + 2*rnorm(1010, 10, 4)*x + rnorm(1010, 20, 100) id <- NULL for(i in 1:10){ id <- c(id, rep(i,101)) } dtfr <- data.frame(x=x,y=y, id=id) library(nlme
    r prediction confidence-interval mixed-models 2013-01-16

    -1热度

    1回答

    如何为混合模型生成残差与预测变量的关系图？

    我的混合模型如下： model <- lme(Cost~1+Units, random=~1+Units|Factory, method="ML", data=A) 我被告知申请下面的代码绘制残差与拟合值和它的工作： plot(fitted(model), resid(model)) 但是，它给我看时，一个错误信息我试图绘制残差与单位，预测变量。 > plot(Units(model),
    r statistics hierarchical-clustering mixed-models 2013-03-08

    1热度

    1回答

    添加abline到augPred情节

    道歉什么可能是一个很基本的问题，我是很新的R. 我期待以平均出的值提供给读出我的augPred情节一段时间之间的预测。 > head(tthm.groupeddata) Grouped Data: TTHM ~ Yearmon | WSZ_Code WSZ_Code Treatment_Code Year Month TTHM CL2_FREE BrO3 Colour PH TURB
    r multi-level mixed-models 2013-05-04
最新问题