/* Do file for RIA, radioimmunoassay analysis */

summarize b100  /* b100 is the 100 pct binding cpm in triplicate (or more) */

generate Bo=_result(3)  /* Compute average 100 pct binding */

label variable Bo"Average 100pct binding;triplicate"

summarize NSB   /* Nonspecific binding in triplicate */

generate N=_result(3) /* Average nonspecific binding*/

label variable N"Average nonspecific binding"

generate B_To=(Bo-N)/(B-N) /* Corrected fraction of total bound */

label variable B_To"Fraction of Total bound"

generate Y =(cpm-N)/(Bo-N) /* Ratio of cpm to 100 pct bound cpm (corrected) */

label variable Y"Fraction of standard bound"

generate logitY= log(Y/(1-Y)) /* Logit transform ln(Y/(1-Y)) */

label variable logitY"Logit transform ln(Y/(1-Y))"

generate logconc= log(stds_pg)/log(10) /* Log (base 10) of concentration */

label variable logconc"Log(base 10) of conc.,pg/ml" 

graph logitY logconc, c(m) border /* View the graph */

quietly regress logitY logconc /* Perform linear regression */

predict hat /* Values for the regression line */

label variable hat"Predicted least squares fit"

graph logitY  hat logconc, c(.l) s(Ti) /* View fit of regression line */

noisily regress /* View the regression table and R squared */