/* Do file for RIA sample unknowns.  Run RIA.DO first, please. */

generate logitS=log(((Scpm-N)/(Bo-N))/(1-(Scpm-N)/(Bo-N)))

label variable logitS"Logit transform; unknown samples"

generate pg_ml = 10^((logitS - _b[_cons])/_b[logconc])

label variable pg_ml"Calculated pg/ml of unknowns"

generate answer=pg_ml/Vol_ml /* Results for unknowns, dilution corrected */

label variable answer"Results, dilution corrected"

format pg_ml answer %8.2f

noisily list smpl_num pg_ml answer /* View results, sample no. and conc. */