楼主: bfzldh
57395 166

[丢掉HLM软件]基于R的多层模型分析教程贴   [推广有奖]

贵宾

院士

78%

还不是VIP/贵宾

-

威望
0
论坛币
368177 个
通用积分
5087.1285
学术水平
334 点
热心指数
381 点
信用等级
274 点
经验
1177 点
帖子
3218
精华
0
在线时间
2774 小时
注册时间
2013-5-30
最后登录
2024-10-12

初级热心勋章 中级热心勋章

楼主
bfzldh 学生认证  发表于 2016-9-29 14:48:17 |只看作者 |坛友微信交流群|倒序 |AI写论文
相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
         还在为找不到破解版HLM软件而头疼?为什么还不放弃HLM软件?尽管HLM软件的图形化操作似乎简单一点,但是
        一、HLM软件需要对多层数据分别建立文件,这种操作方式与我们的一般习惯不符。
        二、HLM软件只能处理一个因变量;
        三、HLM软件不能对中介效应直接进行检验。
        ……
        Mplus和R都可以用来处理多层模型,且各具千秋。本贴将不定时更新,逐步阐述基于R的多层模型分析。更新贴已全部置顶,可点击标题上方的“查看完整内容”定位相应帖子。本贴所介绍的内容主要来自于下面这本书


Finch, W. H., Bolin, J. E., & Kelley, K. (2014). Multilevel modeling using R. Crc Press.

        首先将这本书所涉及的数据奉上
         

本帖隐藏的内容

datasets.zip (1.68 MB) 本附件包括:
  • Achieve.csv
  • Achieve.sav
  • Babies.csv
  • Babies.sav
  • Cassady.csv
  • Cassady.sav
  • Cooking.csv
  • Cooking.sav
  • CoronaryArtery.csv
  • CoronaryArtery.sav
  • Language.csv
  • Language.sav
  • MathFinal.csv
  • MathFinal.sav
  • Politics.csv
  • Politics.sav
  • PrimeTime.csv
  • PrimeTime.sav
  • Rehab.csv
  • Rehab.sav



二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:HLM软件 模型分析 多层模型 HLM Multilevel 软件 破解版 因变量 模型 2014

回帖推荐

Alex_Do 发表于97楼  查看完整内容

非常好的资料,万分感谢

jiandong4388 发表于96楼  查看完整内容

顶起来

tigerlab 发表于75楼  查看完整内容

作者您好,学习了好多,谢谢,非常好的资料。 补充一个,今天刚学的,调用random effects的时候,intervals()只有在nlem中管用,在lme4中可以用rand(). 我还有个问题想问作者,如果想做一个三重交互的简单效应分解,想要检验简单效应系数是否显著,作者能做吗?我的三重交互两个level 1 中的自变量,一个level2中的moderator,写出来应该是x1*x2*w. 非常感谢

bfzldh 发表于62楼  查看完整内容

13.当lmer函数使用REML算法时,如何得到拟合指数。#当lmer函数使用REML算法时,其输出结果中没有直接报告模型拟合信息。#但是,我们可以通过专门的函数来提取这些信息。#基于REML的离差 #AIC #log likelihood,基于REML的离差

bfzldh 发表于61楼  查看完整内容

12.将参数估计的算法改为ML #在nlme程序包lme函数中,需要声明method="ML" #在lme4程序包lmer函数中,需要声明REML =FALSE #使用summary得到模型估计结果,使用intervals得到随机效应的置信区间 #对于lmer函数,只有当算法改为ML之后,summary才会报告模型拟合结果 #对于lmer函数,可使用confint函数得到随机效应的置信区间

bfzldh 发表于60楼  查看完整内容

11.在lme4中将多个随机效应设定为相关的或相互独立的 #在lme4中,我们可以界定多个自变量的随机斜率之间之间是否相关。 #下面呈现了Model3.10和Model3.11的语句, #Model3.10的随机效应写在两个单独的括号里,这表示相应的两个随机效应是不相关的, #Model3.11中,随机效应写在同一个括号里,表示相应的两个随机效应是相关的 #Model3.10的结果如下 Linear mixed model fit by REML t-tests use Satterthwaite approximat ...

bfzldh 发表于59楼  查看完整内容

10.使用lme4(lmerTest)拟合随机斜率模型#下面的模型仅包含一个水平1的自变量gevocab[/backcolor],并设定其斜率在水平2自由变化,只需在括号中声明gevocab[/backcolor]#只有提前中心化,才能成功计算[/backcolor]#结果如下Linear mixed model fit by REML t-tests use Satterthwaite approximations to degrees of freedom [lmerMod]Formula: geread ~ cgevocab + (cgevocab | school) Data: Achieve REML criterion at converg ...

bfzldh 发表于58楼  查看完整内容

9.使用lme4(lmerTest)在模型中加入多个自变量,或者在模型中加入交互作用#下面的例子包含了多个自变量: #结果如下 Linear mixed model fit by REML t-tests use Satterthwaite approximations to degrees of freedom [lmerMod]Formula: geread ~ gevocab + senroll + (1 | school) Data: Achieve REML criterion at convergence: 43152.1 Scaled residuals: Min 1Q Median 3Q Max -3.0834 -0.5 ...

bfzldh 发表于57楼  查看完整内容

8.使用lme4拟合零模型 #6%使用lme4拟合零模型。 #结果如下 Linear mixed model fit by REML ['lmerMod'] Formula: geread ~ (1 | school) Data: Achieve REML criterion at convergence: 46268.3 Scaled residuals: Min 1Q Median 3Q Max -2.3229 -0.6378 -0.2138 0.2850 3.8812 Random effects: Groups Name Variance Std.Dev. school (Intercept) 0 ...

bfzldh 发表于56楼  查看完整内容

7.使用lme4程序包lmer函数拟合随机截距模型 #lme4函数界定模型的语句更简洁,对复杂模型的界定方法也更灵活。我们使用lme4程序包中的lmer函数拟合Model3.1, #不同于lme函数,lmer函数中随机效应是用括号表示出来的,而不是专门用random来声明。上述语句的计算结果如下#6%lme4 1.0以后的版本发生了变化,“[gn]lmer now produces objects of class merMod rather than class mer as before” #6%仅输入命令Model3.7,可得到以下 ...

bfzldh 发表于44楼  查看完整内容

6.随机效应的显著性检验 #使用nlme的lme函数拟合多层模型,通过summary()获得的输出结果中不包含对于随机效应的显著性的检验结果。 #结果如下 Linear mixed-effects model fit by REML Data: Achieve AIC BIC logLik 43175.57 43219.02 -21581.79 Random effects: Formula: ~1 | school (Intercept) Residual StdDev: 0.316492 1.940268 Fixed effects: geread ...

bfzldh 发表于43楼  查看完整内容

5.中心化 #中心化指的是从变量的原始分数中减去该变量的均值。 #在R中可通过创建新变量实现中心化。#首先将变量进行中心化,得到中心化后的自变量#中心化后的变量被创建以后,用中心化后的变量进行建模分析即可: #中心化并不会从本质上改变统计结果(读者可自行比较模型拟合指数),但会改变对系数的解释方式。结果如下 Linear mixed-effects model fit by REML Data: Achieve AIC BIC logLik 43155.49 43198 ...

bfzldh 发表于38楼  查看完整内容

4.使用nlme估计包含交互作用的多层模型 **** 本内容被作者隐藏 ****

bfzldh 发表于36楼  查看完整内容

3.s 去除截距的随机效应,保留斜率的随机效应

bfzldh 发表于34楼  查看完整内容

3.使用nlme拟合随机斜率模型 #在随机斜率模型中(random coeffcients model),自变量对依变量的预测力在水平2随机变化。#就当前问题来说,这意味着gevocab对geread的预测力随着水平2学校的变化而变化。#在界定随机效应时,1指的是截距,如果我仅需要随机截距模型,~1|school就足够了。#如果我们想要水平1自变量的斜率在水平2随机变化,我们需要改变相应的语句#前提是,gevocab已经被包含在固定效应部分。这里需要允许斜率和gevoc ...

bfzldh 发表于25楼  查看完整内容

2.在随机截距模型中加入水平1的自变量,使用AIC,BIC比较模型,并计算R-square #词汇成绩(gevocab)可能对阅读成绩具有预测作用,要将gevocab作为自变量,则相应的R语句为 Model3.1

bfzldh 发表于19楼  查看完整内容

0.用R进行多层分析之前的准备工作。 **** 本内容被作者隐藏 ****

bfzldh 发表于9楼  查看完整内容

1.使用nlme建立零模型,并计算ICC **** 本内容被作者隐藏 ****
已有 6 人评分经验 论坛币 学术水平 热心指数 信用等级 收起 理由
Lyndon-he + 2 + 2 + 2 精彩帖子
jgzjpart + 1 + 1 + 1 精彩帖子
as463632652 + 5 + 5 + 5 精彩帖子
我的素质低 + 100 + 100 + 5 + 5 + 5 精彩帖子
jiandong4388 + 100 精彩帖子
xddlovejiao1314 + 100 + 5 + 5 + 5 精彩帖子

总评分: 经验 + 300  论坛币 + 100  学术水平 + 18  热心指数 + 18  信用等级 + 18   查看全部评分

本帖被以下文库推荐

沙发
xingzhe1204 在职认证  企业认证  发表于 2016-9-29 15:19:01 |只看作者 |坛友微信交流群
靠谱,欢迎多分享,点赞哦
已有 1 人评分论坛币 收起 理由
bfzldh + 6 鼓励积极发帖讨论

总评分: 论坛币 + 6   查看全部评分

使用道具

藤椅
2行者8805 在职认证  发表于 2016-9-29 16:49:19 |只看作者 |坛友微信交流群
非常感谢赛
已有 1 人评分论坛币 收起 理由
bfzldh + 6 鼓励积极发帖讨论

总评分: 论坛币 + 6   查看全部评分

使用道具

板凳
叶孤傲 发表于 2016-9-29 17:07:38 |只看作者 |坛友微信交流群
可以弱弱的问一句:这个基于R的多层模型分析一般用于解决什么问题啊
已有 1 人评分论坛币 收起 理由
bfzldh + 6 多层模型用于处理嵌套数据,例如,学生嵌套.

总评分: 论坛币 + 6   查看全部评分

使用道具

报纸
xddlovejiao1314 学生认证  发表于 2016-9-29 17:28:08 |只看作者 |坛友微信交流群
虽然不用R,支持。一般用stata实现分层线性模型。
已有 1 人评分论坛币 收起 理由
bfzldh + 6 好的意见建议

总评分: 论坛币 + 6   查看全部评分

使用道具

地板
jiandong4388 学生认证  发表于 2016-9-29 17:46:39 |只看作者 |坛友微信交流群
支持支持
已有 1 人评分论坛币 收起 理由
bfzldh + 6 鼓励积极发帖讨论

总评分: 论坛币 + 6   查看全部评分

使用道具

7
孙继伟 发表于 2016-9-29 18:49:53 |只看作者 |坛友微信交流群
支持支持,需要多学习
已有 1 人评分论坛币 收起 理由
bfzldh + 6 鼓励积极发帖讨论

总评分: 论坛币 + 6   查看全部评分

使用道具

8
zlgsx 学生认证  发表于 2016-10-1 20:46:31 |只看作者 |坛友微信交流群
听起来不错的东西,关注一下。
已有 1 人评分论坛币 收起 理由
bfzldh + 6 鼓励积极发帖讨论

总评分: 论坛币 + 6   查看全部评分

使用道具

9
bfzldh 学生认证  发表于 2016-10-5 13:01:07 |只看作者 |坛友微信交流群
1.使用nlme建立零模型,并计算ICC

本帖隐藏的内容


#研究者想要知道词汇成绩(gevocab)对一般阅读成绩(geread)的预测能力。
#由于学生嵌套于学校,标准的线性回归方程是不合适的。
#我们将要建立的第一个模型是零模型,即,模型中不包含自变量。
#该模型可用来估计残差和截距的变异(仅考虑嵌套于学校)。相应的lme语句见下方。

Model3.0 <- lme(fixed = geread~1, random = ~1|school, data =Achieve)

#如果仅输入Model3.0只能得到粗略的结果,我们可通过summary(Model3.0)获得计算结果。

summary(Model3.0)

#结果如下,下面这句话说明,使用的是REML算法
#Linear mixed-effects model fit by REML
#数据为Achieve
# Data: Achieve
#模型拟合指数
#      AIC      BIC    logLik
#  46274.31 46296.03 -23134.15
#随机效应部分
#Random effects:
# Formula: ~1 | school
#        (Intercept) Residual
#StdDev:   0.6257119  2.24611
#固定效应部分
#Fixed effects: geread ~ 1
#               Value  Std.Error    DF t-value p-value
#(Intercept) 4.306753 0.05497501 10160 78.3402       0
#标准化的组内残差,不重要
#Standardized Within-Group Residuals:
#       Min         Q1        Med         Q3        Max
#-2.3229469 -0.6377948 -0.2137753  0.2849664  3.8811630
#样本量
#Number of Observations: 10320
#cluster的数量
#Number of Groups: 160
#根据ICC的公式,ICC=二层方差/(二层方差+一层方差),结合上述随机效应部分,可计算
#ICC=0.6257119*0.6257119/(0.6257119*0.6257119+2.24611*2.24611)=0.0720157351971356

已有 2 人评分经验 论坛币 学术水平 热心指数 信用等级 收起 理由
bdim + 1 + 1 精彩帖子
我的素质低 + 100 + 100 + 5 + 5 + 5 精彩帖子

总评分: 经验 + 100  论坛币 + 100  学术水平 + 6  热心指数 + 6  信用等级 + 5   查看全部评分

使用道具

10
maayunnbiaoo 学生认证  发表于 2016-10-5 14:53:28 |只看作者 |坛友微信交流群

支持支持,需要多学习
已有 1 人评分论坛币 收起 理由
bfzldh + 6 鼓励积极发帖讨论

总评分: 论坛币 + 6   查看全部评分

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-11-5 12:47