data(proact)
## Testing effects of repeated A-B presentations
spec <- mptspec(
.BC = p*q*r,
.Bc = p*q*(1 - r),
.bC = p*(1 - q)*r,
.bc = (1 - p) + p*(1 - q)*(1 - r),
.replicates = 6
)
m1 <- mpt(spec, proact)
m2 <- mpt(update(spec, .restr = list(q2=q1, q3=q1, q5=q4, q6=q4)), proact)
m3 <- mpt(update(spec, .restr = list(r2=r1, r3=r1, r5=r4, r6=r4)), proact)
anova(m2, m1) # q increases with number of A-B presentations
anova(m3, m1) # r remains constant
Run the code above in your browser using DataLab