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
Write Tag