2012-11-10 6 views
0

고정 효과 모델의 상호 작용 용어를 시각화하는 것을 목표로합니다. 나는 'visreg'와 'ggplot2'패키지를 시도했지만, 불가능한 것 같습니다.r에 패널 모델의 상호 작용 용어를 시각화하는 방법은 무엇입니까?

어떻게하면됩니까?

미리 감사드립니다.

예 데이터 세트 및 패널 모델 :

library(plm) 
data(EmplUK) 

m1 <- plm(emp ~ wage + capital + lag(output,1) + capital:lag(output,1), 
      data = EmplUK, model="within", effect="twoway", index=c("firm","year")) 
+0

"조건부 효과 시각화"란 정확히 무엇을 의미합니까? 데이터에 대한 모델의 적합성을 시각화하고 싶습니까? 1 년 고정 효과? 회사 고정 효과? – RoyalTS

+0

@RoyalTS - 불분명해서 죄송합니다. 상호 작용 용어 (자본 : 지연 (출력, 1))를 시각화하고 싶습니다. – tfr

+0

누군가가 관심이있는 경우 : 선형 모델 (lm)을 사용하고 id 변수에 인형을 추가했습니다. 'visreg'를 사용한 시각화는 문제가 없습니다. – tfr

답변

1

현재 CRAN 버전이 sjPlot 인 경우 sjp.int 함수는 plm 개체를 지원하지 않습니다.

그러나 현재 나는 GitHub-version에서 구현했습니다.

sjp.int(m1)

는 생산 플롯 다음 :

, 위에서 다음 코드를 귀하의 예를 복용

enter image description here

나는 디자인/모델의 모든 종류의 장착 여부를 확실하지 않다 plmsjp.int으로 작동하지만 버그가있는 경우 언제든지 문제를 제기하십시오 (at GitHub).