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)
重复测量资料多水平模型
最新推荐文章于 2024-09-13 10:46:15 发布