data(proact)
## Testing hypotheses about the parameters
mpt1 <- mpt(freq ~ list(
p1*q1*r1,
p1*q1*(1 - r1),
p1*(1 - q1)*r1,
(1 - p1) + p1*(1 - q1)*(1 - r1),
p2*q2*r2,
p2*q2*(1 - r2),
p2*(1 - q2)*r2,
(1 - p2) + p2*(1 - q2)*(1 - r2),
p3*q3*r3,
p3*q3*(1 - r3),
p3*(1 - q3)*r3,
(1 - p3) + p3*(1 - q3)*(1 - r3),
p4*q4*r4,
p4*q4*(1 - r4),
p4*(1 - q4)*r4,
(1 - p4) + p4*(1 - q4)*(1 - r4),
p5*q5*r5,
p5*q5*(1 - r5),
p5*(1 - q5)*r5,
(1 - p5) + p5*(1 - q5)*(1 - r5),
p6*q6*r6,
p6*q6*(1 - r6),
p6*(1 - q6)*r6,
(1 - p6) + p6*(1 - q6)*(1 - r6)
), proact, maxit=2000, stepsize=1.175)
mpt2 <- mpt(mpt1$formula, proact,
constr=list(q1 = c("q1","q2","q3"), q2 = c("q4","q5","q6")),
start=coef(mpt1))
mpt3 <- mpt(mpt1$formula, proact,
constr=list(r1 = c("r1","r2","r3"), r2 = c("r4","r5","r6")),
start=coef(mpt1), stepsize=1.9)
anova(mpt2, mpt1) # q increases with number of A-B presentations
anova(mpt3, mpt1) # r remains constant
Run the code above in your browser using DataLab