Multilevel Modeling Using R 第三章

您所在的位置:网站首页 r语言lmer函数 Multilevel Modeling Using R 第三章

Multilevel Modeling Using R 第三章

2023-12-20 11:27| 来源: 网络整理| 查看: 265

从这一章开始我们就来介绍下用R处理双层次的模型

R包

首先介绍一下在多层次模型中,可以用R包nlme的nlme()函数来实现建模,或者R包lme4的lmer()函数,接下来,我们就看一下这两个函数的用法:

lme(fixed, data, random, correlation, weights, subset, method, na.action, control, contrasts = NULL, keep.data = TRUE) lmer(formula, data, family = NULL, REML = TRUE, control = list(), start = NULL, verbose = FALSE, doFit = TRUE, subset, weights, na.action, offset, contrasts = NULL, model = TRUE, x = TRUE,...)

下面我们就来看一下怎么具体使用这些函数

nlme包 1. 简单多层次模型

回顾下第二章的固定效应和随机效应

比如这个模型,β0j代表截距,属于随机效应,那么β0j即为随机截距,因为每一个类的截距都不一样;β1j代表斜率,也属于随机效应,不过是β1j = γ10(每个类斜率的均值),所以β1j是固定的 其中Var(εij) = σ^2组内方差;Var(U0j) = τ^2 组间方差 此例子的固定效应是主要研究阅读分数x与词汇分数y的线性关系;而随机效应则是对Level1的线性回归系数的影响

我们以第二章的例子,以学校聚类,阅读分数与词汇分数的相关性的例子为例

Model3.0


【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3