[SAS] 적합결여검정

!'적합결여검정'에서는 proc reg하지 않아도 된다!

data board;
 input x y;
 cards;
40 66.30
40 64.84
40 64.36
45 69.70
45 66.26
45 72.06
50 73.23
50 71.40
50 68.85
55 75.78
55 72.57
55 76.64
60 78.87
60 77.37
60 75.94
65 78.82
65 77.13
65 77.09
;
run;
proc sort data=board out=sorted;  /* x값에 대해 오름차순으로 정렬*/
 by x;
 run;
 proc print data=sorted;
 run;
data sorted;                            
set sorted;
retain tempx lack;
if _n_ =1 then do;
                     lack=1;     /*x 변수가 첫 수준일때 lack에 1을 저장 */
	   tempx=x;
                     end;
else if x=tempx then tempx=x;
                    else do;
                    lack=lack+1;    /* x 변수의 수준이 달라지면 lack값이 1 증가*/
	  tempx=x;
	  end;
run;
proc print data=sorted;
var y x tempx lack;
run;
proc glm data=sorted;
class lack;
model y= x lack;
run;

표를 보고 각각의 값에 어떤 값이 들어가야 하는지 알아야 한다.(SSLF, SSPE 개념도)

F통계량에 대한 p값을 통해 H0 기각 채택 여부 결정 가능하다.