Skip to menu

logistic regression 으로 구한 predicted prababllity 그래프에서 합치는 방법 문의드립니다.

jh100 2018.11.02 17:25 Views : 616

out1 <- glm(DMR_0m~duration, family = binomial, data=tid1)
out0 <- glm(DMR_0m~duration, family = binomial, data=tid0)

par(mfrow=c(1,2))
plot(DMR_0m~duration, data=tid1, type="n", main="acute-onset T1DM")
curve(predict(out1, data.frame(duration=x), type="resp"), add=T)
curve(
  predict(out1, data.frame(duration=x), type="resp", se.fit=T)$fit +
            predict(out1, data.frame(duration=x), type="resp", se.fit=T)$se.fit*1.96, add=T, col="black", lty=2)
curve(
  predict(out1, data.frame(duration=x), type="resp", se.fit=T)$fit -
    predict(out1, data.frame(duration=x), type="resp", se.fit=T)$se.fit*1.96, add=T, col="black", lty=2)

plot(DMR_0m~duration, data=tid0, type="n", main="slowly-progressive T1DM")
curve(predict(out0, data.frame(duration=x), type="resp"), add=T)
curve(
  predict(out0, data.frame(duration=x), type="resp", se.fit=T)$fit +
    predict(out0, data.frame(duration=x), type="resp", se.fit=T)$se.fit*1.96, add=T, col="black", lty=2)
curve(
  predict(out0, data.frame(duration=x), type="resp", se.fit=T)$fit -
    predict(out0, data.frame(duration=x), type="resp", se.fit=T)$se.fit*1.96, add=T, col="black", lty=2)
------------------------------------------------------
위와 같이 종속 변수는 DMR_0m (binomial) 으로 같으나 두개의 subgroup (tid1, tid0) 으로 나누어서 각각 로지스틱 분석 및 predictive probablitiy 검사를 시행하였습니다.
그렇게 하여 mfrow 기능을 이용하여 나란히 그래프를 그리긴 했는데 이 그래프를 한 그래프로 합치는 방법을 혹시 알수 있을 까요? (첨푸파일 참조)

많은 도움을 받고 있습니다. 미리 감사드립니다.
I want to
List
No. Subject Author Date
1213 교수님 지난 번 문의 드렸던 PSM 재문의 드립니다. [3] secret 긍정적인마음가짐 2023.06.27
1212 안녕하세요 ^^ odds ratio 와 관련된 문의드리고자 합니다 ^^ [2] 성열 2023.06.26
1211 mytable() 함수에서 na처리와 pvalue 계산 문의드립니다. [1] hjk34 2023.06.23
1210 교수님 안녕하세요 ^^ web-R table 작성관련 문의드립니다. [2] file mkw7654 2023.06.21
1209 Methods for continuous variables과 P value corrected for multiple comparisons 방법 궁금합니다. [1] 호호셀 2023.06.19
1208 ordinal logistic regression analysis 문의 [1] 버블 2023.06.18
1207 e-mail 변경 문의 [1] secret kimdo 2023.06.16
1206 안녕하세요 number at risk 관련 문의 드립니다 [1] secret cebukim 2023.06.13
1205 Propensity matching score 에 대해 문의드립니다. [1] saga 2023.06.08
1204 콕스비례위험모델의 적용에 관련해서 질문 드립니다. [2] secret kng1032 2023.06.07
1203 기초통계 계산 [2] secret leesav 2023.06.07
1202 정회원 후원하였습니다. [2] secret toast 2023.06.04
1201 logistic regression analysis 가 에러가 뜹니다. [1] secret Ikarus86 2023.06.02
1200 PMS caliper 에 대해 여쭤봅니다. [2] secret pjhaaa 2023.05.30
1199 ggGam에 대하여 [2] leesav 2023.05.22
1198 R 에서 PFS 의 confidence interval 을 구하고 싶을 때 어떻게 하면 될까요? [1] ddeeny 2023.05.21
1197 메타분석 관련 문의드립니다 secret 화이팅 2023.05.12
1196 교수님, 새로운 앱 개발 문의드립니다! [4] swpapa 2023.05.08
1195 교수님, 안녕하세요? 전용 앱 관련 문의드립니다! [1] swpapa 2023.05.07
1194 PSM 시 범주형 변수가 총 10가지 일 때 어떻게 돌려야 하는지요 [1] 긍정적인마음가짐 2023.05.06
Write Tag