[R-sig-ME] longitudinal analysis when one group switched from control to treatment

Simon Harmel @|m@h@rme| @end|ng |rom gm@||@com
Sun May 17 01:09:17 CEST 2020


Hello All,

I have a 3-year longitudinal dataset (*see link below the table*). Up to
year 2 (coded "1"), 8 schools (4 in Treatment, 4 in Control) cooperated
with the study. But in year 3 (coded "2"), one of the Treatment schools
(named "good") dropped out.

Also in year 3 (coded "2"), we were made to move one of the *Control *schools
(named "*orange*") to the *Treatment *group. The full design of the study
is shown in the Table below.

I want to regress "year" and "group" on "y" (a continuous response) in lme4
package in R. But is there a way to capture the switch of one of the
control schools to the treatment group?

Thank you very much, Simon

·       *Switched from control to treatment*

·       *Out as of year coded 2*

*SCHOOL NAMES*

*Year*

*Codes*

*Control*

*Treatment*

0

har

john

bus

orange

caro

good

bla

carm

1

har

john

bus

*orange*

caro

good

bla

carm

2

har

john

bus

X

caro

*orange*

bla

carm

*library(lme4)*
*dat <- read.csv('https://raw.githubusercontent.com/hkil/m/master/z.csv
<https://raw.githubusercontent.com/hkil/m/master/z.csv>')*

*m1 <- lmer(y~ year*group + (1|stid), data = dat)      #### 'stid' =
student id                m2 <- lmer(y~ year*group + (1|scid/stid), data =
dat) #### 'scid' = school id*

	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list