안녕하세요.
지금 회귀분석을 하는데, 상호작용 효과 모델을 사용합니다.
R 코드를 찾긴 했는데 (첨부파일)
제 데이터의 회귀모델을 넣으면 자꾸 아래와 같은 오류가 뜹니다. 해결 방법이 있을까요?
Error in model$coefficients[[interaction]] : subscript out of bounds
참고로 제 데이터의 관측치가 5000개 정도이고 NA 뜨는 건 상당히 많습니다.
대부분의 값은 log로 그 값을 변형했구요.
interaction ㅡ로 쓰는 변수들 또한 많은 NA를 포함합니다. 상호작용 효과에 쓰이는 변수들의 NA를 모두 제거하니 관측치가 700개 가량 남더군요.
대신 누락되는 관측치를 모두 없애니 환경코너에
Large lm (13elements, ..) 이렇게 표시되던 회귀식의 값이 NA를 모두 제거한 뒤에는 List of 12 이렇게 변하였습니다.
왜 그런건가요? 이것도 궁금하군요. 여튼 이래저래해도 마찬가지로 아래처럼 오류는 뜹니다.
도움주셔서 감사합니다.
p.s. 제 데이터도 첨부했습니다.
Comment 3
-
cardiomoon
2016.01.14 17:17
-
반디12
2016.01.14 18:00
우선 저 코드를 모두 돌리고 회귀모델을 만들어
예제와 같이 하였습니다. 코드의 맨 아랫쪽에 예제가 나와 있더군요.
예제에 쓰이는 제공되는 데이터는 잘 되는데
제 데이터를 집어 넣으면 안 되더라구요. 물론 continous 코맨드에 넣었습니다.
그리고 말씀하신 코맨드를 R에서 실행했는데 아래와 같은 오류가 뜹니다.
Error in mode$coefficents : object of type 'closure' is not subsettable
저 코드를 그대로 쓰고 예제와 같이
m6 <- lm(logmigrant ~ logGlobal + logGDP +logGlobal*logGDP, data=dta)
회귀모델을 만든 뒤interaction_plot_continuous(m6, effect="logGlobal",moderator="logGDP",interaction="logGDP:logGlobal",mean=T,title="",xlabel="",ylabel="")이렇게 돌렸습니다.예제에 첨부된 데이터 파일은 잘 되는데 제 것은 안 됩니다.코드를 발견한 사이트도 함께 첨부합니다.http://causalloop.blogspot.de/2013/06/marginal-effect-plots-for-interaction.html -
반디12
2016.01.17 03:12
혹시 안 되는 건가요?
몇 일 동안 답변을 기다렸는데 답변이 없으시길래..
지금 페이퍼를 쓰고 있는데 이 그래프를 넣으려고 했거든요. ㅜ.ㅜ
No. | Subject | Author | Date |
---|---|---|---|
934 | 교수님 오래간 만입니다. [2] | 에구머니나 | 2020.12.26 |
933 | 저도 생존 분석에 관해 질문 드리고자합니다. [3] | CSJU | 2020.12.23 |
932 | ggplot2 범례 관련하여 도움 요청드립니다. | hamdgogo | 2020.12.22 |
931 | 정회원 후원했습니다 교수님 [1] | 스튜던트 | 2020.12.16 |
930 | 생존분석 관련하여 문의드립니다. [2] | 따거81 | 2020.12.15 |
929 | PSM 질문합니다ㅜ | 카트라이더 | 2020.12.14 |
928 | ggplot으로 지도 그리기 관련 질문 | 양세마리 | 2020.12.14 |
927 | Interval censored data [1] | 김민철 | 2020.12.13 |
926 | 교수님 안녕하세요. 질문드립니다. [1] | JJ88 | 2020.12.09 |
925 | PSM의 matching 방법 기술에 대한 질문입니다. [1] | ng2o321 | 2020.12.01 |
924 | logistic 분석 시 오류 문제 상의드립니다. [1] | britnepak | 2020.12.01 |
923 | mytable을 가지고 ztable시, ztable의 caption, align 등 이 작동을 하지 않습니다. | 김찬식 | 2020.11.30 |
922 | RM-ANOVA 문의드립니다. | 김도 | 2020.11.27 |
921 | 질문드립니다. [3] | JJ1518 | 2020.11.27 |
920 | ggplot2 bar plot 관련 문의입니다. [4] | hamdgogo | 2020.11.26 |
919 | 4 그룹 비교 [2] | 벌레들 | 2020.11.25 |
918 | web-r에서 psm을 수행하였는데 레퍼런스는 어떻게 표시해야 하나요? [1] | 스티브 | 2020.11.24 |
917 | PSM 하려고 하는데 missing value가 없는데 자꾸 에러가 납니다. [1] | 스티브 | 2020.11.18 |
916 | 통계 web R로 돌린 후 통계검증법 확인 [1] | 딸기공주 | 2020.11.17 |
915 | 교수님, wbe-r로 메타분석을 돌렸는데 결과가 보이질 않아요. [2] | hanclan | 2020.11.16 |
위의 코드를 이용하여 어떻게 코드를 작성하셨는지요?
str(mode$coefficents)를 한번 해보시고 결과를 보면 알수 있을 것 같습니다.