- Code: Select all
multigroup.process=process.data(multigroup,model="CJS",time.intervals=c(1.144,10.2,.922,10.467,.822),begin.time=0,
groups=c("sex","trt"))
multigroup.ddl=make.design.data(multigroup.process)
multigroup.ddl$Phi$season=0
multigroup.ddl$Phi$season[multigroup.ddl$Phi$time==1.144]=1
multigroup.ddl$Phi$season[multigroup.ddl$Phi$time==12.266]=1
multigroup.ddl$p$season=0
multigroup.ddl$p$season[multigroup.ddl$p$time== 11.344]=1
multigroup.ddl$p$season[multigroup.ddl$p$time==22.733 ]=1
Thanks,
Amanda