Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- breaks <- V_dep[V_dep >= ... & V_dep <= ...]
- mse <- numeric(length(breaks))
- for(i in 1:length(breaks)){
- piecewise <- lm(V_indep ~ V_dep*(V_dep < breaks[i]) + V_dep*(V_dep >=breaks[i]))
- mse[i] <- summary(piecewise)[6]
- }
- mse <- as.numeric(mse)
- breaks[which(mse==min(mse))]
Add Comment
Please, Sign In to add comment