I have a plot of residual values of a linear model in function of the fitted values where the heteroscedasticity is very clear. However I'm not sure how I should proceed now because as far as I understand this heteroscedasticity makes my linear model invalid right?
So I have been doing some reading about this subjet and I found two suggestions in other stackoverflow threads.
1) Use robut linear fitting using the rlm() function of the MASS package because it's apparently robust to heteroscedasticity.
2) As the standard errors of my coefficients are wrong because of the heteroscedasticity, I can just adjust the standard errors to be robust to the heteroskedasticity? Using the method posted here
Which would be the best method to use to deal with my problem? If I use solution 2 is my predicting capability of my model completely useless?
My residuals in function of the fitted values looks like this http://i.gyazo.com/9407a829a168492b31dfa3d1dd33a21d.png and the Breusch-Pagan test confirmed that the variance is not constant.