Hi,
I wish to perform score tests comparing the fit of a null secr model model (hazard rate detection function) to models with session effects on each of g0 and sigma, and want to be sure I'm using the correct betaindex specification. I think the order of parameters is 1 = g0, 2 = sigma, 3 = z for hazard rate. Beta index numbers should then be as follows, for a data set with 9 sessions?
# SCORE TEST SESSION EFFECTS
st.g0SESS=score.test(g0.sig.fit, g0~session, betaindex = c(1, 10, 11))
st.sigSESS=score.test(g0.sig.fit, sigma~session, betaindex = (1, 2, 11))
#
Thanks,
Eric Howe