capture erase "c:\temp\temp.dta";
global i = 72;
   while $i <=88
     {;
        use $tradedata\imp$i , clear;
        gen year = yr+1900;
        rename name country;
        rename cusvalue valimp;
        gen  tempsic72 = real(substr(sic72,1,4)); drop sic72;
        rename tempsic sic72;
        gen cnum = real(ccode);
        keep cnum ccode country sic72 valimp duty year;
        drop if cnum == 999999 | cnum == 999991;
        gen oecddummy=0;
        replace oecddummy=1 if
              cnum ==    710360 | cnum ==    550400 | cnum ==    530560 | cnum ==    211240
            | cnum ==    532080 | cnum ==    552460 | cnum ==    532500 | cnum ==    532800
            | cnum ==    533000 | cnum ==    553520 | cnum ==    533720 | cnum ==    533800
            | cnum ==    413920 | cnum ==    535280 | cnum ==    715540 | cnum ==    555780
            | cnum ==    536200 | cnum ==    537240 | cnum ==    557520 | cnum ==    557560
            | cnum ==    447920 | cnum ==    538260 ;
        gen valimpLDC = valimp*(oecddummy==0);
        gen valimpmex = valimp*(country=="MEXICO");
        gen valimpcan = valimp*(country=="CANADA");
        egen valimpall = sum(valimp), by(sic72 year);
        keep if sic72 >=2000 & sic72< 4000;
       collapse (sum) valimpmex valimpcan valimpLDC (mean) valimpall, by(sic72 year);
            if $i==72
                  {;
                         save "c:\temp\temp";
                   };
            else if $i <= 88
                   {;
                        append using "c:\temp\temp";
                        save "c:\temp\temp", replace;
                    };        
           display $i;
           global i=$i+1;
         };      
                                       
    save c:\temp\imp7288a, replace;