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 기능을 이용하여 나란히 그래프를 그리긴 했는데 이 그래프를 한 그래프로 합치는 방법을 혹시 알수 있을 까요? (첨푸파일 참조)
많은 도움을 받고 있습니다. 미리 감사드립니다.
Comment 2
-
cardiomoon
2018.11.07 11:32
-
jh100
2018.11.07 21:07
정말 감사드립니다! 항상 배우고 갑니다
No. | Subject | Author | Date |
---|---|---|---|
1254 | 생존분석 [1] | toby | 2015.04.29 |
1253 | RStudio v0.99 Preview: Code Diagnostics [1] | fermat39 | 2015.04.30 |
1252 | Rstudio 사용법` [3] | 방랑자 | 2015.05.02 |
1251 | need finite 'xlim' values [3] | anthleemd | 2015.05.02 |
1250 | 샤이니 앱에서 SQL 명령어 사용 방법 [2] | 김브링이 | 2015.05.06 |
1249 | 표만들기 관련 [2] | 현승협 | 2015.05.07 |
1248 | MS word문서 저장관련 [1] | 아자 | 2015.05.09 |
1247 | R서버 설치 문의드립니다. [4] | 화모기 | 2015.05.13 |
1246 | 웹에서 하는 R통계 - R코드 [2] | 김영진 | 2015.06.07 |
1245 | Plot inset 만들기 | cardiomoon | 2015.06.09 |
1244 | 오류 질문 - 표작성시 [1] | 미소천사 | 2015.06.15 |
1243 | 오류 질문 - PDF 변환 안되는 file입니다. [1] | 미소천사 | 2015.06.15 |
1242 | print ztable [1] | 김영진 | 2015.06.17 |
1241 | 표 작성시 질문 [1] | simsim | 2015.06.17 |
1240 | 나무 분석 질문입니다. | 미소천사 | 2015.06.19 |
1239 | Interactive ggplot2 [2] | cardiomoon | 2015.06.22 |
1238 | Rstudio server 서버 오류 관련해서 문의드립니다. [2] | 노상래 | 2015.07.22 |
1237 | Rstudio server 무한 로딩 관련해서 문의 드립니다. [6] | 노상래 | 2015.07.24 |
1236 | http://www.r-statistics.com/ 번역소식 [1] | fermat39 | 2015.07.25 |
1235 | 생존분석 변수선택 [3] | nash | 2015.08.11 |
제가 만든 ggiraphExtra 패키지의 ggPredict()함수를 써보셔요
감사합니다.
https://cran.r-project.org/web/packages/ggiraphExtra/vignettes/ggPredict.html