메뉴 건너뛰기

웹에서 하는 R 통계

mytable 관련 질문드립니다

2019.02.12 01:21

낡은찻잔 조회 수:18

교수님 저서로 논문 쓰는데, 도움을 받을 기대로 해보는데 이상하게 잘 되지 않아 문의드립니다.


아래와 같은 에러가 뜹니다.


> mytable(ipp~.,rt10)

Error in lm.wfit(x, y, w, offset = offset, singular.ok = singular.ok,  : 

  NA/NaN/Inf in 'x


ipp라는 것으로 두군을 나눠서 보려하는데 이런 에러가 나는 이유는 무엇인지요?


sex나 다른 지표로는 잘 되는데 ipp만 안되는 이유가 무엇인지 궁금합니다. (너무 초보스럽습니다.ㅜㅜ)


참고로 표는 다음과 같습니다. 좋은 책 읽으면서 R에 대한 흥미도 많이 생겼습니다. 고맙습니다. (rt10 엑셀화일을 첨부드립니다.)



> mytable(rt10)


        Descriptive Statistics        

--------------------------------------- 

                       N      Total    

--------------------------------------- 

 Number               30  47.3 ± 34.1 

 Age                  30  54.3 ± 18.4 

 sex                  30               

   - female                  1  (3.3%) 

   - male                  29  (96.7%) 

 ipp                  30               

   - N                     14  (46.7%) 

   - Y                     16  (53.3%) 

 ECOG                 28               

   - 0                       1  (3.6%) 

   - 1                     20  (71.4%) 

   - 2                      7  (25.0%) 

 LN                   30               

   - No                     6  (20.0%) 

   - Yes                   24  (80.0%) 

 Vis                  30               

   - No                    13  (43.3%) 

   - Yes                   17  (56.7%) 

 Nivo                 30               

   - Nivo: 1mg/kg          12  (40.0%) 

   - Nivo: 3mg/kg          18  (60.0%) 

 Cabo                 29               

   - Cabo: 40mg            23  (79.3%) 

   - Cabo: 60mg             6  (20.7%) 

 Ipi                  16               

   - Ipi: 1mg/kg           14  (87.5%) 

   - Ipi: 3mg/kg            2  (12.5%) 

 Baseline             30               

   - Baseline             30  (100.0%) 

 F                    19               

   - Follow up            19  (100.0%) 

 Histology            30               

   - Mucinous Adeno          1  (3.3%) 

   - P-NET                   1  (3.3%) 

   - Penile                  2  (6.7%) 

   - Prostate               5  (16.7%) 

   - Renal Medullary        3  (10.0%) 

   - Testicular              2  (6.7%) 

   - Urothelial            16  (53.3%) 

--------------------------------------- 

> mytable(sex~.,rt10)

Error in var.test.default(x = 85, y = c(6, 7, 9, 10, 14, 15, 16, 18, 21,  : 

  not enough 'x' observations

Error in var.test.default(x = 24, y = c(44, 66, 75, 33, 62, 58, 44, 70,  : 

  not enough 'x' observations


         Descriptive Statistics by 'sex'         

__________________________________________________ 

                       female       male       p  

                        (N=1)      (N=29)   

-------------------------------------------------- 

 Number              85.0 ±   NA 46.0 ± 33.9      

 Age                 24.0 ±   NA 55.3 ± 17.8      

 ipp                                         1.000

   - N                0 ( 0.0%)  14 (48.3%)       

   - Y               1 (100.0%)  15 (51.7%)       

 ECOG                                        0.211

   - 0                0 ( 0.0%)   1 ( 3.7%)       

   - 1                0 ( 0.0%)  20 (74.1%)       

   - 2               1 (100.0%)   6 (22.2%)       

 LN                                          1.000

   - No               0 ( 0.0%)   6 (20.7%)       

   - Yes             1 (100.0%)  23 (79.3%)       

 Vis                                         1.000

   - No               0 ( 0.0%)  13 (44.8%)       

   - Yes             1 (100.0%)  16 (55.2%)       

 Nivo                                        0.836

   - Nivo: 1mg/kg    1 (100.0%)  11 (37.9%)       

   - Nivo: 3mg/kg     0 ( 0.0%)  18 (62.1%)       

 Cabo                                        1.000

   - Cabo: 40mg      1 (100.0%)  22 (78.6%)       

   - Cabo: 60mg       0 ( 0.0%)   6 (21.4%)       

 Ipi                                         0.242

   - Ipi: 1mg/kg      0 ( 0.0%)  14 (93.3%)       

   - Ipi: 3mg/kg     1 (100.0%)   1 ( 6.7%)       

 Baseline                                         

   - Baseline        1 (100.0%)  29 (100.0%)      

 F                                                

   - Follow up       1 (100.0%)  18 (100.0%)      

 Histology                                   0.157

   - Mucinous Adeno   0 ( 0.0%)   1 ( 3.4%)       

   - P-NET            0 ( 0.0%)   1 ( 3.4%)       

   - Penile           0 ( 0.0%)   2 ( 6.9%)       

   - Prostate         0 ( 0.0%)   5 (17.2%)       

   - Renal Medullary 1 (100.0%)   2 ( 6.9%)       

   - Testicular       0 ( 0.0%)   2 ( 6.9%)       

   - Urothelial       0 ( 0.0%)  16 (55.2%)       

-------------------------------------------------- 

> mytable(ipp~.,rt10)

Error in lm.wfit(x, y, w, offset = offset, singular.ok = singular.ok,  : 

  NA/NaN/Inf in 'x'

번호 제목 글쓴이 날짜 조회 수
공지 묻고답하기 게시판을 운영합니다 [1] cardiomoon 2015.04.20 622
420 Survdiff로 log-rank test 시 P-value 소수점 begabung 2019.02.15 3
» mytable 관련 질문드립니다 [3] file 낡은찻잔 2019.02.12 18
418 web-R 기반 분석 후 다운로드 관련 문의점 kato7 2019.02.11 5
417 stepwise logistic 회귀에서 결측치 문제 질문드립니다. 매직하마 2019.02.10 12
416 문의 드립니다. [1] 김날 2019.01.31 20
415 Risk scoring??? file 홍성엽 2019.01.28 16
414 깃허브 적용 bnmy6581 2019.01.24 14
413 짝지은 자료 (PSM 등 매칭을 통한 자료)에서의 COX(생존분석)방법 의학통계 2019.01.23 76
412 통계 초보 질문있습니다(저에게는 매우 어려운 질문입니다.) file 이스라 2019.01.22 94
411 변수를 변경하는데 질문이 있습니다. [2] 최00 2019.01.21 266
410 NRI, IDI... 쑤니 2019.01.14 29
409 안녕하세요. K-M 곡선에서 censored mark에 대해서 문의드립니다. 간간해요 2019.01.14 40
408 로지스틱 회귀분석의 표본 추출에 관하여 문의드립니다 TAOKAKA 2019.01.12 24
407 안녕하세요~ 통계 질문드립니다. [3] 김두영 2019.01.04 59
406 ggplot2에 관하여 쑤니 2018.12.30 110
405 ios bnmy6581 2018.12.28 94
404 Incidence rate 그리는 방법 ? [1] 심심심심 2018.12.25 49
403 Box plot 그래프를 그리고 싶은데요. [2] file 알은어려워 2018.12.25 54
402 knit PDF 및 ztable에서 표 출력에 대하여 질문드립니다. [5] file zepiros15 2018.12.22 35
401 생존분석에서 연속형변수의 cutpoint에 대해 질문드립니다. [2] zepiros15 2018.12.17 48