Multilevel Modeling Using R 第三章 |
您所在的位置:网站首页 › r语言lmer函数 › Multilevel Modeling Using R 第三章 |
从这一章开始我们就来介绍下用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 |