Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- // This is a two-point model. So First im predicting rain or no rain (0,1) using logistic regression then using lm model to predict the quantity of rain
- mod_LR=glm( Rainfall_LR ~ ., data = train_main[,-8] , family = binomial(link = logit),offset = Pressure)
- phat <- predict(mod_LR,test_main , type = "response")
- pred <- data.frame(test_main)
- logols.fit <- lm(log(Rainfall_Next3)~., train_main[train_main$Rainfall_Next3 > 0, -9])
- pred$logols <- phat * exp(predict(logols.fit, test_main))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement