Stefan Th. Gries
unread,Feb 6, 2017, 8:45:03 PM2/6/17Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to StatForLing with R
Some comments (one with some self-promotion, sorry, but the fit is too
close to ignore):
- I don't think you need rescaling: the significance tests and the
predicted values shouldn't be affected by that much;
- I do think that, if you enter any TIME predictor - numeric, ordinal,
categorical - it needs to also interact with probably all other
predictors;
- if you enter TIME as a numeric predictor, using a straight-line
regression is probably way too simplistic because doing that means
you're hypothesizing a linear unchanged increase over time (which I
cannot believe you would want to subscribe to) - thus, I'd recommend
using a structure for TIME that allows for curvature: a simple way is
using a polynomial (2nd or 3rd degree), a complex one is fitting a
GAMM;
- you may consider first grouping together time points into time
periods and now I have to refer to Gries & Hilpert (2010k), where we
have very erratic temporal data and use a method I called VNC to find
structure in the temporal data so that, then, TIME could be entered
into a mixed-effects model as a categorical predictor.