* Chapter 5: Descriptive Statistics using Stata .

* Import data .
insheet using "C:\temp\rat_brain.dat", tab clear

**************************.
* Descriptive Summary #1 .
**************************.

label define reglabels 1 "BST"
label define reglabels 2 "LS", add
label define reglabels 3 "VDB", add
label define trtlabels 1 "Basal"
label define trtlabels 2 "Carbachol", add

label values treatment trtlabels

label values region reglabels

by treatment, sort : tabstat activate, statistics( count mean sd min max) by(region) columns(variables)

**************************.
* Generate Figure 5.1 .
**************************.

* First, restructure data set to make multiple line plots

reshape wide activate, i( treatment region ) j( animal ) string

rename activateR100797 a100797
rename activateR100997 a100997
rename activateR110597 a110597
rename activateR111097 a111097
rename activateR111397 a111397

label var a100797 "100797"
label var a100997 "100997"
label var a110597 "110597"
label var a111097 "111097"
label var a111397 "111397"

* Generate multiple line plots for each treatment

twoway line a100797 a100997 a110597 a111097 a111397 region, xlabel(1(1)3, valuelabel) lcolor(black black black black black) lpattern(solid dash longdash vshortdash shortdash) ytitle(Activation) by(treatment)