*Løsninger til opgaver i "Guide til Stata" *Numrene refererer til numrene i teksten ***************************************************** *********************KAPITEL 2*********************** ***************************************************** *1: Åbner data use "C:\temp\Stata\DifferentGoods.dta", clear *2: beskriver alle variable describe *3: kodebog for alle variable codebook *4: slår -more- fra set more off, permanently *5:kodebog igen codebook *11: describe igen describe *12: noter for alle variable notes *13: Kodebog med noter codebook, notes *14: Åbner databrowseren browse *16: oplistning af data i tabel list *17: tabel for tillid, genbrug og økologi list country gst recycle organic *18: *a) sorterer data stigende efter gst sort gst *b) oplister igen list country gst recycle organic *det fremgår af datamatricen og af tabellen at data nu er sorteret anderledes end før *c) sorterer efter gdpcap sort gdpcap *c) oplister landenavn og gdp per indbygger list country gdpcap *det fremgår at Romania er det fattigste land i datasættet *19: *sorterer gdpcap faldende gsort -gdpcap *opliste igen list country gdpcap *det fremgår at Luxembourg er det rigste land i datasættet *20: sorterer stigende efter gdpcap sort gdpcap *oplister de 10 første/fattigste lande i datasættet list country gdpcap eu in 1/10 *21: oplister EU-landenes BNP per indbygger (og EU-medlemskab) list country gdp eu if eu==1 *22: oplister hhv ikke eu-lande og eu-landes BNP per indbygger by eu,sort: list country gdpcap eu *23: frekvensfordeling for elecsystem tabulate elecsystem *24: frekvensfordeling for elecsystem inkl missing tabulate elecsystem, missing *26: frekvensfordeling for elecsystem uden labels tabulate elecsystem, missing nolabel *27: codebook for elecsystem codebook elecsystem *28: frekvensfordeling for 3 variable tab1 eu seaaccess recycle *29: inspect for alle variable inspect *30: histogram for recycle histogram recycle, frequency *32: summarize for alle variable med 1. option (standard display) summarize *summarize for alle variable med 2. option (additional...) summarize, detail *33: gennemsnitligt BNP for lande med havadgang, hhv. uden by seaaccess, sort : summarize gdpcap *34: Gentager nogle kommandoer *a) summarize, detail *b) inspect *c) histogram for organic histogram organic, frequency *d) kodebog codebook *36: *a) overfører use- og describe-kommandoen use "C:\temp\Stata\DifferentGoods.dta", clear describe *b) Histogram for gdpcap histogram gdpcap, frequency *40:gemmer data med navnet DifferentGoodsKap2.dta save "C:\temp\Stata\DifferentGoodsKap2.dta" *41: a) hjælpen til codebook help codebook *42: hjælpefilen for describe help describe *43: søger efter hjælp til notes search notes ***************************************************** *********************KAPITEL 3*********************** ***************************************************** *44: åbner data fra forrige kapitel use "C:\temp\Stata\DifferentGoodsKap2.dta", clear *46: tildeler labels til urban og female label variable urban "Urbanization" label variable female "Females' share of population" *47b): tildeler labels til education label variable education "Tertiary education (share)" *48: *a) tildeler note til frac notes frac: Ethnic fractionalization (late 1990s) from Alesina et al. 2003 *b) label label variable frac "Ethnic fractionalization" *c) Redigerer noten notes replace education in 1: Share of population above 25 with a tertiary education (in 2001) *49: viser eksisterende value labels label list *50: tilknytter noyesLB til eu label values eu noyesLB *51: opretter random1, random2, text og welfarestate generate random1 = . //a) generate random2 = . //a) generate str text = "" //b) generate byte welfarestate = . //c) browse //d) codebook random1 random2 text welfarestate //d) *52: indtastning af tal og tekst replace random1 = 1.1 in 1 replace random1 = 2.2 in 2 replace random2 = 1.2 in 1 replace random2 = 1.3 in 2 replace text = "test" in 1 replace text = "retest" in 2 *53: oprettlese af value label med navnet welfareLB label define welfareLB 1 "Conservative" 2 "Liberal" 3 "Socialist" //a) label values welfarestate welfareLB // b) tilknytter labelen til variablen welfarestate describe welfarestate // c) Checker med describe *56: *a): indtater data replace welfarestate = 1 in 1 replace welfarestate = 1 in 2 replace welfarestate = 1 in 7 replace welfarestate = 1 in 8 replace welfarestate = 1 in 9 replace welfarestate = 1 in 14 replace welfarestate = 1 in 26 replace welfarestate = 1 in 17 replace welfarestate = 2 in 13 replace welfarestate = 2 in 27 replace welfarestate = 3 in 5 replace welfarestate = 3 in 18 replace welfarestate = 3 in 25 *b) tildeler varibel-label label variable welfarestate "Welfare state regime" *c): tildeler note notes welfarestate: Data are entered 15.01.2011. Source: Esping-Andersen, The three world of welfare capitalism, 1990. *56: sletter variable drop random1 random2 text *57: summarize for gst, git, pm4 og gdpcap ; bestilt med bindestreg summarize gst-gdpcap *58: gemmer datasættet med et nyt navn save "C:\temp\Stata\DifferentGoodsKap3.dta" *59: keep if welfarestate == 3 //a) sletter observationer der ikke er social democratic summarize welfarestate //b) tjekker med summarize ***************************************************** *********************KAPITEL 4*********************** ***************************************************** *60: use "C:\temp\Stata\Valg05.dta", clear // a+b) åbner data describe //c) beskriver data browse //d) åbner datamatricen *61: søger efter oplysninger lookfor sympati socialdemo //sympati for socialdemokratiet *62: generate køn_gen = v310 //a) kopierer v310 til køn_gen codebook v310 køn_gen //b) kodebog *63: clonevar køn = v310 //a) kopierer med clonevar codebook v310 køn* //b) kodebog drop køn_gen //c) sletter køn_gen *64: flytter køn til starten af datasættet order køn *65: åbner dialogboksen til codebook codebook *66: Dialogboks for recode findes under Data > Create or change data > Other variable-transformation commands > Recode categorical variable *67: kodebog for variablen køn codebook køn *68: omkoder køn til mand recode køn (2=0), generate(mand) *69 kodebog for nydannede variabel (mand) codebook mand *70: oplister de 15 første observationers værdier på køn og mand list køn mand in 1/15 *71: list uden value-labels list køn mand in 1/15, nolabel *72: krydstabulering mellem køn og mand tabulate køn mand *73: drop mand //a) sletter mand recode køn (1=1 "Mand") (2=0 "Kvinde"), generate(mand) label(KvMaLB) //b) omkoder og tildeler labels codebook mand //c) *74: codebook v18, tabulate (20) //a) viser oplysninger om v18 recode v18 (1 2 5 7 10 = 0 "Rød") (3 4 6 8 9=1 "Blå") (nonmissing=.), generate(blå) label(RødBlåLB) //a+b) omkoder og tildeler value-label label variable blå "Rød/Blå stemme" //c)tildeler variabel-label list v18 blå in 1/25 //d)tjekker resultatet tabulate v18 blå //e) tjekker omkodning med krydstabularing tabulate blå //f)undersøger stemmefordelingen *75: codebook v18, tabulate (20) //a) viser oplysninger om v18 recode v18 (1=1 "A: Socialdemokraterne") (2= 2 "B: Det Radikale Venstre" ) /// a+b)omkoder og tildeler labels (3=3 "C: De Konservative") (4= 4 "D: Centrumdemokraterne") /// omkoder og tildeler labels (5= 5 "F: Socialistisk Folkeparti") (6= 6 "K: Kristendemokraterne") ///omkoder og tildeler labels (7= 7 "M: Minoritetspartiet") (8= 8 "O: Dansk Folkeparti") ///omkoder og tildeler labels (9 = 9 "V: Venstre") (10 = 10 "Ø: Enhedslisten") (11=11 "Kandidat uden for partierne") ///omkoder og tildeler labels (80 = .a "Husker ikke parti") (81 = .b "81 Vil ikke svare om parti") (82 = .c "Stemte blankt") ///omkoder og tildeler labels (90 = .d "Stemte ikke") (91 = .e "Ikke stemmeret") (92 = .f "Andre svar"), generate(partivalg) label(PativalgLB) //omkoder og tildeler labels tabulate v18 partivalg, missing //c) tjekker med krydstab; missing gør at missing værdierne også vises tabulate partivalg //d) viser stemmefordelingen i stikprøven *76: codebook v351, tabulate (20) //a) viser oplysninger om v351 recode v351 (14 = .a) (15=.b) (16=.c) (17=.d) (18=.e) (19=.f), generate(partivalgnu) //a)omkoder værdierne 13-19 til missing label copy v351 partivalgnuLB //b) kopierer value-label til partivalgnuLB label values partivalgnu partivalgnuLB //c)tilknytter labelen til variablen *77: label define partivalgnuLB 1 "A: Socialdemokraterne" 2 "B: Det Radikale Venstre" 3 "C: De Konservative" /// redigerer value-labelen 4 "D: Centrumdemokraterne" 5 "F: Socialistisk Folkeparti" 6 "K: Kristendemokraterne" 7 "M: Minoritetspartiet" /// 8 "O: Dansk Folkeparti" 9 "V: Venstre" 10 "Z: Fremskridtspartiet" 11 "Ø: Enhedslisten" 13 "Andet parti" /// .a "Stemme blankt" .b "Ved ikke" .c "Vil ikke svare om parti" .d "Ville ikke stemme" .e "Ikke stemmeret" .f "Andet svar", replace //a) tabulate v351 partivalgnu,missing //b) Krydstabulatring inkl. missing *78: codebook v186 v193 v222 v224 //a) inspicerer variablene med kodebogen recode v186 v222 v224 (1=5) (2=4) (4=2) (5=1) (8=.), generate (fremmedskep1 fremmedskep3 fremmedskep4) //b) vender disse variable om og sorterer missing fra recode v193 (8=.), generate(fremmedskep2) //b) sorterer missing fra tab1 fremmed* //c)inspicerer fordelingerne; det fremgår at v222 om udvisning fremkalder flest fremmedskeptiske svar *79: codebook v135-v142 //inspecerer variablene recode v135-v142 (1=0) (2=1) (3=.a) (8=.b), prefix(påstand_) //omkoder til ny_v136 osv *80: omkoder alder til en 4 delt variabel recode v357 (65/100 = 4) (50/65 = 3) (30/50=2) (1/30=1),generate(alder4kat) *81: tabulate v357 //a)frekvensfordeling for alder tabulate partivalgnu if v357 == 17,missing //b)frekvensfordelig for partivalg nu for 17-årrige inkl. missing *82: clonevar partivalgnu20 = partivalgnu //kopierer til ny variabel replace partivalgnu20 = .e if v357 < 20 //omkoder for respondenter under 20 år *83: codebook mand v357 //kodebog for variable clonevar partivalgnu_kønalder = partivalgnu //opretter partivalgnu_kønalder som kopi af partivalgnu replace partivalgnu_kønalder = .e if mand==1 & v357 <30 //erstatter med .e hvis respondenten er mand og under 30 år *84 tabulate partivalg, generate(parti) //a) list partivalg parti1-parti11 in 1/20 //b *88: Gemmer datasættet med et nyt navn save "C:\temp\Stata\Valg05Kap4.dta" ***************************************************** *********************KAPITEL 5*********************** ***************************************************** *89: starter log capture log close log using "C:\temp\Stata\Kap5", text replace *90: Åbner datasættet use "C:\temp\Stata\Valg05Kap4.dta", clear *91: tabulate partivalgnu // frekvensfordeling for partivalgnu *Venstre ville formentligt have fået flest stemmer *92: lagkagediagram i farver graph pie, over(partivalg) note(Kilde: Valgundersøgelsen 2005) //a) lagkagediagram i farver graph pie, over(partivalg) sort(blå) note(Kilde: Valgundersøgelsen 2005) //b) sorteret efter regeringspræference *93: histogram for partivalg histogram partivalg, discrete percent *94: histogram med afstand histogram partivalg, discrete percent gap(10) *95: Ændrer Y-label og indsætter linje ved 2% histogram partivalg, discrete percent gap(10) ytitle(Stemmefordeling (%)) yline(2, lpattern(dash)) *96: Inkluderer partinavne histogram partivalg, discrete percent gap(10) ytitle(Stemmefordeling (%)) yline(2, lpattern(dash)) xlabel(, valuelabel) *97: Inkluderer alle partinavne histogram partivalg, discrete percent gap(10) ytitle(Stemmefordeling (%)) yline(2, lpattern(dash)) xlabel(1(1)11, valuelabel) *98: Drejer labels histogram partivalg, discrete percent gap(10) ytitle(Stemmefordeling (%)) yline(2, lpattern(dash)) xlabel(1(1)11, angle(forty_five) valuelabel) *99: Fjerner label på X-aksen og tilføjer noter histogram partivalg, discrete percent gap(10) ytitle(Stemmefordeling (%)) yline(2, lpattern(dash)) xtitle("") xlabel(1(1)11, angle(forty_five) valuelabel) note(Stiplede linje indikerer spæregrænsen (2%); Kilde: Valgundersøgelsen 2005) *100: opsplitter forrige kommando så den er nemmere at læse histogram partivalg, discrete percent gap(10) ytitle(Stemmefordeling (%)) yline(2, lpattern(dash)) xtitle("") /// xlabel(1(1)11, angle(forty_five) valuelabel) note(Den stiplede linje indikerer spærregrænsen (2%); Kilde: Valgundersøgelsen 2005) *101: danner dummyvariable på baggrund af partivalgnu tabulate partivalgnu, generate(partinu) *102: ci partinu3, binomial //a) 95% konfidensinterval for andelen at observationerne der angivar at de ville stemme på de konservative ci partinu3, binomial level(99) //a) som ovenfor men med 99% bitest partinu6==0.02 //b)Test for om andelen af observationene der angiver at de vil stemme på Kristendemokrateren er lig med 2% *Krinstendemokraterne ville sandsynligvis ikke komme over spærregrænsen *103: codebook v250, tabulate(20) //a)kodebog for v250; tabulate(20) bevirker at der vises en frekvensfordeling selvom der er mere end 10 unikke værdier recode v250 (88=.), generate(vhplacering) //b)omkoder til vhplacering, hvor ved ikke er kodet til missing summarize vhplacering, detail //c) summarize for vhplacering *104: Deskriptive statistikker for vhplacering og v357 tabstat vhplacering v357, statistics( mean iqr ) columns(variables) *105: histogram vhplacering // histogram for vhplacering, Kontinuerlig specifikation histogram vhplacering, discrete // histogram for vhplacering, diskret specifikation *104b: 4 histogrammer med forskellige målestokker på y-aksen histogram vhplacering, discrete histogram vhplacering, discrete fraction histogram vhplacering, discrete frequency histogram vhplacering, discrete percent *107: ci vhplacering //a) 95%konfidensinterval ci vhplacering,level(99) //a) 99%konfidensinterval ttest vhplacering == 5 //b) t-test *108: gemmer log close //b) lukker log-filen save "C:\temp\Stata\Valg05Kap5.dta" //c) gemmer *109: gentager tidligere opgave tabulate partivalgnu ***************************************************** *********************KAPITEL 6*********************** ***************************************************** *113: capture log close log using "C:\temp\Stata\Kap6", text replace use "C:\temp\Stata\Valg05Kap5.dta", clear //Åbner datasættet db tabulate2 // dialogboksen til krydstabulering *114: codebook v178 recode v178 (1=0 skatlet) (2=1 service) (8=.), generate(service) label(skatletserviceLB) //omkoder og tildeler value-labels lab var service "Prioritering af service (vs. skattelettelse)" //tildeler variabel labels tabulate service mand //krydstabulering *115: Krydstabularing med relative (kolonne)frekevnser tabulate service mand, column *116: prtest service, by(mand) //samenligning af to gruppers andele *117: ttest vhplacering, by(mand) //sammenligning af to gruppers gennemsnit *118: varianshomogenitet? robvar vhplacering, by(mand) *119: codebook v67, tabulate(20) //codebog for v67 med tabulering recode v67 (18/19 = .), generate(symsocdem) //omkoder robvar symsocdem,by(v365) //test af varianshomogenitet ttest symsocdem,by(v365) *120: Flere gruppers gennemsnit mean vhplacering, over(priavis) //gennemsnit på vhplacering for grupperne defineret af avisvalg *121: ANOVA oneway vhplacering priavis,tabulate *122: som før men med bonferroni oneway vhplacering priavis,tabulate bonferroni *124: codebook v6,tabulate(20) //codebook med tabularing tab v6 blå,chi2 row //chi2-test tab blå v6,chi2 column //anden opsætning *125: codebook v337 //a) undersøger variablen recode v337 (8=.),generate(urban) //a) omkoder "ved ikke" til missing lab var urban "Urbanisering" //a) tildeler variabel label tab blå urban, gamma taub //b) gamma og taub for blå/urban tab2 fremmedskep*, gamma taub //d)krydstabulering mellem flere variable *126: pearsons r pwcorr vhplacering v357, sig *127: scatterplot scatter vhplacering v357 *128: sunflower-plot sunflower vhplacering v357 * 129: save "C:\temp\Stata\Valg05Kap6.dta" // gemmer datasættet use "C:\temp\Stata\DifferentGoodsKap3.dta", clear //åbner datasæt *130: Scatterplot scatter gst git //a) scatterplot pwcorr gst git,sig //b) korrelation og P-værdi *131: Scatterplot via menuen twoway (scatter gst git) *132: Samme plot, mem med landenavne twoway (scatter gst git, mlabel(country)) *133: twoway (scatter gst git, mlabel(country)) (lfit gst git) //Som før men med OLS linje twoway (scatter gst frac, mlabel(country)) (lfit gst frac) //nye variable *134: med navne twoway (scatter gst git, mlabel(country)) (lfit gst git) , name(scat_gst_git) twoway (scatter gst frac, mlabel(country)) (lfit gst frac) , name(scat_gst_frac) *135: set autotabgraphs on, permanently //fanevisning af grafer twoway (scatter gst git, mlabel(country)) (lfit gst git) , name(scat_gst_git, replace) twoway (scatter gst frac, mlabel(country)) (lfit gst frac) , name(scat_gst_frac, replace) *137: konfidenintervaller med tre stikprøvestørrelser cii 250 0.52 cii 500 0.52 cii 1000 0.52 *138: tabi 30 20 \ 40 10\45 5, gamma *139: h tabulate twoway *140: use "C:\temp\Stata\Valg05Kap6.dta", clear //a tab service mand, gamma taub //b return list //c *141: display r(gamma) *142: display r(gamma)/r(ase_gam) *143: display normal(r(taub)/r(ase_taub)) log close ***************************************************** *********************KAPITEL 7*********************** ***************************************************** *144: log + Åbner data capture log close log using "C:\temp\Stata\Kap7", text replace use "C:\temp\Stata\Valg05Kap6.dta", clear *145: codebook rc_v182 rc_v186 rc_v193 rc_v222 rc_v224 //kodebog *146: recode rc_v182 (1=1) (2=3) (3=5) //omkoder rc_v182 til ny variationsbredde recode rc_v186 rc_v222 rc_v224 (1=5) (2=4) (3=3) (4=2) (5=1) //vender items *147: item-item-analyser med gamma og Pearsons r tab2 rc_*,gamma pwcorr rc_* *148: reliabilitetstest alpha rc_* , item casewise *149: danner tempindeks, kun observationer med mindst to gyldige værdier medtages alpha rc_* , min(2) generate(tempindeks) item *150: su tempindeks //summarize for indekset generate fremmedskepticisme = ((tempindeks-1)/(5-1))*100 //reskalerer *151: sum fremmedskepticisme,detail //summarize for indekset histogram fremmedskepticisme //histogram drop tempindeks //sletter tempindeks *152a: Linjerne fra opgave 152 foreach var in v186 v193{ recode `var' (8=.a),generate(`var'mis) sum `var'mis replace `var'mis = r(mean) if `var'mis ==.a } *153: sletter og gemmer drop v186mis v193mis save "C:\temp\Stata\Valg05Kap7.dta" log close ***************************************************** *********************KAPITEL 8*********************** ***************************************************** *154: startet log+åbner data capture log close log using "C:\temp\Stata\Kap8", text replace use "C:\temp\Stata\Valg05Kap7.dta", clear *155: bivariat regression reg fremmedskepticisme v357 *156: reg fremmedskepticisme v357, beta // som før men med beta pwcorr fremmedskepticisme v357, sig obs //pearsons r *157: reg fremmedskepticisme mand //mand og kvinders fremmedskepticisme med regress ttest fremmedskepticisme, by(mand) // samme med ttest *158: tab urban, generate(durban) //danner dummevariable reg fremmedskepticisme durban2-durban5 //sammenhængen mellem urbanisering og fremmedskepticisme med regress + dummyvariable *159: reg fremmedskepticisme i.urban //som foregående men med faktorvariable oneway fremmedskepticisme urban, bonferroni //som forgående men med oneway *160: codebook urban reg fremmedskepticisme ib(3).urban //som foregående men med kategori 3 som reference *161: codebook v312,tabulate(25) recode v312 (1=1) (nonmis=0),generate(indfødt) *162: reg fremmedskepticisme v357 indfødt mand *163: reg fremmedskepticisme v357 indfødt mand,robust *164: reg fremmedskepticisme v357 indfødt mand i.urban *165: reg fremmedskepticisme v357 indfødt mand i.urban reg fremmedskepticisme v357 indfødt mand if e(sample)==1 *166: save "C:\temp\Stata\Valg05Kap8.dta" *167: Søger efter estout-odulet search estout, all *Nedenstående kommando (uden stjerne) installerer modulet *net install estout.pkg *168: eststo clear //tømmer hukommelsen reg fremmedskepticisme v357 indfødt mand //regression1 eststo //gemme resultaterne reg fremmedskepticisme v357 indfødt mand i.urban //regression2 eststo //gemmer resultaterne esttab //viser resultaterne fra alle estimtioner i hukommelsen *169: gemme tabellen i en txt.fil med navnet fremmedskepticisme esttab using frememdskepticisme, r2 ar2 se replace log close ***************************************************** *********************KAPITEL 9*********************** ***************************************************** *170: capture log close log using "C:\temp\Stata\Kap9", text replace //a) starter log use "C:\temp\Stata\Valg05Kap8",clear //a) åbner data reg fremmedskepticisme v357 mand //b) regression acprplot v357 , lowess //c) acpr-plot for alder *171: leverage-versus-residual-squared-plot lvr2plot *172: predict tempcooksd, cooksd //beregner cook's d display 4/e(N) //udregner kritisk værdi summarize tempcooksd if tempcooksd >(4/e(N)) //opsummerer observationer med kritiske værdier *173: avplots avplots *174: DFBETA predict tempdfbeta, dfbeta(v357) summarize tempdfbeta if abs(tempdfbeta) > 2/sqrt(e(N)) *175: histogram *beregning af residualer, der gemmes i variablen 'tempresid': predict tempresid, residuals *Histogram, inkl. normalfordelingskurve, for residualernes fordeling: histogram tempresid, normal *176: Q-Q-plot qnorm tempresid *177: rvfplot,yline(0) //rvfplot rvpplot v357, yline(0) rvpplot mand, yline(0) *178: white's test estat imtest,white *179: VIF vif *180: drop temp* //rydder op save "C:\temp\Stata\Valg05Kap9.dta" //gemmer datasæt log close ***************************************************** *********************KAPITEL 10********************** ***************************************************** *181: Åbner data capture log close log using "C:\temp\Stata\Kap10", text replace use "C:\temp\Stata\GrowthKap10.dta", clear *182: regressionsmodellen reg meangrowth aidcap goodgov initial forinvest eduyears,robust //uden interaktion reg meangrowth aidcap goodgov c.aidcap#c.goodgov initial forinvest eduyears,robust //med interaktion *183: statistikker for variblene fra forrige estimation baseret på estimationssamplen estat summarize *184: Begegning af marginale effekter margins, dydx(aidcap) at(goodgov=(-1.8(0.2)1.8)) *185: marginsplot marginsplot *186: Illustration af Figur 10.5 marginsplot, /// starter kommandoen marginsplot recast(line) /// beder om en linje fremfor punkter recastci(rline) /// specificerer ar konfidensintervallet også skal tegner med linjer ciopts(lpattern(dash)) /// specificerer at konfidensintervallet skal tegnes med stiplede linjer yline(0) /// indsætter en linje ved y=0 xlabel(#5) /// nedjusterer antallet af tal på x-aksen legend(on /// beder om at få vist forklaring label(1 "95% konfidensinterval") /// specificerer forklaringen label(2 "Marginal effekt af ulandshjælp")) /// specificerer forklaringen xtitle(Institutionel kvalitetet) /// indsætter alternativ titel på x-aksen ytitle("Marginal effekt af ulandshjælp") /// indsætter alternativ titel på y-aksen title("") /// fjerner den overordnde titel scheme(s2mono) // optimerer til sort/hvid-visning *187: regress, coeflegend *188: display _b[aidcap] *************************************************** *kommandoen der producerer Figur 10.4 reg meangrowth aidcap goodgov c.aidcap#c.goodgov initial forinvest eduyears,robust margins, dydx(aidcap) at((p10) goodgov) at((p50) goodgov) at((p90) goodgov) foreach var in initial forinvest eduyears { sum `var' if e(sample)==1 local `var' = r(mean) } sum goodgov if e(sample)==1, detail local pti = r(p10) local pfemti = r(p50) local pniti = r(p90) local level `pti' `pfemti' `pniti' foreach lvl in `level' { local plots `plots' /// (function y = _b[_cons] /// + _b[aidcap]*x + _b[goodgov]*`lvl' + _b[c.aidcap#c.goodgov]*x*`lvl' /// + _b[initial]*`initial' + _b[forinvest]*`forinvest' + _b[eduyears]*`eduyears' /// , range(1 130)) } twoway `plots', /// ytitle("Gennemsnitlig årlig vækst (%)") /// xtitle("Ulandsbistand per capita ($)") /// legend (order(3 2 1) /// label(1 "Lav institutionel kvalitet **") /// label(2 "Middel institutionel kvalitet **") /// label(3 "Høj institutionel kvalitet **")) /// note("**: P< 0,05; Lav, Middel, Høj instititionel kvalitet: 10. 50. og 90. percentil" /// "Kontrolvariable er sat til deres gennensnit") /// scheme(s2mono) ******************************************************** *189: reg meangrowth aidcap goodgov c.aidcap#c.goodgov initial forinvest eduyears,robust margins, dydx(aidcap) at((p25) goodgov) at((p50) goodgov) at((p75) goodgov) foreach var in initial forinvest eduyears { sum `var' if e(sample)==1 local `var' = r(mean) } sum goodgov if e(sample)==1, detail local ptotifem = r(p25) local pfemti = r(p50) local psyvtifem = r(p75) local level `ptotifem' `pfemti' `psyvtifem' foreach lvl in `level' { local plots `plots' /// (function y = _b[_cons] /// + _b[aidcap]*x + _b[goodgov]*`lvl' + _b[c.aidcap#c.goodgov]*x*`lvl' /// + _b[initial]*`initial' + _b[forinvest]*`forinvest' + _b[eduyears]*`eduyears' /// , range(1 130)) } twoway `plots', /// ytitle("Gennemsnitlig årlig vækst (%)") /// xtitle("Ulandsbistand per capita ($)") /// legend (order(3 2 1) /// label(1 "Lav institutionel kvalitet **") /// label(2 "Middel institutionel kvalitet **") /// label(3 "Høj institutionel kvalitet")) /// note("**: P< 0,05; Lav, Middel, Høj instititionel kvalitet: 25. 50. og 75. percentil" /// "Kontrolvariable er sat til deres gennensnit") /// scheme(s2mono) log close