-
Notifications
You must be signed in to change notification settings - Fork 0
/
time_to_trt_graph.do
56 lines (39 loc) · 1.93 KB
/
time_to_trt_graph.do
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
global loc "/Users/christianbaehr/Desktop"
import delimited "$loc/time_to_trt1_trimmed.csv", clear
regress coef time if time < 0
predict a1
twoway (line coef time, lcolor(blue)) ///
(line ci1 time, lcolor(blue) lpattern(dash)) ///
(line ci2 time, lcolor(blue) lpattern(dash)) ///
(lfit coef time if time < 0, lcolor(red)) ///
(line a1 time, lcolor(red) lpattern(dot)) ///
(lfit coef time if time > 0, lcolor(red)), ///
xline(0, lcolor(black)) yline(0, lcolor(black)) xlabel(-4(1)10) legend(off) graphregion(fcolor(white)) ///
title("Time to Treatment (0-1km)") xtitle("Years to Treatment") ytitle("Treatment Effects on NDVI")
graph export "$area/time_to_trt1.png", replace
***
import delimited "$area/time_to_trt2.csv", clear
regress coef time if time <= 0
predict a1
twoway (line coef time, lcolor(blue)) ///
(line ci1 time, lcolor(blue) lpattern(dash)) ///
(line ci2 time, lcolor(blue) lpattern(dash)) ///
(lfit coef time if time <= 0, lcolor(red)) ///
(line a1 time, lcolor(red) lpattern(dot)) ///
(lfit coef time if time >= 0, lcolor(red)), ///
xline(0, lcolor(black)) yline(0, lcolor(black)) xlabel(-4(1)10) legend(off) graphregion(fcolor(white)) ///
title("Time to Treatment (1-2km)") xtitle("Years to Treatment") ytitle("Treatment Effects on NDVI")
graph export "$area/time_to_trt2.png", replace
***
import delimited "$area/time_to_trt3.csv", clear
regress coef time if time <= 0
predict a1
twoway (line coef time, lcolor(blue)) ///
(line ci1 time, lcolor(blue) lpattern(dash)) ///
(line ci2 time, lcolor(blue) lpattern(dash)) ///
(lfit coef time if time <= 0, lcolor(red)) ///
(line a1 time, lcolor(red) lpattern(dot)) ///
(lfit coef time if time >= 0, lcolor(red)), ///
xline(0, lcolor(black)) yline(0, lcolor(black)) xlabel(-4(1)10) legend(off) graphregion(fcolor(white)) ///
title("Time to Treatment (2-3km)") xtitle("Years to Treatment") ytitle("Treatment Effects on NDVI")
graph export "$area/time_to_trt3.png", replace