교수님 안녕하십니까? 훌륭한 책과 자료들을 통해 많은 도움을 주셔서 감사합니다.
mytable 을 통해 작업 시간을 획기적으로 줄일 수 있는 것을 보고 감탄을 했습니다. 이것만으로도 R을 쓸 이유가 충분할 것 같습니다.
제가 가지고 있는 데이터에 두 그룹의 baseline characteristics 를 비교하기 위한 table 을 만들고자 교수님이 만들어주신 mytable 을 사용해 보았는데요,
문제는 변수가 factor 인 것은 결측치까지 포함하여 결측치의 수와 비율, 나머지 데이터의 평균, 표준편차를 잘 구해주지만,
변수가 numeric 일 때에는 결측치가 아예 나타나지 않는다는 것입니다.
제가 의도하는 것은 STROBE recommendation 에도 나와있듯이 각 변수들의 결측치의 수를 table 에 표시하면서 결측치를 제외한 나머지 값들의 평균과 표준편차 (또는 중위수와 사분위수) 를 table 로 만들고 싶은 건데요,
혹시 방법이 있을까요? 결측치를 NULL 같은 문자로 바꾸면 numeric 으로 되어 있는 다른 값들도 같이 문자화되어 분석이 안 되는 문제가 있습니다.
또 하나 여쭙고 싶은 점은, 앞서 언급했던 STROBE recommendation 에서는 표준편차나 significance tests 를 descriptive tables 에 사용하지 말 것을 권고하고 있는데 그럼 mytable 에서와 같은 다른 그룹 간의 특성을 비교하고자 하는 통계분석은 후향적연구에서는 배제되어야 하는 건가요? 혹시 이에 대해 교수님께서는 어떻게 생각하시는지 궁금합니다.
많은 도움 주심에 미리 감사의 인사를 드립니다.
Comment 1
-
cardiomoon
2017.10.31 19:55
No. | Subject | Author | Date |
---|---|---|---|
1034 | 로지스틱 회귀분석 질문드립니다. [2] | 윤이아빠 | 2018.11.06 |
1033 | forest plot 그리는 법 문의 드립니다. [1] | 강민정 | 2019.02.24 |
1032 | Rstudio 서버 미러링 방법 궁금해요 [1] | mjcho | 2018.01.11 |
1031 | 웹에서 클릭만으로 하는 R통계분석에서 궁금한점 [2] | osr33260216 | 2016.02.18 |
1030 | 웹R에서 주석 위치 바꾸는 법 [2] | alse | 2016.09.17 |
1029 | r studio server expired..여서 접속이 안됩니다 [2] | 스튜던트 | 2018.04.22 |
1028 | 보고서/플롯 옵션 관련 문의입니다. [1] | wlsrn | 2016.12.12 |
1027 | 메타분석에서 forest plot 크기 [1] | SumiHan | 2016.05.04 |
1026 | 불충분한 데이터의 생존분석에 관하여 질문드립니다. [2] | fermion503 | 2018.05.14 |
1025 | mytable 문의드립니다. [1] | seyi0918 | 2022.11.30 |
1024 | ORplot x축 표시값 변경문의 [3] | 영진 | 2017.04.28 |
» | mytable 사용 시 결측치의 표시 [1] | twogama2 | 2017.10.30 |
1022 | ggplot2에 관하여 | 쑤니 | 2018.12.30 |
1021 | 변수 변환 결과 해석 관련... [3] | 하늘산책 | 2015.12.11 |
1020 | 세군의 비율에 대한 비교와 사후검정에 대해 질문드립니다. | 통계장님 | 2018.07.11 |
1019 | 상호작용효과를 폴리곤 그래프로 만들 때 오류 [3] | 반디12 | 2016.01.12 |
1018 | PSM caliper 문의드립니다. | elsa | 2019.03.28 |
1017 | 로그변환 후 회귀분석 시 결과물 [1] | 반디12 | 2016.07.19 |
1016 | 변수선택 관련 질문 드립니다. [1] | chamchi | 2019.02.28 |
1015 | stepwise logistic 회귀에서 결측치 문제 질문드립니다. [2] | 매직하마 | 2019.02.10 |
일단 첫번째 질문에 대한 답은 mytable함수를 쓸 때 show.all=TRUE옵션을 주면 해결됩니다. 다음을 참조하세요
> mytable(sex~.,data=acs,show.all=TRUE)
Descriptive Statistics by 'sex'
_____________________________________________________________________________________________________
Female Male p sig p1 p2 p3 class ptest N
(N=287) (N=570)
-----------------------------------------------------------------------------------------------------
age 68.7 ± 10.7 60.6 ± 11.2 0.000 *** 0.001 0.000 0.000 continuous normal 857
cardiogenicShock 0.136 0.136 2.000 categorical chisq.test 857
- No 275 (95.8%) 530 (93.0%)
- Yes 12 ( 4.2%) 40 ( 7.0%)
entry 0.035 ** 0.035 2.000 categorical chisq.test 857
- Femoral 119 (41.5%) 193 (33.9%)
- Radial 168 (58.5%) 377 (66.1%)
Dx 0.012 ** 0.012 3.000 categorical chisq.test 857
- NSTEMI 50 (17.4%) 103 (18.1%)
- STEMI 84 (29.3%) 220 (38.6%)
- Unstable Angina 153 (53.3%) 247 (43.3%)
EF 56.3 ± 10.1 55.6 ± 9.4 0.387 0.000 0.387 0.053 continuous normal 723
height 153.8 ± 6.2 167.9 ± 6.1 0.000 *** 0.000 0.000 0.000 continuous normal 764
weight 57.2 ± 9.3 68.7 ± 10.3 0.000 *** 0.000 0.000 0.000 continuous normal 766
BMI 24.2 ± 3.6 24.3 ± 3.2 0.611 0.000 0.611 0.471 continuous normal 764
obesity 0.580 0.580 2.000 categorical chisq.test 857
- No 194 (67.6%) 373 (65.4%)
- Yes 93 (32.4%) 197 (34.6%)
TC 188.9 ± 51.1 183.3 ± 45.9 0.124 0.000 0.124 0.145 continuous normal 834
LDLC 117.8 ± 41.2 116.0 ± 41.1 0.561 0.000 0.561 0.465 continuous normal 833
HDLC 39.0 ± 11.5 37.8 ± 10.9 0.145 0.000 0.145 0.011 continuous normal 834
TG 119.9 ± 76.2 127.9 ± 97.3 0.195 0.000 0.195 0.804 continuous normal 842
DM 0.077 * 0.077 2.000 categorical chisq.test 857
- No 173 (60.3%) 380 (66.7%)
- Yes 114 (39.7%) 190 (33.3%)
HBP 0.000 *** 0.000 2.000 categorical chisq.test 857
- No 83 (28.9%) 273 (47.9%)
- Yes 204 (71.1%) 297 (52.1%)
smoking 0.000 *** 0.000 3.000 categorical chisq.test 857
- Ex-smoker 49 (17.1%) 155 (27.2%)
- Never 209 (72.8%) 123 (21.6%)
- Smoker 29 (10.1%) 292 (51.2%)
-----------------------------------------------------------------------------------------------------