나는 맞춤형 테스트를 통해 사후 쌍별 기울기 비교를 심사숙고하려고 합니다. 사용자 정의 테스트에 대한 문서를 읽었지만 테스트를 설정하는 데 사용할 계수를 어떻게 결정해야 하는지 여전히 이해하지 못합니다. 저는 세 가지 치료법[PT CS OF]을 가지고 있으며 이러한 치료법과 연속 변수[기간] 간의 상호 작용에 대한 쌍별 테스트를 수행하고 싶습니다. 사용자 정의 테스트를 올바르게 이해하면 각 테스트는 자체 열을 채울 것입니다.
테스트하고 싶은 사항:
PT*기간-CS*기간=0;
PT*기간-OF*기간=0;
CS*기간-OF*기간=0;
아래에는 테스트 계수와 결과를 보여주는 화면 캡처가 첨부되어 있습니다. 내가 뭔가 잘못하고 있다는 뜻이라고 생각되는 경고가 표시됩니다.
도움을 주셔서 미리 감사드립니다.
자유도 문제인 것 같은데 MANOVA 절차의 테스트에 익숙하지 않습니다. 아마도 다른 회원이 귀하에게 조언을 해줄 수 있을 것입니다. 도움이 필요하면 JMP 또는 SAS 소프트웨어 설명서를 참조하는 것이 좋습니다. JMP 기술 지원팀( support@jmp.com ) 귀하의 질문과 함께.
알겠습니다. 맞춤 테스트와 MANOVA 질문입니다.
ANOVA에서 상호작용 기울기에 대해 3개 열 쌍별 테스트를 수행하면 모든 테스트 통계와 대비를 만들 수 없다는 경고가 표시되지만 최소한 모든 추정치와 테스트 통계를 얻을 수 있습니다. 나는 MANOVA에서 같은 것을 시도했습니다. 네 가지 반응 변수 모두에 대한 요인과 상호 작용에 대해 모델을 적합시켰습니다. (중요하지 않은 상위 상호 작용 중 일부를 제거했습니다.) 응답 옵션으로 "각 열도 테스트" 및 "식별"을 선택했습니다. 특정 열에 대해 "전체 제목 변경"이라고 말합니다. 분산 분석에서와 동일한 방식으로 사용자 정의 테스트를 시도했지만 결과는 단일 F-통계량이었고 대비를 테스트할 수 없다는 동일한 경고였습니다. 단일 쌍 테스트를 수행하면 예상한 결과를 얻을 수 있습니다. 내가 관심이 있고 Bonferroni가 올바른 각 열 테스트에서 10억 개의 쌍별 테스트를 수행합니까?
나는 당신이 올바른 길을 가고 있다고 생각합니다!
효과 테스트를 보셨나요? 데이터를 분석하고 모델링하는 방법은 여러 가지가 있지만 저는 전체 모델로 시작하여 효과 테스트 및 기타 정보를 진행하여 모델을 축소한 다음 사용자 정의 테스트를 사용하여 중요한 효과 내에서 흥미로운 비교를 수행하는 것을 좋아합니다.
JMP 사이트에는 유용한 자료가 많이 있습니다. 나는 가능한 모든 것에 대해 아는 척하지 않습니다. 그래도 검색 가능합니다. 예를 들어, 나는 t를 찾았습니다. 그의 페이지 학술 사용자를 위한 튜토리얼과 비디오가 포함되어 있습니다(누구나 사용할 수 있음).
명확히 하고 싶어서 단일 모델에 모든 독립 요인과 중요한 상호 작용 항을 포함합니다. 동일한 모델의 모든 종속 요인을 수집하는 MANOVA를 실행했는데 전체 모델이 중요합니다. 그 외에도 MANOVA를 개별 표준 최소 제곱 분석으로 다시 분해하여 종속 요인의 동작을 분석하라는 지시를 받았습니다. MANOVA를 사용하여 작업을 해석하는 방법을 설명하는 비디오 또는 슬라이드 데크에 대한 링크가 있으면 감사하겠습니다.
나는 다중 비교의 문제가 하나의 연구와 단일 데이터 세트 내에 있다는 것을 이해합니다. 연구를 반복하면 조정 시 첫 번째 연구에서 이루어진 비교를 계산할 필요가 없습니다. 하지만 동일한 연구에 여러 요인이 있는 것처럼 들리므로 누적된 테스트 수에 따라 조정하면 됩니다.
각 요인에 대해 별도의 분석을 수행하는 대신 단일 모델에 모든 요인을 동시에 포함시키지 않는 이유가 있나요?
후속 질문 하나:
저는 동일한 기본 모델을 다양한 요인에 적용하고 있습니다. 그리고 결국에는 서로 다른 상황에서 동일한 상호 작용을 테스트하게 됩니다.
Bonferroni는 모델 내의 공동 테스트(즉, 3개 수준) 내에서 수정해야 합니까? Bonferroni는 단일 모델 내에서 모든 테스트를 수정해야 합니까, 아니면 모든 모델에 걸쳐 이러한 모든 사용자 정의 테스트를 모아서 하나의 대규모 Bonferroni 수정을 수행해야 합니까?
감사해요
선택한 알파에 Bonferroni 조정을 사용하는 경우 귀하가 요청한 의미를 해석합니다. 예, 이는 다중 비교를 통해 제1종 오류율을 제어하는 한 가지 방법입니다.
공동 테스트는 일부 분석가가 높은 오탐률로부터 보호하기 위해 사용합니다. 많은 테스트를 수행할 수 있습니다. 테스트 횟수에 따라 개별적으로 조정되지 않습니다.
경고는 합동 테스트에 관한 것입니다. 조인트 테스트 수에 사용할 수 있는 자유도가 충분하지 않습니다.
감사합니다. 이것이 제가 찾고 있는 솔루션이며 결과가 의미가 있습니다. 유일한 마지막 문제는 여전히 "경고: 테스트할 수 없는 대비"가 표시된다는 것입니다. 공동 테스트가 테스트 불가능하거나 대조 중 하나가 테스트 가능하지 않다는 뜻인가요? 쌍별 대조를 독립적으로 실행할 때 동일한 응답을 얻지 못합니다.
혼란을 드려 죄송합니다. 모든 혼란이 확실하지 않으므로 처음부터 시작하겠습니다.
당신은 당신이 원한다고 말하면서 시작했습니다 ㅏ 사후 사용자 정의 테스트를 통한 쌍별 기울기 비교 . 따라서 각 그룹의 추세선이 다르게 나타나므로 차이를 테스트하려고 합니다. 해보자.
(참고: 기술 지원 부서의 동료로부터 엄청난 도움을 받았습니다. 내 대답은 그의 설명을 바탕으로 한 것입니다.)
사용자 정의 테스트는 한 단계로 수행할 수 있지만 올바른 계수의 개발을 확인할 수 있도록 두 단계로 나누어 보겠습니다.
1단계: 각 처리에 대한 선의 기울기에 대한 사용자 정의 테스트를 사용합니다. 사용자 정의 테스트의 첫 번째 열은 [0, 0, 0, 1, 1, 0, 0]이어야 합니다. 즉, 1이 있는 유일한 매개변수는 펄스 지속시간(Pulse Duration)과 치료[PT]*펄스 지속시간(Pulse Duration)입니다. 그 이유는 선의 기울기가 공통 기울기인 Pulse Duration에 치료[PT]에 대한 차등 기울기를 더한 값이기 때문입니다. 그 합은 0에 대해 테스트하려는 기울기입니다. 펄스 지속 시간과 치료[CS]*펄스 지속 시간을 추가하므로 두 번째 열은 [0, 0, 0, 1, 0, 1, 0]이어야 합니다. 마지막으로 세 번째 처리에 대한 미분 기울기를 추가하지만 명시적으로 보고되지는 않습니다. 그러나 JMP가 범주형 요인을 매개변수화하는 방식으로 인해 이 값이 다른 두 매개변수의 합에 음수라는 것을 알 수 있습니다. 따라서 세 번째 열은 [0, 0, 0, 1, -1, -1, 0]이어야 합니다.
다시 말하지만, 원하는 실제 사용자 지정 테스트에는 이 단계가 필요하지 않습니다. 실제 테스트에 대한 설명에만 도움이 됩니다. 이 시점에서 우리는 개별 경사면의 사용자 정의 테스트에 대한 올바른 계수를 결정했습니다.
2단계: 개별 기울기를 0이 아닌 한 번에 두 개씩 비교하려고 합니다. 1단계에서 사용된 열을 한 번에 두 개씩 빼기만 하면 됩니다! PT와 CS의 기울기를 비교하기 위해 PT에 대한 계수 열에서 CS에 대한 1단계의 계수 열을 뺍니다. 따라서 [ 0, 0, 0, 1, 1, 0, 0 ] - [ 0, 0, 0, 1, 0, 1, 0 ] = [ 0, 0, 0, 0, 1, -1, 0 ]. PT에서 OF까지 테스트하는 열은 [0, 0, 0, 1, 1, 0, 0 ] - [ 0, 0, 0, 1, -1, -1, 0 ] = [ 0, 0, 0, 0 입니다. , 2, 1, 0]. 마지막으로 CS에서 OF까지 테스트하는 열은 [ 0, 0, 0, 1, 0, 1, 0 ] - [ 0, 0, 0, 1, -1, -1, 0 ] = [ 0, 0, 0 , 0, 1, 2, 0].
이전과 마찬가지로 모수 추정치의 합과 차이를 직접 계산하여 JMP가 테스트에 대해 보고하는 값과 비교하여 항상 계수를 확인할 수 있습니다.
아직도 내가 제대로 하고 있다고 생각하지 않습니다.
첫째, 여전히 대조할 수 없는 경고가 표시됩니다(아래 참조).
둘째, CS*Duration -(-OF*Duration)=0 교차에 대해 두 번째 계수를 -1로 만들어 내가 원하는 것이 아닌 두 값을 추가하는 것으로 보입니다. 결국 -0.0447이 됩니다.
셋째, PT*Duration을 OF*Duration(OF*Duration-PT*Duration=0)과 비교하려는 경우 값은 0.00777이어야 하며 대신 OF*Duration(0.0174926)의 값입니다. 따라서 계수를 0으로 두는 것은 작동하지 않는 것 같습니다.
첫 번째 추정값의 계수를 1로 설정하고 나머지는 0으로 두었다고 생각합니다. 결과를 확인할 수 있습니다. 값은 차이가 0인 null과 비교하는 추정치 간의 차이여야 합니다.
설명해 주셔서 감사합니다. 아래에 다른 창을 첨부했습니다. 내가 올바르게 이해했다면 첫 번째 대비(CS*Duration-OF*Duration=0)를 올바르게 설정한 것입니다. 나는 음수 계수를 사용하여 세 번째 범주형 효과(PT*기간)에 대해 테스트할 수 있다고 생각했습니다. 확장된 추정치에서는 볼 수 있지만 사용자 정의 테스트에 대한 설정에서는 볼 수 없습니다. PT*Duration에 대해 어떻게 테스트합니까?
이러한 사용자 정의 테스트는 다른 곳에서는 대조라고도 합니다.
사양이 누락된 이유는 첫 번째 테스트에서 두 수준 모두에 +1을 입력했기 때문입니다. 계수의 합이 0인지 확인합니다. 그런 다음 "=" 뒤의 값이 귀무가설인지 확인하세요. 일반적으로 유의미한 차이를 검정하는 경우에는 0을 입력하는 것이 적합하지만 일부 상황에서는 그렇지 않을 수도 있습니다.
