교육통계/Mplus

Mplus 5. 자기회귀교차지연모형

개노답이글스 2022. 11. 18. 17:17

1.자기회귀교차지연 모형 명령문

자기회귀교차지연모형 분석을 위해선 좀 많은 명령문이 필요합니다.

1)무조건모형

2)측정동일성

3)경로동일성

4) 오차공분산동일

 

모형에대한 간단한 설명은 여기서 봐주시고..ㅎ

 

https://gaenodapeagles.tistory.com/23

 

교육통계 R랑가몰라 9. 구조방정식(SEM) 5) 자기회귀교차지연모형

잠재성장모형에 이어 ​ 또 다른 종단 분석방법을 살펴볼까요. ​ 자기회귀교차지연모형입니다. ​ 1. 간단한 이론..? ​ ​ 예시 모형에서 보다시피 여러 시점에서 측정된 변수들 간의 관계를

gaenodapeagles.tistory.com

 

 

또 너무 많으면 토나오니까

저렇게 4개만...

역시 종단모형이니까

종단 파일로 합치셔야 되고,

그다음 명령문으로 들어가야 됩니다.

====무조건모형=======================================

data: file = auto.dat;

format = 80F12.0;

variable:

names = HID PID hl1w1 hl2w1 hl3w1 hl4w1 hl5w1 hl6w1 hl7w1 hl8w1 hl9w1 hl10w1

hl11w1 hl12w1 hl13w1 hl14w1 hl15w1 hl16w1

sad1w1 sad2w1 sad3w1 sad4w1 sad5w1 sad6w1 sad7w1 sad8w1 sad9w1 sad10w1

hl1w2 hl2w2 hl3w2 hl4w2 hl5w2 hl6w2 hl7w2 hl8w2 hl9w2 hl10w2

hl11w2 hl12w2 hl13w2 hl14w2 hl15w2 hl16w2

sad1w2 sad2w2 sad3w2 sad4w2 sad5w2 sad6w2 sad7w2 sad8w2 sad9w2 sad10w2

hl1w3 hl2w3 hl3w3 hl4w3 hl5w3 hl6w3 hl7w3 hl8w3 hl9w3 hl10w3

hl11w3 hl12w3 hl13w3 hl14w3 hl15w3 hl16w3

sad1w3 sad2w3 sad3w3 sad4w3 sad5w3 sad6w3 sad7w3 sad8w3 sad9w3 sad10w3;

usevariable = hl1w1 hl2w1 hl3w1 hl4w1 hl5w1 hl6w1 hl7w1 hl8w1 hl9w1 hl10w1

hl11w1 hl12w1 hl13w1 hl14w1 hl15w1 hl16w1

sad1w1 sad2w1 sad3w1 sad4w1 sad5w1 sad6w1 sad7w1 sad8w1 sad9w1 sad10w1

hl1w2 hl2w2 hl3w2 hl4w2 hl5w2 hl6w2 hl7w2 hl8w2 hl9w2 hl10w2

hl11w2 hl12w2 hl13w2 hl14w2 hl15w2 hl16w2

sad1w2 sad2w2 sad3w2 sad4w2 sad5w2 sad6w2 sad7w2 sad8w2 sad9w2 sad10w2

hl1w3 hl2w3 hl3w3 hl4w3 hl5w3 hl6w3 hl7w3 hl8w3 hl9w3 hl10w3

hl11w3 hl12w3 hl13w3 hl14w3 hl15w3 hl16w3

sad1w3 sad2w3 sad3w3 sad4w3 sad5w3 sad6w3 sad7w3 sad8w3 sad9w3 sad10w3;

model:

help1 by hl1w1 hl2w1 hl3w1 hl4w1 hl5w1 hl6w1 hl7w1 hl8w1 hl9w1 hl10w1

hl11w1 hl12w1 hl13w1 hl14w1 hl15w1 hl16w1;

help2 by hl1w2 hl2w2 hl3w2 hl4w2 hl5w2 hl6w2 hl7w2 hl8w2 hl9w2 hl10w2

hl11w2 hl12w2 hl13w2 hl14w2 hl15w2 hl16w2;

help3 by hl1w3 hl2w3 hl3w3 hl4w3 hl5w3 hl6w3 hl7w3 hl8w3 hl9w3 hl10w3

hl11w3 hl12w3 hl13w3 hl14w3 hl15w3 hl16w3;

sad1 by sad1w1 sad2w1 sad3w1 sad4w1 sad5w1 sad6w1 sad7w1 sad8w1 sad9w1 sad10w1;

sad2 by sad1w2 sad2w2 sad3w2 sad4w2 sad5w2 sad6w2 sad7w2 sad8w2 sad9w2 sad10w2;

sad3 by sad1w3 sad2w3 sad3w3 sad4w3 sad5w3 sad6w3 sad7w3 sad8w3 sad9w3 sad10w3;

help2 on help1;

help3 on help2;

sad2 on sad1;

sad3 on sad2;

help2 on sad1;

help3 on sad2;

sad2 on help1;

sad3 on help2;

help1 with sad1;

help2 with sad2;

help3 with sad3;

output:

standardized;

끔찍하죠..?ㅋㅋ

그래도 뜯어보면 간단합니다.

앞에 'by'로 요인 설정 해주고

'on'으로 어떤 변수가 영향을 주는지 설정해준다.

다만, 교차지연이니까

영향을 서로 주고 받게 설정을 해줍니다.

그리고 하나 더 새로 추가된 것은

with 을 써서 상관을 표시해준다는 것을 확인해주면 되겠습니다.

그리고 일단 돌리면,

다른 것과 마찬가지로 모형 적합도 봐주시고.

 

비표준화계수는 Model results에서

 

표준화계수는 standardized model resurlts - STDYX 에서 봐주면 됩니다.

 

다음에 이제 측정동일성 가야겠죠

측정동일성, 각 시점별로 동일하게 측정되었다는 명령어를 설정되어야 겠죠

나머지 동일하고

model:

help1 by hl1w1 hl2w1 hl3w1 hl4w1 hl5w1 hl6w1 hl7w1 hl8w1 hl9w1 hl10w1

hl11w1 hl12w1 hl13w1 hl14w1 hl15w1 hl16w1(1);

help2 by hl1w2 hl2w2 hl3w2 hl4w2 hl5w2 hl6w2 hl7w2 hl8w2 hl9w2 hl10w2

hl11w2 hl12w2 hl13w2 hl14w2 hl15w2 hl16w2(1);

help3 by hl1w3 hl2w3 hl3w3 hl4w3 hl5w3 hl6w3 hl7w3 hl8w3 hl9w3 hl10w3

hl11w3 hl12w3 hl13w3 hl14w3 hl15w3 hl16w3(1);

sad1 by sad1w1 sad2w1 sad3w1 sad4w1 sad5w1 sad6w1 sad7w1 sad8w1 sad9w1 sad10w1(2);

sad2 by sad1w2 sad2w2 sad3w2 sad4w2 sad5w2 sad6w2 sad7w2 sad8w2 sad9w2 sad10w2(2);

sad3 by sad1w3 sad2w3 sad3w3 sad4w3 sad5w3 sad6w3 sad7w3 sad8w3 sad9w3 sad10w3(2);

여기서 다른점 (1), (2) 가 들어왔습니다.

이렇게 괄호 숫자가 들어가면서, 고정을 해주는 의미로 보면 됩니다.

매시점 학업무기력(우울) 동일하게 고정을 해주는 것으로

결과를 보면,

*명령어가 너무 길어서 제대로 안 먹히네요...

그래서 이렇게 돌렸습니다.

model:

help1 by hl1w1-hl16w1(1);

help2 by hl1w2-hl16w2(1);

help3 by hl1w3-hl16w3(1);

sad1 by sad1w1 sad2w1 sad3w1 sad4w1 sad5w1 sad6w1 sad7w1 sad8w1 sad9w1 sad10w1(2);

sad2 by sad1w2 sad2w2 sad3w2 sad4w2 sad5w2 sad6w2 sad7w2 sad8w2 sad9w2 sad10w2(2);

sad3 by sad1w3 sad2w3 sad3w3 sad4w3 sad5w3 sad6w3 sad7w3 sad8w3 sad9w3 sad10w3(2);

 

 

결과를 확인해보면,

 

 

이렇게 모든 시점에서 

 

요인부하량이 동일하게 들어가도록 해주었습니다.

 

다음으로 경로동일성을 보면

마찬가지로 (3), (4), (5), (6)을 써서 고정해줍니다.

model:

help1 by hl1w1-hl16w1(1);

help2 by hl1w2-hl16w2(1);

help3 by hl1w3-hl16w3(1);

sad1 by sad1w1 sad2w1 sad3w1 sad4w1 sad5w1 sad6w1 sad7w1 sad8w1 sad9w1 sad10w1(2);

sad2 by sad1w2 sad2w2 sad3w2 sad4w2 sad5w2 sad6w2 sad7w2 sad8w2 sad9w2 sad10w2(2);

sad3 by sad1w3 sad2w3 sad3w3 sad4w3 sad5w3 sad6w3 sad7w3 sad8w3 sad9w3 sad10w3(2);

help2 on help1(3);

help3 on help2(3);

sad2 on sad1(4);

sad3 on sad2(4);

help2 on sad1(5);

help3 on sad2(5);

sad2 on help1(6);

sad3 on help2(6);

help1 with sad1;

help2 with sad2;

help3 with sad3;

그리고 돌리면,

 

이렇게 계수값들이 다 동일한 숫자로 고정된 것이 확인 됩니다.

마지막으로 오차공분산동일성

help1 with sad1;

help2 with sad2(7);

help3 with sad3(7);

상관에다(1차년도 빼고) 또 동일성 제약을 걸어줍니다.

돌리면~~

 

 

이렇게 계수값들이 고정된 것을 확인할 수 있습니다.

진짜진짜 마지막으로

그림은 보고 가야겠죠..?

 

 

세상 끔찍합니다..ㅋㅋ

 

마찬가지로 필요한 값들 잘 확인하시고

정리하셔서 사용하시면 될 것 같습니다.

이상 끝!

 

 

'교육통계 > Mplus' 카테고리의 다른 글

Mplus 4. 잠재성장모형  (1) 2022.11.18
Mplus 3. 경로모형  (0) 2022.11.18
Mplus 2. 탐색적 요인분석  (0) 2022.11.17
Mplus 1. 확인적 요인분석  (0) 2022.11.16