안녕하세요.

지금 회귀분석을 하는데, 상호작용 효과 모델을 사용합니다. 

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. 제 데이터도 첨부했습니다.