replace primaryprovince = subinstr(primaryprovince, "Federally Administer", "Federally Administered", .)
replace primaryprovince = subinstr(primaryprovince, "Jowzjan", "Jawzjan", .)
replace primaryprovince = subinstr(primaryprovince, "KKhost", "Khost", .)
replace primaryprovince = subinstr(primaryprovince, "SGhazni", "Ghazni", .)
replace primaryprovince = subinstr(primaryprovince, "S Ghazni", "Ghazni", .)
replace primaryprovince = subinstr(primaryprovince, "Sar-e Pul", "Sari Pul", .)
replace primaryprovince = subinstr(primaryprovince, "a Herat", "Hirat", .)
replace primaryprovince = subinstr(primaryprovince, "aPaktika", "Paktika", .)
replace primaryprovince = subinstr(primaryprovince, "dPaktika", "Paktika", .)
replace primaryprovince = subinstr(primaryprovince, "dParwan", "Parwan", .)
replace primaryprovince = subinstr(primaryprovince, "eKhost", "Khost", .)
replace primaryprovince = subinstr(primaryprovince, "gSamangan", "Samangan", .)
replace primaryprovince = subinstr(primaryprovince, "hKapisa", "Kapisa", .)
replace primaryprovince = subinstr(primaryprovince, "hLogar", "Logar", .)
replace primaryprovince = subinstr(primaryprovince, "hWardak", "Wardak", .)
replace primaryprovince = subinstr(primaryprovince, "iHelmand", "Helmand", .)
replace primaryprovince = subinstr(primaryprovince, "lBaghlan", "Baghlan", .)
replace primaryprovince = subinstr(primaryprovince, "lSamangan", "Samangan", .)
replace primaryprovince = subinstr(primaryprovince, "nBalkh", "Balkh", .)
replace primaryprovince = subinstr(primaryprovince, "nHerat", "Herat", .)
replace primaryprovince = subinstr(primaryprovince, "uHerat", "Herat", .)
replace primaryprovince = subinstr(primaryprovince, "Day Kundi", "Daykundi", .)
replace primaryprovince = subinstr(primaryprovince, "Hirat", "Herat", .)
replace primaryprovince = subinstr(primaryprovince, ") Herat", "Herat", .)
replace primaryprovince = subinstr(primaryprovince, "JGhazni", "Ghazni", .)
replace primaryprovince = subinstr(primaryprovince, "Paktya", "Paktiya", .)

* create year-month  date
gen modate = ym(year1, month)
format modate %tm

* aggregate
contract primaryprovince modate, zero
rename _freq sigacts
label var sigacts "Significant Activities (SIGACTs) Totals"
label var primaryprovince "Province"
label var modate "Time"
label var modate "Period"

* remove missing
drop if primaryprovince==""
drop if modate==.


* Define a list of provinces
local primaryprovince "Badakhshan" "Badghis" "Baghlan" "Balkh" "Bamyan" "Daykundi" "Farah" "Faryab" "Federally Administered Tribal Areas" "Ghazni" "Ghor" "Helmand" "Herat" "Jawzjan" "Kabul" "Kandahar" "Kapisa" "Khost" "Khyber Pakhtunkhwa" "Kunar" "Kunduz" "Laghman" "Logar" "Nangarhar" "Nimroz" "Nuristan" "Paktika" "Paktiya" "Panjshir" "Parwan" "Samangan" "Sar-e Pol" "Takhar" "Uruzgan" "Wardak" "Zabul"



* Set default font to Times New Roman
set scheme s1mono

* Loop over each province
foreach primaryprovince of local primaryprovince {

    * Generate a line plot of sigacts vs modate for the current province
    twoway line sigacts modate if primaryprovince == "`primaryprovince'",
        title("Monthly SIGACTs in `primaryprovince', 2005-2014")
        xtitle("") ytitle("Total")
        legend(off)
        tlabel(2005m1(12)2014m12)
        xlabel(, angle(vertical))

        * Shade the july - January periods in each year
        || recast(area)
        line(, y(0) lpattern(shortdash))
        if primaryprovince == "`primaryprovince'" &
        month(modate) >= 7 &
        (month(modate) <= 12 | month(modate) == 1),
        color(gs12)

        * Add a July marker on the X axis for each year
        || scatteri
        7, ypos(0) mcolor(black) msize(large) mlab("July")
        if primaryprovince == "`primaryprovince'" &
        month(modate) == 7,

        * Add a demarcation line on the Prvoince "Kabul" for August 2003 labeled "NATO deployment"
        || line
        21.5, yscale(alt) lcolor(blue) lpattern(dash)
        text(, position(11 21.5) color(blue) size(large))
        if primaryprovince == "`primaryprovince'" &
        month(modate) == 8 &
        year(modate) == 2003,
        yline(, lpatten(shortdash))

    * Save the plot as an image
    graph export "sigacts_`primaryprovince'.png"
}