diff --git a/NEED/analyse-NEED-EULF-2014-models-v2.0.do b/NEED/analyse-NEED-EULF-2014-models-v2.0.do
index 8b9ba19eb4786258bc52e26d090ff57210ebbd9e..4ba6603dbee47f0da82c19c63c2dbfdcd0e25145 100644
--- a/NEED/analyse-NEED-EULF-2014-models-v2.0.do
+++ b/NEED/analyse-NEED-EULF-2014-models-v2.0.do
@@ -182,13 +182,14 @@ foreach v of local vars {
 	
 	est store `v'
 	* test a variable
-	predict r, resid
+	predict `v'_r, resid
 	foreach testv of local generic_hvarsg {
 		di "Testing residuals against `v'"
 		* based on http://www.ats.ucla.edu/stat/stata/webbooks/reg/chapter2/statareg2.htm
 		* can't use factor variables in acprplot...
-		graph box r, over(`testv') name(rtest_`testv')
+		graph hbox `v'_r, over(`testv') name(rtest_`testv')
 		graph export "`rpath'/graphs/rplot_`v'_`testv'-`version'.png", replace
+		* possibly floor area - heteroskedasticity ?
 	}