[SAS] 패널 분석 매크로(2) (이원 임의효과 고정효과 하우즈만 검정)

안녕하세요,

지난 포스팅에서 사회과학 분야의 학위논문 작성에 많이 활용되는 패널 회귀분석을 위한 SAS 매크로 코드를 소개하였습니다. 아래의 포스팅은 개체의 효과만 고려하는 일원 효과 모형이었습니다. 

[SAS] 패널 분석 매크로(1) (일원 임의효과 고정효과 하우즈만 검정)

이번 포스팅에서는 SAS ODA를 활용하여 클릭만으로 이원 효과 패널회귀분석을 하는 방법을 살펴보겠습니다. 이원 효과 모형은 개체효과와 시간의 효과를 함께 고려하는 방법입니다. 아래의 코드를 우선 실행하고.

/* 이원 */;
%MACRO PANELTWO(DATA=, CROSS=, TIME=, DV=, IV=, CLASS=, CV=, MV=, INTER=);
PROC SORT DATA=&DATA;
	BY &CROSS &TIME;
RUN;

TITLE "확률효과모형";
PROC PANEL DATA=&DATA;
	CLASS &CLASS;
	ID &CROSS &TIME;
	MODEL &DV=&IV &CLASS &CV &MV &INTER/RANTWO;
RUN;

TITLE "고정효과모형";
PROC PANEL DATA=&DATA;
	CLASS &CLASS;
	ID &CROSS &TIME;
	MODEL &DV=&IV &CLASS &CV &MV &INTER/FIXTWO;
RUN;

TITLE "Pooled";
PROC PANEL DATA=&DATA;
	CLASS &CLASS;
	ID &CROSS &TIME;
	MODEL &DV=&IV &CLASS &CV &MV &INTER/POOLED;
RUN;
%MEND;

이후 아래의 코드를 필요에 맞게 변형하면 됩니다. 조절변수(MV)와 상호작용 변수(INTER)를 입력하면 조절효과를 고려한 패널분석이 가능합니다.

 

%PANELTWO(DATA=데이터셋 이름, CROSS=횡단면(지역, 기업 등), TIME=시간(연도 등), DV=종속변수, IV=독립변수, CLASS=범주형 변수, CV=통제변수, MV=조절변수, INTER=상호작용항);

 

SAS PROC PANEL 프로시저 예제에서 제공하는 airline 데이터셋을 활용해서 코드를 연습해 볼 수 있습니다. airline데이터셋은 SAS Help Center에서 제공하는 코드를 이용하여 만들 수 있는데 아래의 링크를 참고해 주세요.

SAS Help Center

아래의 코드는 이원효과 모형을 사용하는 매크로를 실행하려고 하는데 ID는 airlineid, 시간은 year인 데이터셋 airlined를 분석하는데 독립변수는 lq, 종속변수는 lq로 하는 모형을 분석하라는 의미입니다.

%PANELTWO(DATA=airline, CROSS=airlineid, TIME=year, DV=lc, IV=lq, CLASS=, CV=, MV=, INTER=);

이후 분석결과는 일원 효과 모형과 동일하게 확률효과모형, 고정효과모형, Pooled OLS 결과가 차례대로 출력되고 하우즈만 검정, F-검정 통계량을 확인할 수 있습니다. 그리고 분석결과에 제목이 함께 출력되기 때문에 결과를 파악하는데 용이합니다. 

 

이상으로 개체효과와 시간효과를 함께 고려하는 패널회귀분석 매크로 코드를 살펴보았습니다.

 

궁금하신 점은 언제든 댓글로 남겨주세요:)

감사합니다!

☆ 제 코드는 정답이 아니며 틀린 부분이 있을 수 있으니 주의하여 활용하시기 바랍니다.