DATA ALL; SET STATS;
RETAIN MS_ERR DF_ERR MS_LAB DF_LAB S_R S_RR ;
IF _TYPE_ = ‘SS1’ THEN DELETE;
IF _SOURCE_ = ‘ERROR’ THEN DO;
MS_ERR = SS/DF;
DF_ERR = DF;
END;
IF _SOURCE_ = ‘LAB’ THEN DO;
MS_LAB = SS/DF;
DF_LAB = DF;
IF MS_LAB < = MS_ERR THEN SIGMA2_L = 0;
ELSE;
SIGMA2_L = (MS_LAB - MS_ERR)/&REPS;
S_R = SQRT(MS_ERR);
S_RR = SQRT(MS_ERR + SIGMA2_L);
OUTPUT;
END;
RUN;
PROC MEANS NOPRINT DATA=FSIM;
BY I;
VAR X ;
OUTPUT OUT=A N=N MEAN= XBAR;
RUN;
DATA AB; SET A;
N_LABS = &N_LABS;
REPS = &REPS;
DROP _TYPE_ _FREQ_;
RUN;
DATA VV; MERGE AB ALL;
RSD_R = ROUND(100*(S_R/XBAR),.01);
RSD_RR = ROUND(100*(S_RR/XBAR),.01);
THETA1 = S_R/S_RR;
RUN;
PROC SORT DATA=VV;
BY RSD_RR;
RUN;
PROC FREQ DATA = VV;
TABLES RSD_RR;
RUN;
DATA D; SET VV;
LOG10_MU = LOG10(&C);
POP_THETA = &THETA;
POP_RSD = &XI_R.;
KEEP N_LABS REPS PCTILE POP_RSD LOG10_MU RSD_RR POP_THETA;
DO PCTILE = .99, .95, .90, .80, .70, .60, .50, .40, .30, .20, .10, .05, .01;
J=CEIL(PCTILE*&TEST);
SET VV POINT=J;
RSD_RR=RSD_RR;
OUTPUT D;
END;STOP;
RUN;
PROC PRINT DATA=D NOOBS;
VAR POP_RSD LOG10_MU POP_THETA N_LABS REPS RSD_RR PCTILE;
RUN;
M
C
C
LURE
& L
EE
: J
OURNAL OF
AOAC I
NTERNATIONAL
V
OL
. 89, N
O
. 3, 2006
803
Recommended to OMB by Committee on Statistics: 07-17-2013
Reviewed and approved by OMB: 07-18-2013
34
65
