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 |
---|---|---|---|
1274 | logistic regression analysis 가 에러가 뜹니다. [1] | Ikarus86 | 2023.06.02 |
1273 | PSM 문의 [3] | starwars | 2023.04.13 |
1272 | 에러 문의 [2] | 화이팅 | 2023.04.19 |
1271 | moonBook 오류 관련 재문의드립니다. [1] | holyhyun | 2023.04.24 |
1270 | 통계마당 [2] | ejo | 2023.04.27 |
1269 | 연속형 변수인데 범주형으로 인지됩니다. [2] | shararaku | 2023.04.28 |
1268 | 교수님 안녕하세요~~ | 메시12 | 2023.04.29 |
1267 | 메타분석 U3 그래프 관련 문의드립니다 | 화이팅 | 2023.05.02 |
1266 | regression 결과 ppt 다운로드 오류 [2] | JuEuna | 2023.01.26 |
1265 | 반복 분석 자료 분석에 대해서 [2] | ng2o321 | 2023.01.30 |
1264 | 안녕하세요. 정회원인데, 정회원 서버에서 Table 만들기 하면 [2] | brainok | 2023.02.03 |
1263 | Cluster analysis in Web-R [3] | 깨수깡 | 2023.02.07 |
1262 | logistic regression 결과 다운로드 오류 [2] | 남세이짱 | 2023.03.14 |
1261 | mytable을 사용했는데 표가 console에서 보이지 않고 오류도 없는데, 뭐가 문제일까요? [3] | 난둥 | 2023.02.19 |
1260 | 결과표가 다운에러가 납니다. [2] | sophiachoi | 2023.02.21 |
1259 | PSM 문의 | qubic | 2023.02.23 |
1258 | PSM 결과표 문의입니다. | 오공이공 | 2023.03.02 |
1257 | 안녕하세요. [4] | 유령 | 2023.03.12 |
1256 | 보고서 다운로드 에러가 납니다.. | 세렌 | 2023.03.16 |
1255 | group-based trajectory modeling | swpapa | 2023.03.19 |
제가 만든 ggiraphExtra 패키지의 ggPredict()함수를 써보셔요
감사합니다.
https://cran.r-project.org/web/packages/ggiraphExtra/vignettes/ggPredict.html