/* Rebuild Data file for PEPS Project */ /* - using PWT 6.1, Polity 4 (2004), Freedom House & Vanhanen */ /* */ /********************************************************************/ /* 1. Merge Freedom House with Polity 4 */ /********************************************************************/ /* Begin with STATA Freedom House file previously converted during /// filebuild for ISA 2003 project with Bill freehouse.dat->democ.dta*/ /*adding 2001 & 2002 & 2003 data from FH website summer 2004 */ insheet using "C:\StataData\FHaddtoStata.csv", comma clear // read Excel file of 2001-2003 additions destring pol civ, replace force gen free = 14 - pol -civ summarize *tab ccode tab year sort ccode year merge ccode year using freehouse.dta // merge additions with previous FH data 1972-2000 tab _merge sort ccode year *use freehouse, clear replace name = "Colombia" if name=="Columbia", nopromote replace name = "Romania" if name=="Rumania", nopromote replace name = "Cameroon" if name=="Cameroun", nopromote /* Check for name changes */ gen ccodelag = ccode[_n-1] replace name = name[_n-1] if name=="" & ccode==ccodelag /// replace blank names because // 2001-2003 additions do not have a country name gen namelag = name[_n-1] gen namechange = match(name,namelag) list if ccode==ccodelag & namechange==0 rename pol FHpolrights rename civ FHcivliberties rename free FHfree label variable FHfree "Sum with perfect 12" /* keep ccode year name free */ /* Stupid ccode assignments require fixing... */ drop if ccode==565 & year==1989 /* drop Namibia to correspond to Polity */ * drop if ccode==265 & year==1990 /* drop E Germy to correspond to Polity */ replace ccode=678 if ccode==679 /* Yemen, reconciliation */ replace ccode=232 if ccode==233 /* Andorra, reconciliation */ replace ccode=316 if ccode==315 & year>1992 /* Czech Repulic, reconciliation */ drop if ccode==403 & year < 1975 /*Sao Tome & Principe suspect data */ summarize drop ccodelag namelag namechange _merge generate nFH = 1 sort ccode year summarize save tmp/freehouse2004, replace /********************************************************************/ /* Clean up Polity 4 versions 2000 & 2002 */ *use p4v2000, clear /* Temporarily access 2000 version of Polity from Bill */ /* just long enough to save the records from Luxemourg & Iceland */ /* that are deleted from the 2002 version */ use p4v2000, clear keep if ccode==212 | ccode==395 summarize save tmp/addLuxIce, replace /* Now read Polity 4, version 2002, acquired from Polity site */ insheet using "C:\Data\p4v2002.csv", clear summarize append using "C:\StataData\tmp\addLuxIce.dta", nolabel nonotes // Add in Lux & Ice summarize rename democ Pdemocracy rename autoc Pautocracy rename polity Polity1 rename polity2 Polity2 rename durable Pdurable rename change Pchange rename regtrans Pregtrans keep ccode country year Pdemocracy Pautocracy Polity1 Polity2 Pdurable Pchange Pregtrans sort ccode year save tmp/p4v2002temp, replace generate nP = 1 /* Stupid Polity ccode assignments require fixing... */ drop if ccode==255 & year==1990 /* duplicate case */ drop if ccode==347 & year==1991 /* duplicate case */ drop if ccode==365 & year==1922 /* duplicate case */ drop if ccode==530 & year==1993 /* duplicate Ethiopia */ replace ccode=260 if ccode==255 /* Germany */ replace ccode=345 if ccode==347 /* Yugoslavia */ * replace ccode=315 if ccode==316 /* Czechoslovakia */ replace ccode=365 if ccode==364 /* Russia */ * drop if ccode==531 /* Eritrea */ replace ccode=530 if ccode==529 /* Ethiopia */ drop if ccode==679 & year==1990 /* Yemen duplicate */ * drop if ccode==678 & year > 1989 /* Yeman */ replace ccode=678 if ccode==679 /* Yemen */ replace ccode=770 if ccode==769 /* Pakistan */ drop if ccode==818 & year==1976 /* Vietnam */ replace ccode=816 if ccode==818 /* Vietnam */ * Correct Polity2 scores which have a generic missing data code replace Polity2=Polity1 if Polity2>10 * Create a Polity3 that does not recode -77 to 0 gen Polity3 = Polity2 replace Polity3 = Polity1 if Polity1 == -77 * Interpolate Pdemocracy so that PEPS3 can be computed for cases * in which Polity2&3 were interpolated by the Polity coders replace Pdemocracy = 3 if ccode==40 & year==1952 replace Pdemocracy = 2 if ccode==40 & year==1953 replace Pdemocracy = 1 if ccode==40 & year==1954 replace Pdemocracy = 0 if ccode==41 & year==1986 replace Pdemocracy = 0 if ccode==41 & year==1987 replace Pdemocracy = 4 if ccode==41 & year==1999 replace Pdemocracy = 3 if ccode==42 & year==1961 replace Pdemocracy = 6 if ccode==42 & year==1962 replace Pdemocracy = 2 if ccode==90 & year==1985 replace Pdemocracy = 2 if ccode==91 & year==1980 replace Pdemocracy = 5 if ccode==91 & year==1981 replace Pdemocracy = 1 if ccode==92 & year==1979 replace Pdemocracy = 2 if ccode==92 & year==1980 replace Pdemocracy = 3 if ccode==92 & year==1981 replace Pdemocracy = 4 if ccode==92 & year==1982 replace Pdemocracy = 5 if ccode==92 & year==1983 replace Pdemocracy = 2 if ccode==135 & year==1978 replace Pdemocracy = 5 if ccode==135 & year==1979 replace Pdemocracy = 2 if ccode==140 & year==1964 replace Pdemocracy = 1 if ccode==145 & year==1952 replace Pdemocracy = 1 if ccode==145 & year==1953 replace Pdemocracy = 1 if ccode==145 & year==1954 replace Pdemocracy = 1 if ccode==145 & year==1955 replace Pdemocracy = 1 if ccode==160 & year==1955 replace Pdemocracy = 2 if ccode==160 & year==1956 replace Pdemocracy = 5 if ccode==170 & year==1971 replace Pdemocracy = 2 if ccode==170 & year==1972 replace Pdemocracy = 2 if ccode==230 & year==1975 replace Pdemocracy = 4 if ccode==230 & year==1976 replace Pdemocracy = 6 if ccode==230 & year==1977 replace Pdemocracy = 3 if ccode==235 & year==1974 replace Pdemocracy = 6 if ccode==235 & year==1975 replace Pdemocracy = 0 if ccode==310 & year==1956 replace Pdemocracy = 6 if ccode==310 & year==1989 replace Pdemocracy = 0 if ccode==315 & year==1968 replace Pdemocracy = 3 if ccode==344 & year==1999 replace Pdemocracy = 0 if ccode==345 & year==1951 replace Pdemocracy = 0 if ccode==345 & year==1952 replace Pdemocracy = 4 if ccode==350 & year==1974 replace Pdemocracy = 2 if ccode==360 & year==1989 replace Pdemocracy = 3 if ccode==432 & year==1991 replace Pdemocracy = 1 if ccode==433 & year==1962 replace Pdemocracy = 2 if ccode==434 & year==1963 replace Pdemocracy = 1 if ccode==434 & year==1964 replace Pdemocracy = 3 if ccode==434 & year==1990 replace Pdemocracy = 4 if ccode==436 & year==1991 replace Pdemocracy = 2 if ccode==437 & year==1999 replace Pdemocracy = 3 if ccode==439 & year==1977 replace Pdemocracy = 2 if ccode==452 & year==1969 replace Pdemocracy = 3 if ccode==452 & year==1978 replace Pdemocracy = 0 if ccode==452 & year==1991 replace Pdemocracy = 0 if ccode==461 & year==1991 replace Pdemocracy = 1 if ccode==461 & year==1992 replace Pdemocracy = 4 if ccode==475 & year==1978 replace Pdemocracy = 2 if ccode==475 & year==1998 replace Pdemocracy = 2 if ccode==461 & year==1998 replace Pdemocracy = 0 if ccode==481 & year==1990 replace Pdemocracy = 3 if ccode==484 & year==1991 replace Pdemocracy = 3 if ccode==500 & year==1966 replace Pdemocracy = 1 if ccode==516 & year==1965 replace Pdemocracy = 0 if ccode==517 & year==1993 replace Pdemocracy = 7 if ccode==552 & year==1979 replace Pdemocracy = 7 if ccode==560 & year==1992 replace Pdemocracy = 8 if ccode==560 & year==1993 replace Pdemocracy = 5 if ccode==580 & year==1991 replace Pdemocracy = 4 if ccode==625 & year==1964 replace Pdemocracy = 5 if ccode==625 & year==1969 replace Pdemocracy = 3 if ccode==625 & year==1970 replace Pdemocracy = 4 if ccode==625 & year==1985 replace Pdemocracy = 2 if ccode==630 & year==1953 replace Pdemocracy = 1 if ccode==630 & year==1954 replace Pdemocracy = 0 if ccode==630 & year==1979 replace Pdemocracy = 0 if ccode==630 & year==1980 replace Pdemocracy = 0 if ccode==630 & year==1981 replace Pdemocracy = 7 if ccode==640 & year==1960 replace Pdemocracy = 6 if ccode==165 & year==1971 replace Pdemocracy = 3 if ccode==165 & year==1972 replace Pdemocracy = 0 if ccode==678 & year==1991 replace Pdemocracy = 0 if ccode==678 & year==1992 replace Pdemocracy = 0 if ccode==690 & year==1990 replace Pdemocracy = 4 if ccode==732 & year==1987 replace Pdemocracy = 1 if ccode==790 & year==1957 replace Pdemocracy = 3 if ccode==790 & year==1958 replace Pdemocracy = 2 if ccode==800 & year==1968 replace Pdemocracy = 2 if ccode==800 & year==1973 replace Pdemocracy = 2 if ccode==800 & year==1977 *replace Pdemocracy = 4 if ccode==840 & year==1968 replace Pdemocracy = 4 if ccode==840 & year==1986 replace Pdemocracy = 6 if ccode==950 & year==2000 replace Pdemocracy = 0 if ccode==811 & year==1970 replace Pdemocracy = 0 if ccode==811 & year==1971 label variable Polity1 "Raw Polity score MD included" label variable Polity3 "Preferred Polity score" /* drop if year < 1950 */ /* Better keep years prior to 1950 */ /* Occupations (-66) did not change level of democ in these states: */ * replace Polity2=-7 if ccode==652 & year==1958 /* Syria */ * replace Polity2=-7 if ccode==652 & year==1959 /* Syria */ * replace Polity2=-7 if ccode==652 & year==1960 /* Syria */ replace Polity2=-7 if ccode==310 & year==1956 /* Hungary */ replace Polity2=-7 if ccode==315 & year==1968 /* Czechoslovakia */ * replace Polity2=-10 if ccode==690 & year==1990 /* Kuwait */ * replace Polity3=-7 if ccode==652 & year==1958 /* Syria */ * replace Polity3=-7 if ccode==652 & year==1959 /* Syria */ * replace Polity3=-7 if ccode==652 & year==1960 /* Syria */ replace Polity3=-7 if ccode==310 & year==1956 /* Hungary */ replace Polity3=-7 if ccode==315 & year==1968 /* Czechoslovakia */ * replace Polity3=-10 if ccode==690 & year==1990 /* Kuwait */ *mvdecode Pdemocracy Pautocracy Polity1, mv(-66 -77 -88) mvdecode Pdemocracy Pautocracy Polity1 Polity2 Polity3, mv(-66=.a \ -77=.b \ -88=.c) sort ccode year /********************************************************************/ /* Now merge FH (democ) with cleaned Polity */ merge ccode year using tmp/freehouse2004 /* summarize */ drop nP nFH replace country=name if country=="" rename _merge mergePFH drop name sort ccode year save tmp/PolFH, replace /********************************************************************/ /* 3. Merge Vanhanen */ /* (vanhanen STATA file created from xls file acquired on Web) */ *Step 1 : Read revised Vanhanen data (.csv) into Stata & save note: .csv file created by Patrick Schmid from original Vanhanen Access file, bypassing weights ** insheet using "C:\STATADATA\vanhanen_ALL.csv",clear insheet using vanhanendataexpanded.csv, clear /********************************************************************/ sort ccode year election_type *Correct Vanhanen population error replace population=population*2 if ccode==481& year<1979 //Gabon *Calculate toV in place of vote_share, which contains only two significant digits gen toV = vote_share replace toV = votes/population if votes > 0 & population > 0 replace toV = vote_share if toV >1000 pwcorr summarize *Remove duplicate records when multiple elections are held in same year duplicates tag ccode year election_type, generate(dup) //find duplicates *summarize *duplicates list ccode year election_type, separator(5) //list duplicates *list if dup>0 drop if dup>0 & population==0 & year==year[_n-1] //drop the second of two zero duplicates drop dup duplicates tag ccode year election_type, generate(dup) //find remaining duplicates drop if dup>0 & population==0 & year==year[_n+1] //drop the first of two zero duplicates *list if dup>0 drop dup duplicates tag ccode year election_type, generate(dup) //find remaining duplicates *duplicates list ccode year election_type, separator(5) //list remaining duplicates *list if dup>0 drop if dup>0 & year==year[_n-1] & toV<= toV[_n-1] //drop the second of duplicates if smaller drop dup duplicates tag ccode year election_type, generate(dup) *list if dup>0 drop if dup>0 & year==year[_n+1] & toV< toV[_n+1] //drop the first of duplicates if smaller drop dup duplicates tag ccode year election_type, generate(dup) duplicates list ccode year election_type, separator(5) drop votes vote_share population dup summarize *Reshape data set by pulling A & B elections into the same record, * then choosing the larger turnout of the two reshape wide toV, i( ccode year ) j( election_type ) string summarize mvencode toVA toVB, mv(.=-1) gen toV = toVA replace toV = toVB if toVB>toVA list if toV<=0 & (toVA>0 | toVB>0) replace toV=toV*100 /66.67 replace toV=1 if toV>1 & toV<1000 summarize pwcorr drop toVA toVB *Fix country codes replace ccode=305 if ccode==300 replace ccode=315 if ccode==316 replace ccode=678 if ccode==679 replace ccode=260 if ccode==255 replace ccode=342 if ccode==345 & year<=1920 // Recode from Yugoslavia to Serbia replace ccode=99 if ccode==100 & year<=1831 // Recode from Colombia to GranColombia replace ccode=316 if ccode==315 & year>=1993 // Recode from Czechoslovakia to Czech Rep drop if ccode==260 & year==1945 & toV==0 duplicates tag ccode year, generate(dup) duplicates list ccode year , separator(5) list if dup>0 drop dup sort ccode year tsset ccode year, yearly **drop if year>2002 drop if year>2004 rename toV toVnew sort ccode year **save tmp/vanhanen_june , replace save tmp/vanhanen_july , replace // // Now merge Vanhanen with Polity & Freedom House merge ccode year using tmp/PolFh, _merge(_mergeV) *merge ccode year using "C:\STATADATA\vanhanen.dta", _merge(_mergeV) tabulate _mergeV summarize *tab ccode if _mergeV==1 *drop if _mergeV==1 & (part==0 | part>1000) *replace country=vabb if country=="" *rename comp Vcompete *rename id Vindex *rename part Vparticipate *replace Vparticipate = 66.67 if Vparticipate > 66.67 & Vparticipate < 1000 rename toVnew toV generate PEPS1v = (Pdemocracy * toV) - Pautocracy replace PEPS1v = (0 - Pautocracy) if Pdemocracy == 0 generate PEPS2v = (((toV/.05)-10) + Polity3) / 2 summarize drop if ccode==260 & year==1945 & country=="Germany West" duplicates tag ccode year, generate(dup) duplicates list ccode year , separator(5) list if dup>0 sort ccode year tsset ccode year, yearly compress save tmp/PolFhVan2, replace //4. Fix IDEA // * insheet using "C:\My Download Files\IDEA.zip\ExpandedIDEA_Code-Type-Year.csv",clear * insheet using "C:\STATADATA\tblfinalIDEA.315.csv",clear **insheet using IDEA-van48.csv, clear insheet using qryidea_un.csv, clear * Fix some codes */ drop if ccode==680 & year>1990 do later if necessary */ *drop if ccode ==678 & year > 1992 *replace ccode = 680 if ccode==678 *replace ccode = 680 if ccode==679 *replace ccode = 484 if ccode==490 & year > 1991 * Set year range NO KEEP THEM ALL */ drop if year>2000 */ */ drop if year<1950 */ * Drop calculated variables so they can later be re-calculated **drop touta toutb tout vapp van_vap vanhanen idea_tout van_tout drop touta toutb tout unvapp drop if year==2005 rename a_totalvote Avote rename b_totalvote Bvote rename totalvote Ivote rename vap IVAP **rename van_totalvote Vvote rename unvap UNVAP rename unpopulation UNpop rename a_vap AVAP rename b_vap BVAP **format %20s country format %9.0fc UNVAP AVAP Avote mvencode AVAP Avote BVAP Bvote, mv(. = -99) *Various corrections involving obvious IDEA misstatements of VAP replace AVAP = AVAP*2 if ccode==223 & year<1986 //count women in Liectenstein replace AVAP=UNVAP if ccode==225 & year<1975 //count women in Switzerland replace BVAP=UNVAP if ccode==225 & year<1975 replace AVAP=UNVAP if ccode==350 & year<1956 //count women in Greece replace BVAP=UNVAP if ccode==350 & year<1956 replace AVAP=UNVAP if ccode==690 & year>1991 //count women in Kuwait replace BVAP=UNVAP if ccode==690 & year>1991 replace AVAP=UNVAP if ccode==404 & year>1998 //Guinea-Bissau replace BVAP=UNVAP if ccode==404 & year>1998 *drop if ccode==317 & year==1993 *To diagnose source of differences between Vanhanen & IDEA data *gen pctVAP = UNVAP / UNpop *gen VAPdif = abs(UNVAP-BVAP)/BVAP *gen votediff = abs(Vvote-Avote)/Vvote if Vvote>0 & Avote>0 *format %5.2f VAPdif pctVAP *mvdecode Avote AVAP Bvote BVAP Ivote IVAP, mv(-1=.a \-2=.b \ -3=.c \ -4=.d \-5= .e \ -9=.f) gen toutA = Avote / AVAP replace toutA = -9 if AVAP<=0 replace toutA = Avote if Avote<=0 gen toutB = Bvote / BVAP replace toutB = -9 if BVAP<=0 replace toutB = Bvote if Bvote<=0 *gen toutV = Vvote / UNVAP *replace toutV = -9 if UNVAP<=0 *replace toutV = Vvote if Vvote<=0 * Create xtout, the larger of the two touts (Presidential or Parliamentary) * in a given year, as an alternative to tout. compare toutA toutB * 1. It should be coded missing if there is no IDEA data gen xtout = -9 replace xtout = toutA if toutA<0 & toutB<0 * 2. If either is missing, use the other replace xtout = toutA if toutA>=0 & toutB<0 replace xtout = toutB if toutB>=0 & toutA<0 * 3. if both are present, use the larger replace xtout = toutA if toutA>=0 & toutB >= 0 & toutA>=toutB replace xtout = toutB if toutA>=0 & toutB >= 0 & toutA< toutB replace xtout=1 if xtout>1 & xtout<1000 *Replace IDEA missing records with appropriate missing data codes replace xtout= -1 if ccode==40 & year<1961 replace xtout = -2 if ccode==40 & year>=1961 replace xtout = -2 if ccode==265 & year<1991 replace xtout = -1 if ccode==345 & year>1991 replace xtout = -2 if ccode==678 & year<1990 replace xtout = -2 if ccode==630 & year<1979 replace xtout = -2 if ccode==411 replace xtout = -9 if ccode==450 replace xtout = -2 if ccode==517 replace xtout = -9 if ccode==560 & year<1960 replace xtout = 0 if ccode==572 replace xtout = -1 if ccode==620 & year<1969 replace xtout = -2 if ccode==620 & year>1968 replace xtout = -1 if ccode==645 & year<1959 replace xtout = -2 if ccode==645 & year>1958 replace xtout = 0 if ccode==670 replace xtout = 0 if ccode==692 replace xtout = 0 if ccode==694 replace xtout = 0 if ccode==696 replace xtout = 0 if ccode==698 replace xtout = 0 if ccode==700 & year<1965 replace xtout = -2 if ccode==700 & year>1964 replace xtout = -2 if ccode==701 replace xtout = -2 if ccode==710 replace xtout = -2 if ccode==731 replace xtout = -2 if ccode==760 replace xtout = -2 if ccode==812 replace xtout = -2 if ccode==816 replace xtout = -2 if ccode==817 **list ccode iso country year AVAP BVAP xtout toutV if xtout==-99 & ccode!=ccode[_n-1] list ccode year AVAP BVAP xtout if xtout==-99 & ccode!=ccode[_n-1] tab ccode if xtout== -99 **mvdecode Avote AVAP toutA Bvote BVAP IVAP Ivote toutB xtout toutV, mv(-1=.a \-2=.b \ -3=.c \ -4=.d \-5= .e \ -9=.f \-99=.f) mvdecode Avote AVAP toutA Bvote BVAP IVAP Ivote toutB xtout , mv(-1=.a \-2=.b \ -3=.c \ -4=.d \-5= .e \ -9=.f \-99=.f) rename xtout toI label variable toutA "IDEA Votes/VAP, Parliamentary" label variable toutB "IDEA Votes/VAP, Presidential" label variable Avote "IDEA Votes, Parliamentary" label variable Bvote "IDEA Votes, Presidential" label variable AVAP "IDEA Voting age population, Parliamentary" label variable BVAP "IDEA Voting age population, Presidential" sort ccode year duplicates tag ccode year, generate(dup) duplicates list ccode year , separator(5) list if dup>0 // Now merge IDEA with Vanhanen, Polity & Freedom House merge ccode year using tmp/PolFhVan2, _merge(_mergeI) generate PEPS1i = (Pdemocracy * toI) - Pautocracy replace PEPS1i = (0 - Pautocracy) if Pdemocracy == 0 generate PEPS2i = (((toI/.05)-10) + Polity3) / 2 *gen toutV = Vvote / UNVAP *replace toutV = -9 if UNVAP<=0 *replace toutV = Vvote if Vvote<=0 summarize sort ccode year tsset ccode year, yearly compress save tmp/AllDem2, replace // /********************************************************************/ /* 2. Merge Summers & Heston, "Penn World Tables, 6.1" (1950-2000) */ /* (pwt61 STATA file created from Bill's zipfile) */ /* clean up file first */ /********************************************************************/ use pwt61, clear replace ccode=316 if ccode==315 & year>= 1993 sort ccode year /*Now merge */ merge ccode year using tmp/AllDem2 * mvencode ccode, mv(0) drop if ccode < 2 *drop if country=="Germany" & year==1945 drop xrat cc cg ci kc kg ki csave rgdpl cgdp rgdptt y p pc pg pi cgnp rgdpeqa rgdpwok /* mvencode ccode, mv(1) override */ tab _merge drop if _merge==1 & rgdpch==. tab _merge replace country=isocode if country=="" drop _merge sort ccode year tsset ccode year, yearly compress label data "PEPS1_1800-2003 built from Polity4 (2004), IDEA web, VanHanen, PWT6.1 " save tmp/PEPS1_1800-2003v2, replace *Create lags of Polity3, toutI, PEPS & PEPS2 gen P3L1 = L1.Polity3 gen P3L2 = L2.Polity3 gen P3L3 = L3.Polity3 gen P3L5 = L5.Polity3 gen P3L10 = L10.Polity3 gen P3L20 = L20.Polity3 gen P3FD1 = Polity3 - L1.Polity3 gen tVL1 = L1.toV gen tVL2 = L2.toV gen tVL3 = L3.toV gen tVL4 = L4.toV gen tVL5 = L5.toV gen tVL10 = L10.toV gen tVL20 = L20.toV gen tIL1 = L1.toI gen tIL2 = L2.toI gen tIL3 = L3.toI gen tIL4 = L4.toI gen tIL5 = L5.toI gen tIL10 = L10.toI gen tIL20 = L20.toI gen PEPS1vL1 = L1.PEPS1v gen PEPS1vL2 = L2.PEPS1v gen PEPS1vL3 = L3.PEPS1v gen PEPS1vL5 = L5.PEPS1v gen PEPS1iL1 = L1.PEPS1i gen PEPS1iL2 = L2.PEPS1i gen PEPS1iL3 = L3.PEPS1i gen PEPS1iL5 = L5.PEPS1i gen FHL1 = L1.FHfree gen PEPS2vL1 = L1.PEPS2v gen PEPS2vL5 = L5.PEPS2v gen PEPS2vL10 = L10.PEPS2v gen PEPS2vL20 = L20.PEPS2v gen PEPS2iL1 = L1.PEPS2i gen PEPS2iL5 = L5.PEPS2i gen PEPS2iL10 = L10.PEPS2i gen PEPS2iL20 = L20.PEPS2i exit *Random effects GLS regression of Polity on lagged toutI xtreg Polity3 tIL1 if P3L1 < 1000 & tIL1 < 1000, re xtreg Polity3 tIL2 if P3L2 < 1000 & tIL2 < 1000, re xtreg Polity3 tIL3 if P3L3 < 1000 & tIL3 < 1000, re xtreg Polity3 tIL5 if P3L5 < 1000 & tIL5 < 1000, re xtreg Polity3 tIL10 if P3L10 < 1000 & tIL10 < 1000, re xtreg Polity3 tIL20 if P3L20 < 1000 & tIL20 < 1000, re xtreg Polity3 tIL20 P3L20 if P3L20 < 1000 & tIL20 < 1000, re /********************************************************************/ /* 4. Merge with country_codes & save */ /********************************************************************/ merge ccode using "C:\STATADATA\country_codes.dta", drop if _merge==2 replace country=alpha if country=="" drop vnum vabb _mergeV presence presence2 pwtorder pwt61order pwtcode pwtname politycodes wdi2002 wa alpha _merge summarize sort ccode year tsset ccode year, yearly compress save ir127a, replace *drop if nP*nWT*nV!=1 *drop presence presence2 presence3 _merge drop if year<1945 tsset ccode year, yearly compress save ir127b, replace exit order ccode isocode country name year nP nWT nV exit