重复测量资料多水平模型

setwd("C:/Users/86153/Desktop/RData/")
bc<-read.csv("C:/Users/86153/Desktop/EX4.csv",sep=',',header=TRUE)
head(bc)
class(bc$treat)
class(bc$time)
class(bc$y)
as.character(bc$treat)
as.character(bc$time)
as.numeric(bc$y)
as.character(bc$id)
as.numeric(bc$age)
library(Matrix)
library(lme4)
install.packages("lmerTest")
library(lmerTest)
model1 <- lmerTest::lmer(y ~treat + (1 | id), data = bc,REML = FALSE)
summary(model1)
model2 <- lmerTest::lmer(y ~ treat+time+treat * time + (1 | id), data = bc,REML = FALSE)
summary(model2)
model3 <- lmerTest::lmer(y ~ treat+time+treat * time + (1+time| id), data = bc,REML = FALSE)
summary(model3)
 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值