평행선 모델에서 상대적 효능을 계산하는 방법은 무엇입니까?안녕,
적합 곡선에 평행선 모델 LogY=Sample+Log_Dose를 적용했습니다.
병렬성 테스트와 모델 매개변수는 괜찮습니다. R이나 PLA 소프트웨어에서 얻은 결과와 정확히 같습니다. 하지만 상대적인 효능은 완전히 다릅니다. 계산에 대한 기술적인 세부 정보를 찾을 수 있는 곳이 있을까요? 감사합니다!
고급 통계 모델링기본 데이터 분석 및 모델링Re: 평행선 모델에서 상대적 효능을 계산하는 방법은 무엇입니까?작은 테스트 데이터 표, 사용된 다양한 코드, 결과 개요를 통해 조사가 용이해집니다.
Re: 평행선 모델에서 상대적 효능을 계산하는 방법은 무엇입니까?예측 및 특수 모델링 > 곡선 맞춤 > 모델 맞춤 옵션 > 병렬 테스트에는 일부 정보가 포함되어 있습니다.
Re: 평행선 모델에서 상대적 효능을 계산하는 방법은 무엇입니까?PLA의 예를 사용했습니다. R의 RP는 2.170981로 PLA 소프트웨어에서 얻은 값과 같습니다. 하지만 JMP의 결과는 5.9입니다.
추정된 매개변수가 동일하므로 둘 다 동일한 모델을 사용하고 있는 것으로 생각됩니다. 다만 JMP의 RP가 너무 큽니다.
데이터 <- read_excel("PLA Example_parallel line.xlsx", sheet = "Sheet1")
데이터$Log_Dose = log(데이터$Dose)
데이터$LogY = log(데이터$Y)
fit <- lm(LogY ~ Sample+Log_Dose , data = data) #평행선 모델: 공통 기울기
요약(적합)
계수 <- coef(fit)
db <- 계수["SampleT"] #
베타 <- 계수["Log_Dose"]
LogRP <- db / 베타
RP<-exp(logRP)
JMP에서 특수 모델링->곡선 적합, 다항식->선형 적합, 병렬성 테스트를 통해 모델을 적합했습니다.
Re: 평행선 모델에서 상대적 효능을 계산하는 방법은 무엇입니까?안녕하세요 @액톤 ,
차이가 보이는 이유는 JMP의 평행선 검정에서 추정치가 S의 '기준선'과의 차이를 사용하지 않고 대신 각 항을 표시하기 때문입니다. 스크린샷(표본 그룹에 대한 더미 코딩 사용)과 동일한 값을 얻을 수 있습니다.
- 절편 = 5.71347 (S의 경우)
- SampleT = 0.70423 → T 절편 = 5.71347 + 0.70423 = 6.4177 (JMP 테이블과 일치)
- 샘플U = 0.51262 → U 절편 = 6.2261
- 샘플V = 0.61602 → V 절편 = 6.3295
표준 오차 추정치는 JMP가 공분산 행렬에서 추정치를 가져오는 반면, 스크린샷에서 얻은 추정치는 원시 추정치처럼 보이기 때문에 다릅니다.
감사해요,
벤
Re: 평행선 모델에서 상대적 효능을 계산하는 방법은 무엇입니까?안녕하세요 벤,
JMP의 절편과 기울기 추정값은 괜찮은데 상대 효능은 정확하지 않은 것 같습니다. '평행성 테스트'를 클릭했더니 상대 효능 부분이 나왔습니다. 평행선 모형의 상대 효능인지 확실하지 않습니다. 이 분석에 대한 JMP 도움말 파일을 찾을 수 없습니다.
감사해요,
액톤
