# using hardman.m02, only Trial as predictor
# in fixed part and in random part at Listener, Talker, Item levels
# 6 talkers
# no lines, only axes
curve( ranef(hardman.m03)$talker[1,1]+ranef(hardman.m03)$talker[1,2]*(x-30.5),
from=1, to=60, ylim=c(-2,+2), xlab="Trial", ylab="log odds", type="n" )
for (j in 1:6) {
curve( ranef(hardman.m03)$talker[j,1]+ranef(hardman.m03)$talker[j,2]*(x-30.5),
add=TRUE, lwd=2 )
}
for (j in 1:6) {
text( x=30.5, y=ranef(hardman.m03)$talker[j,1]+0.1,
row.names(ranef(hardman.m03)$talker)[j] )
}
title( main="Trial effect varies between talkers (Hardman data)" )