문제는 한계 계수 (X의 Y에 대한)에 관한 것입니다. 개별 계수를 해석하는 것에 대해서는별로 생각하지 않습니다. 사람들이 유용하게 지적했듯이, 이들은 선형 및 가산 관계가있을 때 효과 크기로만 식별 할 수 있습니다.
그것이 초점이라면 문제에 대해 생각하는 (개념적으로, 실용적이지는 않지만) 가장 간단한 방법은 다음과 같습니다.
어떤 상호 작용 선형 일반 회귀 모형에서 Y에서 X의 한계 효과를 얻으려면, 당신은 할 수 있습니다 만 그것은하지 알려진 추정되기 때문에 매우 충분하지 않습니다 X. 그러나에 계수 봐. 어쨌든, 한계 효과를 위해 실제로 원하는 것은 X 값 범위에 대한 Y에 대한 예측과 불확실성의 척도를 제공하는 일종의 플롯 또는 요약입니다. 일반적으로 예측 평균 Y와 신뢰 구간을 원하지만 X에 대한 Y의 전체 조건부 분포에 대한 예측을 원할 수도 있습니다.이 분포는 모형 계수에 대한 불확실성을 고려하기 때문에 적합 모형의 시그마 추정치보다 넓습니다. .
이와 같은 간단한 모델을위한 다양한 폐쇄 형 솔루션이 있습니다. 현재의 목적을 위해 우리는 그것들을 무시하고 대신에 임의의 복잡한 모델을 다루는 방식으로 시뮬레이션에 의해 한계 효과 그래프를 얻는 방법에 대해 더 일반적으로 생각할 수 있습니다.
X의 변화가 Y의 평균에 미치는 영향을 원하고 다른 모든 변수를 의미있는 값으로 고치면 행복하다고 가정합니다. X의 새로운 각 값에 대해 모형 계수 분포에서 크기 B 표본을 가져옵니다. R에서 그렇게하는 쉬운 방법은 그것이 평균 coef(model)
과 공분산 행렬을 갖는 보통이라고 가정하는 것입니다 vcov(model)
. 각 계수 세트에 대해 새로운 예상 Y를 계산하고 간격으로 로트를 요약합니다. 그런 다음 X의 다음 값으로 이동하십시오.
이 방법은 각 샘플링 단계에서 변수 (또는 그 역수)를 적용하는 경우 변수에 적용된 멋진 변환의 영향을받지 않는 것 같습니다. 따라서 적합 모형에 log (X)가 예측 변수 인 경우 샘플링 된 계수에 곱하기 전에 새 X를 기록합니다. 적합 모형에 종속 변수로 sqrt (Y)가있는 경우 표본으로 예측 된 각 평균을 제곱하여 구간으로 요약합니다.
요컨대, 프로그래밍은 많지만 확률 계산은 적으며 결과적으로 임상 적으로 이해 가능한 한계 효과입니다. 이 '방법'은 때때로 정치 과학 문헌에서 CLARIFY로 언급되지만 매우 일반적입니다.