interaction with dummy variables in lm() in R -
ref: http://www.r-bloggers.com/r-for-ecologists-putting-together-a-piecewise-regression/
in paper, confused argument:
y ~ x*(x < breaks[i]) + x*(x>=breaks[i]) in lm().
i know * in lm means interactions , main effects mean predictors x x (x < breaks[i]) (x < breaks[i]) , interactions?
this method of doing "segmented" regression. creating 2 different models, 1 section x < breaks[i] , opposite true. webpage seems pretty nice job of illustrating it's unclear missing. model formula might more clear if written as:
y ~ x*i(x < breaks[i]) + x*i(x>=breaks[i]) it means there 2 predictors: first 1 beingness x , sec 1 beingness logical vector 1 in part less breaks[i] , 0 in other region. in fact not need 2 terms in model if used:
y ~ x*i(x < breaks[i]) i thought predictions same, different, perhaps because 2 term model implicitly allowed independent intercepts.
r
No comments:
Post a Comment