R |mixed effect model

最新推荐文章于 2023-10-18 15:39:47 发布

gavin_cdc

最新推荐文章于 2023-10-18 15:39:47 发布

阅读量1.9k

点赞数

分类专栏： lme4 R语言

本文链接：https://blog.csdn.net/gavin_cdc/article/details/99682309

版权

本文深入探讨了R语言中如何使用混合效应模型处理数据，包括模型构建、参数估计及结果解释，阐述了该方法在处理重复测量和相关数据中的优势。

摘要由CSDN通过智能技术生成

---
title: "mixed effect model"
author: "gavin"
date: "2019年8月16日"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

```{r}
library(lme4)
politeness=read.csv("/Users/Gavin/Downloads/politeness_data.csv")
# head(politeness), tail(politeness),summary(politeness), str(politeness), colnames(politeness)
which(is.na(politeness)==T)
boxplot(frequency ~ attitude*gender,
        col=c("white","lightgray"),politeness)
```

```{r}
NEWdata <- data.frame(t(politeness))
NEWdata <- NEWdata[,-39]
NEW_MODEL = lmer(as.numeric(NEWdata[5,]) ~ as.numeric(NEWdata[6,]) + as.numeric((1|NEWdata[7,])))
NEW_MODEL
```


```{r}
politeness.model = lmer(frequency ~ attitude +
                          (1|subject) + (1|scenario), data=politeness)
politeness.model

pol