#delimit ;

clear;

capture log close;
log using lab1, replace;

infile logvlch tx codon10 codon71 codon90 logvl0 using "333data.raw", clear;

label variable logvlch "change in log-10 HIV viral load";
label variable tx "treatment assignment";
label variable codon10 "mutation at codon 10";
label variable codon71 "mutation at codon 71";
label variable codon90 "mutation at codon 90";
label variable logvl0 "log-10 HIV viral load at baseline";

label define mutation 0 "absent" 1 "present";
label values codon10 mutation;
label values codon71 mutation;
label values codon90 mutation;

label define tx 0 "SQVhc" 1 "SQVsgc" 2 "IDV";
label values tx tx;

summarize;
sum logvl0, detail;

tabulate tx;
tab codon10;
tab codon71;
tab codon90;

histogram logvlch;
graph box logvlch;
qnorm logvl0;

gen rannor = invnorm(uniform()); qnorm rannor; drop rannor;
gen rannor = invnorm(uniform()); qnorm rannor; drop rannor;
gen rannor = invnorm(uniform()); qnorm rannor; drop rannor;
gen rannor = invnorm(uniform()); qnorm rannor; drop rannor;

histogram logvlch, title(Histogram of Changes in Log-10 HIV Viral Load);

compress;
save lab1, replace;

log close;