|
|
|
||||
|
|
|||||
page 4, Figure 1.1
use "c:\vizdata\barley.dta", clear
drop rownames
reshape wide yield, i(variety site) j(year) (note: j = 1931 1932)
Data long -> wide ----------------------------------------------------------------------------- Number of obs. 120 -> 60 Number of variables 4 -> 4 j variable (2 values) year -> (dropped) xij variables: yield -> yield1931 yield1932 -----------------------------------------------------------------------------
graph dot yield1931 yield1932, /// over(variety, gap(20) label(labsize(vsmall))) /// by(site, compact cols(2) b1title("Barely Yield (bushels/acre)") /// note("")) ylabel(15(15)60) marker(1,msymbol(+)) marker(2,msymbol(Oh)) /// legend(label(1 1931) label(2 1932))
You can use the code below to have the graph ordered as it is in the text.
use "c:\vizdata\barley.dta", clear
drop rownames
reshape wide yield, i(variety site) j( year) (note: j = 1931 1932)
Data long -> wide ----------------------------------------------------------------------------- Number of obs. 120 -> 60 Number of variables 4 -> 4 j variable (2 values) year -> (dropped) xij variables: yield -> yield1931 yield1932 -----------------------------------------------------------------------------
gen site2 = 1 if site=="Waseca" (50 missing values generated)
replace site2 = 2 if site=="Crookston" (10 real changes made)
replace site2 = 3 if site=="Morris" (10 real changes made)
replace site2 = 4 if site=="University Farm" (10 real changes made)
replace site2 = 5 if site=="Duluth" (10 real changes made)
replace site2 = 6 if site=="Grand Rapids" (10 real changes made)
label define sitelbl 1 "Waseca" 2 "Crookston" 3 "Morris" /// 4 "University Farm" 5 "Duluth" 6 "Grand Rapids"
label values site2 sitelbl
gen var = 1
replace var = 2 if variety == "Wisconsin No. 38" (6 real changes made)
replace var = 3 if variety == "No. 457" (6 real changes made)
replace var = 4 if variety == "Glabron" (6 real changes made)
replace var = 5 if variety == "Peatland" (6 real changes made)
replace var = 6 if variety == "Velvet" (6 real changes made)
replace var = 7 if variety == "No. 475" (6 real changes made)
replace var = 8 if variety == "Manchuria" (6 real changes made)
replace var = 9 if variety == "No. 462" (6 real changes made)
replace var = 10 if variety == "Svansota" (6 real changes made)
label define varlbl 1 "Trebi" 2 "Wisconsin No. 38" 3 "No. 457" /// 4 "Glabron" 5 "Peatland" 6 "Velvet" 7 "No. 475" 8 "Manchuria" /// 9 "No. 462" 10 "Svansota"
label values var varlbl
graph dot yield1931 yield1932, /// over(var, gap(20) label(labsize(vsmall))) /// by(site2, compact cols(2) b1title("Barely Yield (bushels/acre)") /// note("")) ylabel(15(15)60) marker(1,msymbol(+)) marker(2,msymbol(Oh)) /// legend(label(1 1931) label(2 1932))
page 7, Figure 1.2
use "c:\vizdata\singer.dta", clear
histogram height, by(voice, col(2) note("")) bin(20) /// ylabel(0 .20 .40, nogrid) xtitle("Height (inches)") /// ytitle("Percent of Total") xlabel(60(5)75) xsize(1) ysize(1.2) /// ylabel( , angle(0))
To have the panels ordered as they are in the text, you need to create a new variable and label its values.
gen voice1 = 1
replace voice1 = 2 if voice_part == "Soprano 1" (36 real changes made)
replace voice1 = 3 if voice_part == "Alto 2" (27 real changes made)
replace voice1 = 4 if voice_part == "Alto 1" (35 real changes made)
replace voice1 = 5 if voice_part == "Tenor 2" (21 real changes made)
replace voice1 = 6 if voice_part == "Tenor 1" (21 real changes made)
replace voice1 = 7 if voice_part == "Bass 2" (26 real changes made)
replace voice1 = 8 if voice_part == "Bass 1" (39 real changes made)
label define voicelbl 1 "Soprano 2" 2 "Soprano 1" 3 "Alto 2" /// 4 "Alto 1" 5 "Tenor 2" 6 "Tenor 1" 7 "Bass 2" 8 "Bass 1"
label values voice1 voicelbl
histogram height, by(voice1, col(2) note("")) bin(20) /// ylabel(0 .20 .40, nogrid) xtitle("Height (inches)") /// ytitle("Percent of Total") xlabel(60(5)75) xsize(1) ysize(1.2) /// ylabel( , angle(0))
page 9, Figure 1.3
use "c:\vizdata\polarization.dta", clear
scatter babinet concentration, xlabel(0(40)120) /// ylabel(14(4)26, nogrid angle(0)) xtitle("Concentration (ug/m3)") /// ytitle("Babinet Point (degrees)") msymbol(Oh)
page 11, Figure 1.4
use "c:\vizdata\ethanol.dta", clear
scatter nox e, ylabel(1(1)4, nogrid angle(0)) xlabel(.6(.2)1.2) /// ylabel(1(1)4) ymtick(.5(1)3.5) msymbol(Oh) xmtick(.7(.2)1.2) /// ytitle("NO x (ug/J)")
page 11, Figure 1.5
scatter nox c, ylabel(1(1)4, nogrid angle(0)) xlabel(10(4)18) ///
ylabel(1(1)4) ymtick(.5(1)3.5) msymbol(Oh) xmtick(8(2)12) xtitle("C") ///
ytitle("NO x (ug/J)")
page 13, Figure 1.6
use "c:\vizdata\barley.dta", clear
graph dot yield, over(variety, label( labsize(small))) /// by(site year, title("Barely Yield (bushels/acre)") /// note("") cols(2)) ylabel(15(15)60) ysize(4) xsize(2) ymtick(20(5)65)
You can use the code below to have the graph ordered as it is in the text.
use "c:\vizdata\barley.dta", clear
gen site2 = 1 if site=="Waseca" & year == 1932 (110 missing values generated)
replace site2 = 2 if site=="Waseca" & year == 1931 (10 real changes made)
replace site2 = 3 if site=="Crookston" & year == 1932 (10 real changes made)
replace site2 = 4 if site=="Crookston" & year == 1931 (10 real changes made)
replace site2 = 5 if site=="Morris" & year == 1932 (10 real changes made)
replace site2 = 6 if site=="Morris" & year == 1931 (10 real changes made)
replace site2 = 7 if site=="University Farm" & year == 1932 (10 real changes made)
replace site2 = 8 if site=="University Farm" & year == 1931 (10 real changes made)
replace site2 = 9 if site=="Duluth" & year == 1932 (10 real changes made)
replace site2 = 10 if site=="Duluth" & year == 1931 (10 real changes made)
replace site2 = 11 if site=="Grand Rapids" & year == 1932 (10 real changes made)
replace site2 = 12 if site=="Grand Rapids" & year == 1931 (10 real changes made)
label define sitelbl 1 "Waseca" 2 "Waseca" 3 "Crookston" 4 "Crookston" /// 5 "Morris" 6 "Morris" 7 "University Farm" 8 "University Farm" /// 9 "Duluth" 10 "Duluth" 11 "Grand Rapids" 12 "Grand Rapids"
label values site2 sitelbl
gen var = 1
replace var = 2 if variety == "Wisconsin No. 38" (12 real changes made)
replace var = 3 if variety == "No. 457" (12 real changes made)
replace var = 4 if variety == "Glabron" (12 real changes made)
replace var = 5 if variety == "Peatland" (12 real changes made)
replace var = 6 if variety == "Velvet" (12 real changes made)
replace var = 7 if variety == "No. 475" (12 real changes made)
replace var = 8 if variety == "Manchuria" (12 real changes made)
replace var = 9 if variety == "No. 462" (12 real changes made)
replace var = 10 if variety == "Svansota" (12 real changes made)
label define varlbl 1 "Trebi" 2 "Wisconsin No. 38" 3 "No. 457" /// 4 "Glabron" 5 "Peatland" 6 "Velvet" 7 "No. 475" 8 "Manchuria" /// 9 "No. 462" 10 "Svansota"
label values var varlbl
graph dot yield, over(var, label(labsize(small))) /// by(site2 year, title("Barely Yield (bushels/acre)") note("") cols(2)) /// ymtick(20(5)65) ylabel(15(15)60) ysize(4) xsize(2)
UCLA Researchers are invited to our Statistical Consulting Services
We recommend others to our list of Other Resources for Statistical Computing Help
These pages are Copyrighted (c) by UCLA Academic Technology Services
