Introduction to Stata 17 871 Spring 2012 The

  • Slides: 34
Download presentation
Introduction to Stata 17. 871 Spring 2012

Introduction to Stata 17. 871 Spring 2012

The role of statistical packages in research • Obvious answer – Manage data –

The role of statistical packages in research • Obvious answer – Manage data – Carry out appropriate statistical tests – Assist in displaying data • Less obvious answer – Channel the type of research you are likely to do • Limitations as to variables and cases • Types of analysis is sometimes guided by choice of package

Analysis -> Packages • Baby exercises – Minitab, spreadsheets • Time series – TSP

Analysis -> Packages • Baby exercises – Minitab, spreadsheets • Time series – TSP • Cross-sectional – SPSS, SAS • Time series & cross-sectional – Stata, R

Logic of quant research in this class

Logic of quant research in this class

Logic of data setup: V 1 Obs 2 … Obsi V 2 … Vj

Logic of data setup: V 1 Obs 2 … Obsi V 2 … Vj

Example, VRS Data HRHHID 199960521980910 160916068405549 941159210626002 GESTCEN 63 63 63 PES 1 2

Example, VRS Data HRHHID 199960521980910 160916068405549 941159210626002 GESTCEN 63 63 63 PES 1 2 2 2 PES 8 4 -3 6 6 6

Example, House Elections

Example, House Elections

Using Stata to Analyze Data in Matrix Form • Question: Did Ron Paul do

Using Stata to Analyze Data in Matrix Form • Question: Did Ron Paul do better in Iowa in 2012, compared to 2008 in counties with college students? • Data sources: – 2008: Des Moines Register web site, (http: //caucuses. desmoinesregister. com/data/iowacaucus/caucus-history-gop/) – 2012: Iowa Republican Party, Google Doc (https: //www. google. com/fusiontables/Data. Source? d srcid=2475248)

2008

2008

2012

2012

Switch over to Stata run-through

Switch over to Stata run-through

Return from Stata run-through • Why would you use different input commands?

Return from Stata run-through • Why would you use different input commands?

insheet • Data is output from a spreadsheet into “csv” or “comma-delimited” format •

insheet • Data is output from a spreadsheet into “csv” or “comma-delimited” format • Data is a simple I x J matrix, and all the variables are separated either by a tab or comma • Stata is now smart enough to figure out that the first line of the file contains the variable names

insheet Assume the following file was created by outputting a file from Excel in

insheet Assume the following file was created by outputting a file from Excel in csv format: HRHHID 199960521980910 160916068405549 941159210626002 insheet using filename GESTCEN 63 63 63 PES 1 2 2 2 PES 8 4 -3 6 6 6

infile • Data is not in Stata format, is in an ASCII file, but

infile • Data is not in Stata format, is in an ASCII file, but is not separated only by a tab or comma (e. g. , by a space)

insheet Assume the following file was created using an ASCII text editor (e. g.

insheet Assume the following file was created using an ASCII text editor (e. g. , EMACS), and that spaces separate the variables: 199960521980910 160916068405549 941159210626002 63 63 63 2 2 2 4 -3 6 6 6 infile HRHHID GESTCEN PES 1 PES 8 using filename Or infile str HRHHID GESTCEN PES 1 PES 8 using filename

infix • Data is in an ASCII file, but you cannot rely on spaces,

infix • Data is in an ASCII file, but you cannot rely on spaces, commas, or other standard “delimeters” to separate variables • Datasets may have observations on more than one line

infix Assume the following file was created using an ASCII text editor: 1 2

infix Assume the following file was created using an ASCII text editor: 1 2 12345678901 ----------19996052198091063 2 4 16091606840554963 2 -3 94115921062600263 2 6 Handy label, not in dataset Dataset infix HRHHID 1 -15 GESTCEN 16 -17 PES 1 18 -19 PES 8 20 -21 using filename Or infile str 15 HRHHID 1 -15 GESTCEN 16 -17 PES 1 18 -19 PES 8 20 -21 using filename

House Roll Call votes in the 27 th Cong. 01 R 327031200290003401 ADAMS 02

House Roll Call votes in the 27 th Cong. 01 R 327031200290003401 ADAMS 02 R 327031200290003401 ADAMS 03 R 327031200290003401 ADAMS 04 R 327031200290003401 ADAMS 05 R 327031200290003401 ADAMS 06 R 327031200290003401 ADAMS 07 R 327031200290003401 ADAMS 08 R 327031200290003401 ADAMS 09 R 327031200290003401 ADAMS 10 R 327031200290003401 ADAMS 11 R 327031200290003401 ADAMS 12 R 327031200290003401 ADAMS 13 R 327031200290003401 ADAMS 14 R 327031200290003401 ADAMS 15 R 327031200290003401 ADAMS 16 R 327031200290003401 ADAMS 17 R 327031200290003401 ADAMS 18 R 327031200290003401 ADAMS 19 R 327031200290003401 ADAMS 20 R 327031200290003401 ADAMS 21 R 327031200290003401 ADAMS 01 R 327449800320009111 ALFORD 02 R 327449800320009111 ALFORD 16555561661111222226261116611966116116116666 66616111666116666166111166116116191611666666 6611666111166666119166116991611161616116661616611916991111661611661616616611611 16666661611161916616116166661611116666161111 166666161116161166111111661666661126611661666666 6966616166666111691116111111611666111111616 11916666166666611166666999991161661169999161 6666161116666966161611166661611119 6116166161666116116111116161119919966 116191666161161166696616111616661161166911691666 61116669966161666116616111616611116611666661 611661616166661661696166661166661666611 11616111116116661161116666166616616616661166 6116161111611611111616611166111666166 161116619116666616611616166661966661611616616611161611116116666166661111616666611111 1116669916166611116616166116116161666166111611611611161666666111911611 61661661611916166196666116616661111161 65555599661691616555525651111611611199199999 916916661169611661661161999911611611111161169999

1 2 3 4 5 6 7 8 1234567890123456789012345678901234567890 01 R 327031200290003401 ADAMS 16555561661111222226261116611966116116116666

1 2 3 4 5 6 7 8 1234567890123456789012345678901234567890 01 R 327031200290003401 ADAMS 16555561661111222226261116611966116116116666 02 R 327031200290003401 ADAMS 66616111666116666166111166116116191611666666 03 R 327031200290003401 ADAMS 661166611116666611916611699161116161611 04 R 327031200290003401 ADAMS 16116661616611916991111661611661616616611611 05 R 327031200290003401 ADAMS 16666661611161916616116166661611116666161111 06 R 327031200290003401 ADAMS 166666161116161166111111661666661126611661666666 07 R 327031200290003401 ADAMS 6966616166666111691116111111611666111111616 08 R 327031200290003401 ADAMS 11916666166666611166666999991161661169999161 09 R 327031200290003401 ADAMS 6666161116666966161611166661611119 10 R 327031200290003401 ADAMS 6116166161666116116111116161119919966 11 R 327031200290003401 ADAMS 116191666161161166696616111616661161166911691666 12 R 327031200290003401 ADAMS 61116669966161666116616111616611116611666661 13 R 327031200290003401 ADAMS 611661616166661661696166661166661666611 14 R 327031200290003401 ADAMS 11616111116116661161116666166616616616661166 15 R 327031200290003401 ADAMS 6116161111611611111616611166111666166 16 R 327031200290003401 ADAMS 161116619116666616611616166661966661611616616611 17 R 327031200290003401 ADAMS 6611161611116116666166661111616666611111 18 R 327031200290003401 ADAMS 1116669916166611116616166116116161666 19 R 327031200290003401 ADAMS 1666166111611611611161666666111911611 20 R 327031200290003401 ADAMS 6166166161191616619666611661666111116 21 R 327031200290003401 ADAMS 61111161 01 R 327449800320009111 ALFORD 65555599661691616555525651111611611199199999 02 R 327449800320009111 ALFORD 916916661169611661661161999911611611111161169999 VAR # 0004 WIDTH = 0002 MD=0 DK 01 COL 07 -08 H 27 STATE: . . . NEW ENGLAND. . . 01. CONNECTICUT 02. MAINE 03. MASSACHUSETTS BORDER STATES. . . 51. KENTUCKY 52. MARYLAND 53. OKLAHOMA

1 2 3 4 5 6 7 8 1234567890123456789012345678901234567890 01 R 327031200290003401 ADAMS 16555561661111222226261116611966116116116666

1 2 3 4 5 6 7 8 1234567890123456789012345678901234567890 01 R 327031200290003401 ADAMS 16555561661111222226261116611966116116116666 02 R 327031200290003401 ADAMS 66616111666116666166111166116116191611666666 03 R 327031200290003401 ADAMS 661166611116666611916611699161116161611 04 R 327031200290003401 ADAMS 16116661616611916991111661611661616616611611 05 R 327031200290003401 ADAMS 16666661611161916616116166661611116666161111 06 R 327031200290003401 ADAMS 166666161116161166111111661666661126611661666666 07 R 327031200290003401 ADAMS 6966616166666111691116111111611666111111616 08 R 327031200290003401 ADAMS 11916666166666611166666999991161661169999161 09 R 327031200290003401 ADAMS 6666161116666966161611166661611119 10 R 327031200290003401 ADAMS 6116166161666116116111116161119919966 11 R 327031200290003401 ADAMS 116191666161161166696616111616661161166911691666 12 R 327031200290003401 ADAMS 61116669966161666116616111616611116611666661 13 R 327031200290003401 ADAMS 611661616166661661696166661166661666611 14 R 327031200290003401 ADAMS 11616111116116661161116666166616616616661166 15 R 327031200290003401 ADAMS 6116161111611611111616611166111666166 16 R 327031200290003401 ADAMS 161116619116666616611616166661966661611616616611 17 R 327031200290003401 ADAMS 6611161611116116666166661111616666611111 18 R 327031200290003401 ADAMS 1116669916166611116616166116116161666 19 R 327031200290003401 ADAMS 1666166111611611611161666666111911611 20 R 327031200290003401 ADAMS 6166166161191616619666611661666111116 21 R 327031200290003401 ADAMS 61111161 01 R 327449800320009111 ALFORD 65555599661691616555525651111611611199199999 02 R 327449800320009111 ALFORD 916916661169611661661161999911611611111161169999 VAR # 0005 WIDTH = 0002 MD=0 DK 01 COL 09 -10 H 27 DISTRICT NUMBER: . . . . CODED BLANK FOR SENATE. AT-LARGE DISTRICTS ARE CODED 98, 97, 96, ACCORDING TO ALPHABETICAL ORDER OF NAMES OF OCCUPANTS. NO DISTINCTION BETWEEN THE VARIOUS KINDS OF AT-LARGE DISTRICTS IS MADE. DUE TO REPLACEMENTS WITHIN A CONGRESS, TWO MEMBERS MAY LEGITIMATELY HAVE THE SAME DISTRICT NUMBER WITHIN A STATE.

1 2 3 4 5 6 7 8 1234567890123456789012345678901234567890 01 R 327031200290003401 ADAMS 165555616

1 2 3 4 5 6 7 8 1234567890123456789012345678901234567890 01 R 327031200290003401 ADAMS 165555616 661661111222226261116611966116116116666 02 R 327031200290003401 ADAMS 66616111666116666166111166116116191611666666 03 R 327031200290003401 ADAMS 661166611116666611916611699161116161611 04 R 327031200290003401 ADAMS 16116661616611916991111661611661616616611611 05 R 327031200290003401 ADAMS 16666661611161916616116166661611116666161111 06 R 327031200290003401 ADAMS 166666161116161166111111661666661126611661666666 07 R 327031200290003401 ADAMS 6966616166666111691116111111611666111111616 08 R 327031200290003401 ADAMS 11916666166666611166666999991161661169999161 09 R 327031200290003401 ADAMS 6666161116666966161611166661611119 10 R 327031200290003401 ADAMS 6116166161666116116111116161119919966 11 R 327031200290003401 ADAMS 116191666161161166696616111616661161166911691666 12 R 327031200290003401 ADAMS 61116669966161666116616111616611116611666661 13 R 327031200290003401 ADAMS 611661616166661661696166661166661666611 14 R 327031200290003401 ADAMS 11616111116116661161116666166616616616661166 15 R 327031200290003401 ADAMS 6116161111611611111616611166111666166 16 R 327031200290003401 ADAMS 161116619116666616611616166661966661611616616611 17 R 327031200290003401 ADAMS 6611161611116116666166661111616666611111 18 R 327031200290003401 ADAMS 1116669916166611116616166116116161666 19 R 327031200290003401 ADAMS 1666166111611611611161666666111911611 20 R 327031200290003401 ADAMS 6166166161191616619666611661666111116 21 R 327031200290003401 ADAMS 61111161 01 R 327449800320009111 ALFORD 655555996 61691616555525651111611611199199999 02 R 327449800320009111 ALFORD 916916661169611661661161999911611611111161169999 VAR # 0020 SESSION 1 WIDTH = 0001 MD=0 DK 01 COL 42 -42 H 27 G-10 - -27 A J 27 -1 -39 JUNE 7, 1841 H 271004 Y=66 N=149 MALLORY, VA. TO ADJOURN, IN ORDER TO END DEBATE ON THE ADOPTION OF THE HOUSE RULES. ADOPTION OF THE RULES WOULD PREVENT RECEIVING ANY ABOLITION PETITIONS. (P. 27 -2)

Enter data yourselves

Enter data yourselves

Return again to Stata run-through

Return again to Stata run-through

merge command • Used when you want to add data to a preexisting data

merge command • Used when you want to add data to a preexisting data set, or you have more than one dataset that has all the variables you need for analysis. • Most important thing: each dataset must have (at least) one identifier that links observations, and allows merging. • Second thing: both datasets must be sorted on the common identifier(s)

Example: one-for-one match Election results, election_results. dta county cand 1 cand 2 A 10

Example: one-for-one match Election results, election_results. dta county cand 1 cand 2 A 10 20 30 B 40 50 60 C 70 80 90 Z 500 40 30 Demographics, demographics. dta county income educ catholic A 10, 000 . 2 . 3 B 40, 000 . 5 . 6 C 70, 000 . 8 . 9 Z 5, 000 . 95 . 3

merge command results • [assume both datasets have previously been sorted on county, by

merge command results • [assume both datasets have previously been sorted on county, by typing the command sort county] • use election_results. dta • merge county using demographics. dta OR • merge 1: 1 county using demographics. dta

Voila! county cand 1 cand 2 income educ catholic A 10 20 30 10,

Voila! county cand 1 cand 2 income educ catholic A 10 20 30 10, 000 . 2 . 3 B 40 50 60 40, 000 . 5 . 6 C 70 80 90 70, 000 . 8 . 9 Z 500 40 30 5, 000 . 95 . 3

many-to-one merge Demographic data, demographic_data. dta county_code town income education A Aville 50000 .

many-to-one merge Demographic data, demographic_data. dta county_code town income education A Aville 50000 . 3 A Bobville 60000 . 4 B Candiceville 70000 . 5 B Dogville 80000 . 5 C Catville 100000 . 5 County code mapping, county_code_mapping. dta county_code county_name A Adams B Brooks C Calhoun

merge command • [make same sorting assumptions as before] • use demographic_data. dta •

merge command • [make same sorting assumptions as before] • use demographic_data. dta • merge m: 1 county_code using county_code_mapping. dta

Voila! county_code town income education county_name A Aville 50000 . 3 Adams A Bobville

Voila! county_code town income education county_name A Aville 50000 . 3 Adams A Bobville 60000 . 4 Adams B Candiceville 70000 . 5 Brooks B Dogville 80000 . 5 Brooks C Catville 100000 . 5 Calhoun

collapse command county District. Name-en Adair - 1 NW ADAIR voters 46 Paul 7

collapse command county District. Name-en Adair - 1 NW ADAIR voters 46 Paul 7 Bachmann 4 Johnson Gingrich 11 Santorum 10 0 8 6 0 Adair - 2 NE STUART 51 8 5 3 15 1 0 0 6 13 0 Adair - 3 SW FONTANELLE 55 9 6 16 14 0 0 0 3 7 0 Adair - 4 SE ORIENT 50 4 6 6 15 0 0 0 13 6 0 Adair - 5 GF GREENFIELD 67 14 5 8 12 0 0 0 13 15 0 Adams – Carbon 28 7 0 5 12 0 0 0 3 1 0 Adams - Corning 1 A 19 7 0 1 6 0 0 0 4 1 0 Adams - Corning 1 B 3 3 0 0 0 0 0 Adams - Corning 2 A 9 2 0 0 0 0 5 0 0 Adams - Corning 2 B 8 5 1 0 0 0 0 1 0 Adams - Corning 3 A 12 4 0 0 0 6 2 0 Adams - Corning 3 B 19 9 0 1 6 0 0 0 1 2 0 Adams – Nodaway 10 1 1 5 0 0 2 1 0 Adams – Prescott 32 21 3 2 1 0 0 0 3 2 0 Adams – Quincy 22 7 0 2 8 0 0 0 3 2 0 Adams - SE Adams 38 8 4 6 13 0 0 0 5 2 0 Allamakee - FV/TL/HF CITY 28 7 0 6 9 0 0 0 6 0 0 Allamakee - LF/CN/LS/LS CITY 64 20 2 21 7 0 0 0 4 10 0 Allamakee - PC/LT/WV CITY 42 20 0 7 9 0 0 0 5 1 0 Allamakee - PO/FK 20 4 1 5 3 0 0 0 6 1 0 Allamakee - PV CITY 35 7 1 2 3 0 0 0 21 1 0 Allamakee - UC/IA/NA CITY 31 4 1 3 4 0 0 0 16 3 0 Allamakee - UP/MK/FC/JF/LL 122 53 2 18 14 0 0 0 28 7 0 Allamakee - WK 1 CITY 33 8 1 6 12 0 0 0 5 1 0 collapse (sum) voters-Cain, by(county) Huntsman Other Roemer Romney Perry Cain

county Adair Adams Allamakee voters Paul Bachmann Johnson Gingrich Santorum Huntsman Other Roemer Romney

county Adair Adams Allamakee voters Paul Bachmann Johnson Gingrich Santorum Huntsman Other Roemer Romney Perry Cain 269 200 518 42 74 157 26 9 18 0 0 0 44 24 82 66 47 77 1 0 0 0 0 43 32 155 47 14 28 0 0 0 Appanoose Audubon Benton 537 223 1042 77 41 202 25 17 66 0 0 0 71 32 121 174 54 290 1 0 5 0 0 1 12 0 0 87 48 184 90 31 168 0 0 4 Black Hawk Boone Bremer Buchanan 3642 1344 933 459 870 276 194 66 262 104 57 40 0 0 596 160 98 77 783 400 215 133 29 4 14 1 0 0 2 2 4 0 0 0 835 230 246 78 259 170 105 62 1 0 0 0 Buena Vista Butler Calhoun Carroll Cass Cedar 716 552 435 716 674 711 169 99 75 133 116 188 26 41 31 32 32 34 0 0 0 128 71 54 145 147 84 157 131 168 170 167 3 4 2 2 2 4 0 0 2 0 0 124 92 69 146 141 165 110 87 71 85 66 67 2 0 0 1 0 0 Cerro Gordo Cherokee 1571 537 304 95 100 20 0 0 235 78 345 155 5 0 1 0 0 0 408 126 170 63 2 0 Chickasaw Clarke Clayton Clinton Crawford 443 367 733 625 1384 437 142 98 150 205 295 72 14 42 40 28 62 22 0 0 0 53 46 137 72 149 84 72 51 165 122 354 101 3 1 4 1 9 0 0 2 2 0 0 0 0 0 85 65 149 116 437 93 74 62 75 81 73 64 0 0 5 0

Do-files • Do-files are the Stata scripting language to automate analysis. • Here is

Do-files • Do-files are the Stata scripting language to automate analysis. • Here is how the first five lines of the Iowa exercise would look in a do-file: #delimit; insheet using iowa_example_csv. dat; list; generate paulpct 08=paul 08/tvotes 08; generate paulpct 12=paul 12/tvotes 12;