- Code: Select all
w<-c(test.ddl$Psi$fix[test.ddl$Psi$stratum=="B"& test.ddl$Psi$tostratum=="b"],
test.ddl$Psi$fix[test.ddl$Psi$stratum=="C"& test.ddl$Psi$tostratum=="c"])
test.ddl$Psi$fix[test.ddl$Psi$stratum=="A"& test.ddl$Psi$tostratum=="a"] <- mean(w, na.rm = FALSE)
I also tried fixing PsiAa to PsiBb using the code below which also didn't give the same estimate for both transitions.
- Code: Select all
test.ddl$Psi$fix[test.ddl$Psi$stratum=="A"& test.ddl$Psi$tostratum=="a"] <-test.ddl$Psi$fix[test.ddl$Psi$stratum=="B"& test.ddl$Psi$tostratum=="b"]
Thank you for your help!
-RB