Skip to menu

mytable 관련 질문드립니다

낡은찻잔 2019.02.12 01:21 Views : 1062

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


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


> 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'

I want to
List
No. Subject Author Date
1073 Bonferroni's adjustment 관련 문의드립니다. [1] kjptrn 2022.03.23
1072 파일 업로드 문의드립니다. [1] ulysses11 2022.03.07
1071 x 축 변수가 어려개인 그래프를 그릴 수 있을까요? [1] 나현둥 2022.02.20
1070 독립변수가 이상한 (?) 경우의 로지스틱 회귀분석 [4] papahong 2022.02.20
1069 survival tree 분석을 하였는데 node가 한 개만 나옵니다 [1] 뒤통수조심 2022.02.10
1068 예후예측 점수 만든후 그것으로 survival 예측하는 데에 문제 겪고있습니다. [1] 권혜미 2022.02.08
1067 기초적인 통계에서 평균, 중간값, 최빈값에 대한 질문입니다. 우엉우왕 2022.02.04
1066 메타분석 plot 출력 문제와 관련해 답변해 주신 내용에 따라 다시 글을 올립니다. [2] secret endostat 2022.01.27
1065 메타분석 후 plot을 다운로드했을 때 일부 plot이 나오지 않습니다. [1] endostat 2022.01.27
1064 에러가 났는데 어떻게 해야 하나요? [1] 세렌 2022.01.21
1063 multipleROC package 질문 [2] 구현아빠 2022.01.13
1062 교수님, 맥에서 R studio 설치해서 엑셀을 불러오는데 어려움이 있습니다. [1] hslee0511 2022.01.12
1061 multiple imputation 에 관한 질문 [2] 열공빡공 2022.01.08
1060 혹시 median, IQR은 확인할수없나요? [1] ulysses11 2022.01.07
1059 Survival outcome 예측하는 모델 비교 [2] mjsong8705 2022.01.06
1058 출간오류에 대해 문의드립니다. file 구오 2021.12.27
1057 데이터 업로드 오류 [1] 보라도리 2021.12.27
1056 안녕하십니까, ROC curve AUC값 비교 (single predictors, regression models)에 대해서 여쭙습니다. [2] kjg 2021.12.25
1055 ggplot2에서 p값구하기문의 [2] MCB 2021.12.14
1054 기초적인 질문이지만 .. 여쭤봅니다. [2] 나현둥 2021.12.13
Write Tag