다중 회귀 모델을 통해 두 예측이 크게 다른지 테스트합니다.저는 JMP 16에서 다음을 기반으로 뇌졸중 환자의 임상 결과(" 퇴원 mRS "라는 수치 척도 형식)를 예측하는 다중 선형 회귀 모델을 작업하고 있습니다.
- 1개의 연속 예측변수( 입원 시 NIHSS , 병원 입원 시 뇌졸중 심각도를 측정함)
- 1개의 범주형 예측 변수(뇌졸중 환자가 받을 수 있는 치료인 EVT 타이밍 )는 다음 3가지 값 중 하나를 가질 수 있습니다.
- 1 = 병원에 도착 하자마자 EVT(치료)가 투여되었습니다.
- 2 = 지연 후 EVT가 시행 되었습니다.
- 0 = EVT가 관리되지 않았습니다 .
- 상호 작용 기간( 입학 시 NIHSS * EVT 타이밍 )
현재 JMP에서 모델은 다음과 같습니다.
내 질문/문제는:
입원 시 NIHSS 의 특정 값에 대해 환자가 즉시 EVT를 받는 경우(예: EVT 타이밍 = 1 )와 예측된 평균 결과(이 모델을 기반으로 함)가 크게 다른지 (p 값 사용) 여부를 계산하고 싶습니다. . EVT 없음 (예: EVT 타이밍 = 0 ) - 즉, 범주형 예측 변수의 서로 다른 2개 값 사이입니다.
JMP에서 이러한 분석을 수행할 수 있는 방법이 있나요?
( 개념적으로 내 목표는 임상적으로 유용한 방식으로 이 모델의 2가지 예측변수를 활용할 수 있는 것입니다. 예를 들어 모델이 어떤 환자(이 경우 뇌졸중의 심각도에 따라)를 결정하는 데 도움이 될 수 있기를 바랍니다. 치료를 받을지 여부에 따라 해당 치료의 결과를 예측합니다.)
당신의 도움을 주셔서 감사합니다!
기본 데이터 분석 및 모델링윈도우Re: 다중 회귀 모델을 통해 두 예측이 크게 다른지 테스트합니다.안녕하세요 @MRB3855 님 ,
매우 감사합니다! 나는 그것이 효과가 있다고 생각합니다. 분명히 나는 테스트의 출력을 직접 확인하기 위해 수학을 수행하지 않았지만 (방법을 몰랐습니다) 이제 사용자 정의 테스트의 예측 값이 내 모델의 예측 값과 일치하고 테스트 출력이 의미가 있습니다.
도움을 주셔서 정말 감사합니다!
Re: 다중 회귀 모델을 통해 두 예측이 크게 다른지 테스트합니다.안녕하세요 @broccoli : 죄송합니다. 실수를 했습니다. 시간=24(내 데이터)에서의 예측은 아래를 참조하세요. 상호 작용 기간을 입력하려면 (NIHSS 입학시 - 14.7059). 따라서 입학 = 20에 대해 NIHSS에서 예측하려면 상호 작용을 위해 20-14.7059 = 5.2941을 입력해야 합니다.
그리고 나열되지 않은 배치(4개의 배치가 있음)의 경우 다음과 같이 매개변수화되어야 합니다.
이상해 보이는데, 나도 알아요; 이는 JMP가 모델을 매개변수화하는 방식 때문입니다.
https://www.jmp.com/support/help/en/18.0/index.shtml#page/jmp/nominal-factors.shtml#
그리고 테스트할 수 없는 대비는 모든 대비(열)에 대한 결합 테스트를 제공하는 출력 하단의 제곱합을 나타냅니다. 이 경우 대비가 직교하지 않기 때문에 의미가 없습니다(비공식적으로는 독립적이라고 생각할 수 있음). 그러나 세 가지 대비 각각에 대한 테스트는 개별적으로 유효하고 정확합니다.
직교성에 대한 자세한 내용은 여기를 참조하세요. https://stats.stackexchange.com/questions/171324/what-is-the-relationship-between-orthogonal-correlation-and-independentence
Re: 다중 회귀 모델을 통해 두 예측이 크게 다른지 테스트합니다.안녕하세요 @broccoli 절편을 잊어버렸습니다(처음 두 열에 1을 입력하세요).
그리고 입학 시 NIHSS가 무엇인지가 중요하지 않다고 확신하시나요?
편집: 흠, 공변량과 관련된 상호 작용 용어가 있을 때 예상대로 작동하는지 잘 모르겠습니다.
계속 지켜봐 주시기 바랍니다.
Re: 다중 회귀 모델을 통해 두 예측이 크게 다른지 테스트합니다.안녕하세요 @MRB3855 님 , 이 목적을 위해 사용자 정의 테스트를 설정하려고 시도했지만 올바르게 수행할 수 없는 것 같습니다.
예를 들어 아래 그림에서는 EVT Timing = 0 (열 1)과 EVT Timing = 1 (열 2)을 비교하려고 하며, NIHSS는 입학 허가 = 0 입니다.
테스트 출력의 처음 두 열( 1.467 및 2.100 )의 " 값 " 행에 나열된 숫자는 내 모델이 해당 두 조건에서 예측하는 값(각각 0.0563 및 3.568 - 아래 예측 프로파일러 참조)과 동일하지 않습니다. .
마찬가지로, 입학 = 20 대신 NIHSS를 사용하여 이 사용자 정의 테스트를 반복하려고 시도하면(아래) 출력된 값 의 처음 두 열( 2.329 및 2.962 )이 다시 내 모델의 예측(각각 5.105 및 5.089 )과 일치하지 않습니다.
내가 여기서 무엇을 놓치고 있는 걸까요? 분명히 그것은 중요한 것입니다. 테스트 " 결과 " - 내가 올바르게 이해하고 있다면 세 번째 출력 열에 나열되어야 하기 때문입니다 ( 값 예측의 차이, Std Error , t Ratio , p-value 및 SS ). - 입학 시 NIHSS를 설정한 것과 상관없이 똑같은 것 같습니다(절대 그렇지 않습니다).
불편을 끼쳐드려 죄송하고 다시 한번 감사드립니다!
ps. 참고로 " 경고: 테스트할 수 없는 대비 "는 무엇을 의미합니까?
Re: 다중 회귀 모델을 통해 두 예측이 크게 다른지 테스트합니다.안녕하세요 @broccoli : 네, 가능합니다. 출력 상단의 빨간색 삼각형을 클릭하고 Estimates/Custom Test를 선택합니다(아래 그림 참조). 여기에 예가 있습니다.
사용자 정의 테스트의 예(jmp.com)
하지만 까다로울 수 있습니다. 살펴보고 도움이 필요하면 다시 오세요. 두 개의 열을 추가해야 합니다. 그런 다음 관심 있는 EVT 타이밍(및 상호 작용)에 1을 입력하고 처음 두 열에 관심 있는 입학 시 NIHSS 값을 선택합니다. 그런 다음 세 번째 열은 처음 두 행의 차이입니다.
다음은 시간=3.14(첫 번째 열 101.2)의 배치 1_11과 시간=3.14(두 번째 열 101.4)의 배치 3_13에 대한 예측을 비교하는 예입니다. 세 번째 열은 예측의 차이입니다(-0.202, p-값 0.6758).
