project 2: BSE
Binh Thang
2022-10-25
library(haven)dd2=read_sav("C:/Users/binht/Dropbox/Hue/data/2022_data_consultant/2 dd2/R/spss_r.sav")class(dd2)## [1] "tbl_df" "tbl" "data.frame"data.label.table <- attr(dd2, "label.table")
missings <- attr(dd2, "missings")names(dd2)## [1] "MaID" "NhomNC" "A1" "A2" "A3" "A4" "A5"
## [8] "A6.1" "A6.1_1" "A6.2" "A6.2_1" "A6.3" "A6.3_1" "A7.1"
## [15] "A7.1_1" "A7.2" "A7.2_1" "A7.3" "A8.1" "A8.2" "A8.3"
## [22] "A8.3_31" "A8.3_32" "A8.4" "B1" "B2" "B3" "B4"
## [29] "B5" "B6" "B7" "B8" "B9" "B10" "B11"
## [36] "B12" "B13" "B14" "B15" "B16" "B17" "B18"
## [43] "B19" "B20" "B21" "B22" "B23" "B24" "B25"
## [50] "B26" "B27" "B28" "B29" "B30" "B31" "B32"
## [57] "B33" "B34" "B35" "B36" "B37" "B38" "B39"
## [64] "B40" "B41" "B42" "TB1_5" "TB6_12" "TB13_18" "TB"
## [71] "TB25_31" "TB32_42" "C0" "C1" "C2" "C3.1" "C3.2"
## [78] "C3.3" "C4" "C5" "C6" "C7" "C8" "C9.1"
## [85] "C9.2" "C9.3" "C9.4" "C10.1" "C10.2" "C9" "TB19_24"library(compareGroups)
library(table1)##
## Attaching package: 'table1'## The following objects are masked from 'package:base':
##
## units, units<-dd2$NhomNC=as.factor(dd2$NhomNC)
dd2$C10.1=as.factor(dd2$C10.1)
dd2$A2=as.factor(dd2$A2)
dd2$A3=as.factor(dd2$A3)
dd2$A4=as.factor(dd2$A4)
dd2$A5=as.factor(dd2$A5)
dd2$A6.1=as.factor(dd2$A6.1)
dd2$A6.2=as.factor(dd2$A6.2)
dd2$A6.3=as.factor(dd2$A6.3)
dd2$A7.1=as.factor(dd2$A7.1)
dd2$A7.2=as.factor(dd2$A7.2)
dd2$A7.3=as.factor(dd2$A7.3)
dd2$A8.1=as.factor(dd2$A8.1)
dd2$A8.2=as.factor(dd2$A8.2)
dd2$A8.3=as.factor(dd2$A8.3)
dd2$A8.4=as.factor(dd2$A8.4)
dd2$C10.2=as.factor(dd2$C10.2)0. Description
a0=table1(NhomNC~C10.1 + A2 + A3 + A4 + A5 + A6.1 + A6.2 + A6.3 + A7.1 + A7.2 + A7.3 + A8.1 + A8.2 + A8.3 + A8.4 + C10.1 + C10.2 + A1 + TB1_5 + TB6_12 + TB13_18 + TB + TB25_31 + TB32_42 + TB19_24, data=dd2)## Warning in table1.formula(NhomNC ~ C10.1 + A2 + A3 + A4 + A5 + A6.1 + A6.2 + :
## Unexpected LHS in formula ignored (table1 expects a 1-sided formula)a0| Overall (N=352) |
|
|---|---|
| NhomNC | |
| 1 | 142 (40.3%) |
| 2 | 104 (29.5%) |
| 3 | 106 (30.1%) |
| C10.1 | |
| 0 | 111 (31.5%) |
| 1 | 95 (27.0%) |
| Missing | 146 (41.5%) |
| A2 | |
| 1 | 2 (0.6%) |
| 2 | 249 (70.7%) |
| 3 | 10 (2.8%) |
| 4 | 91 (25.9%) |
| A3 | |
| 1 | 8 (2.3%) |
| 2 | 330 (93.8%) |
| 3 | 9 (2.6%) |
| 4 | 3 (0.9%) |
| 5 | 2 (0.6%) |
| A4 | |
| 1 | 7 (2.0%) |
| 2 | 96 (27.3%) |
| 3 | 155 (44.0%) |
| 4 | 1 (0.3%) |
| 5 | 63 (17.9%) |
| 6 | 10 (2.8%) |
| 8 | 20 (5.7%) |
| A5 | |
| 1 | 53 (15.1%) |
| 2 | 48 (13.6%) |
| 3 | 112 (31.8%) |
| 4 | 139 (39.5%) |
| A6.1 | |
| 1 | 13 (3.7%) |
| 2 | 339 (96.3%) |
| A6.2 | |
| 1 | 38 (10.8%) |
| 2 | 314 (89.2%) |
| A6.3 | |
| 1 | 11 (3.1%) |
| 2 | 341 (96.9%) |
| A7.1 | |
| 1 | 119 (33.8%) |
| 2 | 233 (66.2%) |
| A7.2 | |
| 1 | 83 (23.6%) |
| 2 | 269 (76.4%) |
| A7.3 | |
| 1 | 229 (65.1%) |
| 2 | 13 (3.7%) |
| 3 | 16 (4.5%) |
| 4 | 10 (2.8%) |
| 5 | 6 (1.7%) |
| 6 | 9 (2.6%) |
| 7 | 69 (19.6%) |
| A8.1 | |
| 1 | 92 (26.1%) |
| 2 | 38 (10.8%) |
| 3 | 222 (63.1%) |
| A8.2 | |
| 1 | 47 (13.4%) |
| 2 | 13 (3.7%) |
| 3 | 292 (83.0%) |
| A8.3 | |
| 1 | 275 (78.1%) |
| 2 | 9 (2.6%) |
| 3 | 67 (19.0%) |
| Missing | 1 (0.3%) |
| A8.4 | |
| 1 | 125 (35.5%) |
| 2 | 58 (16.5%) |
| 3 | 169 (48.0%) |
| C10.2 | |
| 0 | 139 (39.5%) |
| 1 | 65 (18.5%) |
| Missing | 148 (42.0%) |
| Tuoi | |
| Mean (SD) | 45.9 (9.77) |
| Median [Min, Max] | 47.5 [2.00, 73.0] |
| Tong diem nhan thuc nhay cam | |
| Mean (SD) | 11.2 (5.42) |
| Median [Min, Max] | 10.0 [5.00, 25.0] |
| Tong diem muc do nghiem trong | |
| Mean (SD) | 23.0 (7.00) |
| Median [Min, Max] | 24.0 [7.00, 35.0] |
| Tong diem loi ich TKV | |
| Mean (SD) | 21.6 (8.13) |
| Median [Min, Max] | 24.0 [6.00, 30.0] |
| Missing | 1 (0.3%) |
| Tong diem rao can | |
| Mean (SD) | 12.8 (6.33) |
| Median [Min, Max] | 12.0 [6.00, 30.0] |
| Tong diem dong luc TKV | |
| Mean (SD) | 27.5 (6.50) |
| Median [Min, Max] | 28.0 [7.00, 35.0] |
| Tong diem hieu qua ban than | |
| Mean (SD) | 26.2 (13.6) |
| Median [Min, Max] | 22.0 [11.0, 55.0] |
| Tong diem raocan TKV | |
| Mean (SD) | 12.9 (6.60) |
| Median [Min, Max] | 12.0 [6.00, 48.0] |
1. Regular BSE & Factors
t1=compareGroups(C10.1~NhomNC + A2 + A3 + A4 + A5 + A6.1 + A6.2 + A6.3 + A7.1 + A7.2 + A7.3 + A8.1 + A8.2 + A8.3 + A8.4 + C0+ C1+ C2 + C3.1 + C3.2 + C3.3 + C4 + C5 + C6 + C7 +C8 +C9.1+C9.2+C9.3+C9.4+C10.1+C10.2+C9+C10.1 + C10.2 + A1 + TB1_5 + TB6_12 + TB13_18 + TB + TB25_31 + TB32_42 + TB19_24 , data=dd2)## Warning in compare.i(X[, i], y = y, selec.i = selec[i], method.i = method[i], :
## Some levels of 'A3' are removed since no observation in that/those levels## Warning in compare.i(X[, i], y = y, selec.i = selec[i], method.i = method[i], :
## Some levels of 'A4' are removed since no observation in that/those levels## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect## Warning in compare.i(X[, i], y = y, selec.i = selec[i], method.i = method[i], :
## Some levels of 'C6' are removed since no observation in that/those levelst1##
##
## -------- Summary of results by groups of 'C10.1'---------
##
##
## var N p.value
## 1 NhomNC 206 <0.001**
## 2 A2 206 0.770
## 3 A3 206 0.169
## 4 A4 206 0.405
## 5 A5 206 0.007**
## 6 A6.1 206 0.791
## 7 A6.2 206 0.396
## 8 A6.3 206 0.475
## 9 A7.1 206 0.001**
## 10 A7.2 206 0.003**
## 11 A7.3 206 0.001**
## 12 A8.1 206 0.030**
## 13 A8.2 206 0.225
## 14 A8.3 205 0.012**
## 15 A8.4 206 <0.001**
## 16 Thuc hanh TKV 206 <0.001**
## 17 Thoi gian th TKV 206 0.885
## 18 Tu the coi ao boc lo nguc truoc guong hoac ke goi duoi vai 206 <0.001**
## 19 Quan sat hinh the 206 <0.001**
## 20 Quan sat da vung vu 206 0.002**
## 21 Quan sat num vu 206 0.004**
## 22 So sanh 2 vu 205 <0.001**
## 23 Danh gia thay doi huong num vu 206 <0.001**
## 24 Chong tay quan sat 19 .
## 25 Kiem tra tiet dich 206 <0.001**
## 26 Kham hach 205 <0.001**
## 27 ep ngon tay de len mo vu 206 <0.001**
## 28 Luc de thay doi 206 <0.001**
## 29 Dùng tay T khám vú P và tay P khám vú T 206 <0.001**
## 30 Xoay deu ko bo sot phan mo vu 206 <0.001**
## 31 C10.1 206 <0.001**
## 32 C10.2 204 <0.001**
## 33 thuc hien so nan vu 11 .
## 34 C10.1 206 <0.001**
## 35 C10.2 204 <0.001**
## 36 Tuoi 206 0.606
## 37 Tong diem nhan thuc nhay cam 206 <0.001**
## 38 Tong diem muc do nghiem trong 206 0.980
## 39 Tong diem loi ich TKV 206 0.003**
## 40 Tong diem rao can 206 <0.001**
## 41 Tong diem dong luc TKV 206 0.128
## 42 Tong diem hieu qua ban than 206 <0.001**
## 43 Tong diem raocan TKV 206 <0.001**
## method selection
## 1 categorical ALL
## 2 categorical ALL
## 3 categorical ALL
## 4 categorical ALL
## 5 categorical ALL
## 6 categorical ALL
## 7 categorical ALL
## 8 categorical ALL
## 9 categorical ALL
## 10 categorical ALL
## 11 categorical ALL
## 12 categorical ALL
## 13 categorical ALL
## 14 categorical ALL
## 15 categorical ALL
## 16 categorical ALL
## 17 categorical ALL
## 18 categorical ALL
## 19 categorical ALL
## 20 categorical ALL
## 21 categorical ALL
## 22 categorical ALL
## 23 categorical ALL
## 24 categorical ALL
## 25 categorical ALL
## 26 categorical ALL
## 27 categorical ALL
## 28 categorical ALL
## 29 categorical ALL
## 30 categorical ALL
## 31 categorical ALL
## 32 categorical ALL
## 33 categorical ALL
## 34 categorical ALL
## 35 categorical ALL
## 36 continuous normal ALL
## 37 continuous normal ALL
## 38 continuous normal ALL
## 39 continuous normal ALL
## 40 continuous normal ALL
## 41 continuous normal ALL
## 42 continuous normal ALL
## 43 continuous normal ALL
## -----
## Signif. codes: 0 '**' 0.05 '*' 0.1 ' ' 1t1a0 <- createTable(t1, hide.no = "no", show.n = TRUE)
print(t1a0)##
## --------Summary descriptives table by 'C10.1'---------
##
## _________________________________________________________________________________________________
## 0 1 p.overall N
## N=111 N=95
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯
## NhomNC: <0.001 206
## 1 26 (23.4%) 49 (51.6%)
## 2 69 (62.2%) 24 (25.3%)
## 3 16 (14.4%) 22 (23.2%)
## A2: 0.770 206
## 1 1 (0.90%) 0 (0.00%)
## 2 82 (73.9%) 73 (76.8%)
## 3 2 (1.80%) 3 (3.16%)
## 4 26 (23.4%) 19 (20.0%)
## A3: 0.169 206
## 1 0 (0.00%) 1 (1.05%)
## 2 108 (97.3%) 89 (93.7%)
## 3 3 (2.70%) 2 (2.11%)
## 4 0 (0.00%) 3 (3.16%)
## A4: 0.405 206
## 1 2 (1.80%) 2 (2.11%)
## 2 27 (24.3%) 22 (23.2%)
## 3 51 (45.9%) 48 (50.5%)
## 5 19 (17.1%) 18 (18.9%)
## 6 2 (1.80%) 3 (3.16%)
## 8 10 (9.01%) 2 (2.11%)
## A5: 0.007 206
## 1 28 (25.2%) 8 (8.42%)
## 2 19 (17.1%) 13 (13.7%)
## 3 26 (23.4%) 34 (35.8%)
## 4 38 (34.2%) 40 (42.1%)
## A6.1: 0.791 206
## 1 5 (4.50%) 6 (6.32%)
## 2 106 (95.5%) 89 (93.7%)
## A6.2: 0.396 206
## 1 12 (10.8%) 15 (15.8%)
## 2 99 (89.2%) 80 (84.2%)
## A6.3: 0.475 206
## 1 3 (2.70%) 5 (5.26%)
## 2 108 (97.3%) 90 (94.7%)
## A7.1: 0.001 206
## 1 31 (27.9%) 49 (51.6%)
## 2 80 (72.1%) 46 (48.4%)
## A7.2: 0.003 206
## 1 24 (21.6%) 40 (42.1%)
## 2 87 (78.4%) 55 (57.9%)
## A7.3: 0.001 206
## 1 75 (67.6%) 40 (42.1%)
## 2 1 (0.90%) 5 (5.26%)
## 3 2 (1.80%) 10 (10.5%)
## 4 2 (1.80%) 5 (5.26%)
## 5 2 (1.80%) 2 (2.11%)
## 6 1 (0.90%) 5 (5.26%)
## 7 28 (25.2%) 28 (29.5%)
## A8.1: 0.030 206
## 1 25 (22.5%) 32 (33.7%)
## 2 11 (9.91%) 16 (16.8%)
## 3 75 (67.6%) 47 (49.5%)
## A8.2: 0.225 206
## 1 13 (11.7%) 19 (20.0%)
## 2 4 (3.60%) 4 (4.21%)
## 3 94 (84.7%) 72 (75.8%)
## A8.3: 0.012 205
## 1 94 (85.5%) 66 (69.5%)
## 2 3 (2.73%) 3 (3.16%)
## 3 13 (11.8%) 26 (27.4%)
## A8.4: <0.001 206
## 1 25 (22.5%) 47 (49.5%)
## 2 10 (9.01%) 19 (20.0%)
## 3 76 (68.5%) 29 (30.5%)
## Thuc hanh TKV: <0.001 206
## khong TKV 47 (42.3%) 0 (0.00%)
## co TKV 64 (57.7%) 95 (100%)
## Thoi gian th TKV: 0.885 206
## Không ?úng th?i ?i?m 85 (76.6%) 71 (74.7%)
## khám ?úng th?i ?i?m 7-10 ngay 26 (23.4%) 24 (25.3%)
## Tu the coi ao boc lo nguc truoc guong hoac ke goi duoi vai: <0.001 206
## khong 69 (62.2%) 33 (34.7%)
## co 42 (37.8%) 62 (65.3%)
## Quan sat hinh the: <0.001 206
## khong 65 (58.6%) 29 (30.5%)
## co 46 (41.4%) 66 (69.5%)
## Quan sat da vung vu: 0.002 206
## khong 64 (57.7%) 33 (34.7%)
## co 47 (42.3%) 62 (65.3%)
## Quan sat num vu: 0.004 206
## khong 63 (56.8%) 34 (35.8%)
## co 48 (43.2%) 61 (64.2%)
## So sanh 2 vu: <0.001 205
## Khong 73 (65.8%) 17 (18.1%)
## Co 38 (34.2%) 77 (81.9%)
## Danh gia thay doi huong num vu: <0.001 206
## Khong 98 (88.3%) 36 (37.9%)
## Co 13 (11.7%) 59 (62.1%)
## Chong tay quan sat: Co 7 (100%) 12 (100%) . 19
## Kiem tra tiet dich: <0.001 206
## Khong 93 (83.8%) 48 (50.5%)
## Co 18 (16.2%) 47 (49.5%)
## Kham hach: <0.001 205
## khong 97 (88.2%) 55 (57.9%)
## Co 13 (11.8%) 40 (42.1%)
## ep ngon tay de len mo vu: <0.001 206
## Khong 96 (86.5%) 28 (29.5%)
## Co 15 (13.5%) 67 (70.5%)
## Luc de thay doi: <0.001 206
## khong thay doi luc de 104 (93.7%) 28 (29.5%)
## co thay doi luc de 7 (6.31%) 67 (70.5%)
## Dùng tay T khám vú P và tay P khám vú T: <0.001 206
## khong dung tay nay kham vu kia 77 (69.4%) 38 (40.0%)
## co dung tay nay kham vu kia 34 (30.6%) 57 (60.0%)
## Xoay deu ko bo sot phan mo vu: <0.001 206
## khong xoay deu cac mo vu 100 (90.1%) 57 (60.0%)
## co xoay deu cac mo vu 11 (9.91%) 38 (40.0%)
## C10.1: <0.001 206
## 0 111 (100%) 0 (0.00%)
## 1 0 (0.00%) 95 (100%)
## C10.2: <0.001 204
## 0 109 (100%) 30 (31.6%)
## 1 0 (0.00%) 65 (68.4%)
## thuc hien so nan vu: . 11
## khong TH 9 (81.8%) 0 (.%)
## Co TH 2 (18.2%) 0 (.%)
## C10.1: <0.001 206
## 0 111 (100%) 0 (0.00%)
## 1 0 (0.00%) 95 (100%)
## C10.2: <0.001 204
## 0 109 (100%) 30 (31.6%)
## 1 0 (0.00%) 65 (68.4%)
## Tuoi 46.5 (9.64) 47.1 (8.41) 0.606 206
## Tong diem nhan thuc nhay cam 9.75 (5.27) 13.0 (5.52) <0.001 206
## Tong diem muc do nghiem trong 22.8 (7.43) 22.9 (7.26) 0.980 206
## Tong diem loi ich TKV 20.1 (9.89) 23.4 (5.59) 0.003 206
## Tong diem rao can 10.9 (6.39) 14.1 (5.93) <0.001 206
## Tong diem dong luc TKV 27.0 (7.27) 28.4 (5.41) 0.128 206
## Tong diem hieu qua ban than 22.1 (13.7) 30.6 (13.5) <0.001 206
## Tong diem raocan TKV 10.9 (6.39) 14.1 (5.93) <0.001 206
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯1. Doctors BSE & Factors
t1a=compareGroups(C10.2~C10.1+NhomNC + A2 + A3 + A4 + A5 + A6.1 + A6.2 + A6.3 + A7.1 + A7.2 + A7.3 + A8.1 + A8.2 + A8.3 + A8.4 + C10.1 + C10.2 + A1 + TB1_5 + TB6_12 + TB13_18 + TB + TB25_31 + TB32_42 + TB19_24 , data=dd2)## Warning in compare.i(X[, i], y = y, selec.i = selec[i], method.i = method[i], :
## Some levels of 'A3' are removed since no observation in that/those levels## Warning in compare.i(X[, i], y = y, selec.i = selec[i], method.i = method[i], :
## Some levels of 'A4' are removed since no observation in that/those levels## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrect
## Warning in chisq.test(xx, correct = FALSE): Chi-squared approximation may be
## incorrectt1a##
##
## -------- Summary of results by groups of 'C10.2'---------
##
##
## var N p.value method selection
## 1 C10.1 204 <0.001** categorical ALL
## 2 NhomNC 204 0.001** categorical ALL
## 3 A2 204 0.757 categorical ALL
## 4 A3 204 0.554 categorical ALL
## 5 A4 204 0.485 categorical ALL
## 6 A5 204 0.182 categorical ALL
## 7 A6.1 204 0.746 categorical ALL
## 8 A6.2 204 0.625 categorical ALL
## 9 A6.3 204 0.269 categorical ALL
## 10 A7.1 204 0.473 categorical ALL
## 11 A7.2 204 0.078* categorical ALL
## 12 A7.3 204 0.058* categorical ALL
## 13 A8.1 204 0.106 categorical ALL
## 14 A8.2 204 0.528 categorical ALL
## 15 A8.3 203 0.003** categorical ALL
## 16 A8.4 204 0.007** categorical ALL
## 17 C10.1 204 <0.001** categorical ALL
## 18 C10.2 204 <0.001** categorical ALL
## 19 Tuoi 204 0.183 continuous normal ALL
## 20 Tong diem nhan thuc nhay cam 204 0.102 continuous normal ALL
## 21 Tong diem muc do nghiem trong 204 0.873 continuous normal ALL
## 22 Tong diem loi ich TKV 204 0.045** continuous normal ALL
## 23 Tong diem rao can 204 <0.001** continuous normal ALL
## 24 Tong diem dong luc TKV 204 0.210 continuous normal ALL
## 25 Tong diem hieu qua ban than 204 <0.001** continuous normal ALL
## 26 Tong diem raocan TKV 204 <0.001** continuous normal ALL
## -----
## Signif. codes: 0 '**' 0.05 '*' 0.1 ' ' 1t1a1=createTable(t1a)
print(t1a1)##
## --------Summary descriptives table by 'C10.2'---------
##
## _______________________________________________________________
## 0 1 p.overall
## N=139 N=65
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯
## C10.1: <0.001
## 0 109 (78.4%) 0 (0.00%)
## 1 30 (21.6%) 65 (100%)
## NhomNC: 0.001
## 1 40 (28.8%) 33 (50.8%)
## 2 76 (54.7%) 17 (26.2%)
## 3 23 (16.5%) 15 (23.1%)
## A2: 0.757
## 1 1 (0.72%) 0 (0.00%)
## 2 102 (73.4%) 51 (78.5%)
## 3 3 (2.16%) 2 (3.08%)
## 4 33 (23.7%) 12 (18.5%)
## A3: 0.554
## 1 1 (0.72%) 0 (0.00%)
## 2 134 (96.4%) 61 (93.8%)
## 3 3 (2.16%) 2 (3.08%)
## 4 1 (0.72%) 2 (3.08%)
## A4: 0.485
## 1 3 (2.16%) 1 (1.54%)
## 2 35 (25.2%) 14 (21.5%)
## 3 60 (43.2%) 37 (56.9%)
## 5 28 (20.1%) 9 (13.8%)
## 6 3 (2.16%) 2 (3.08%)
## 8 10 (7.19%) 2 (3.08%)
## A5: 0.182
## 1 30 (21.6%) 6 (9.23%)
## 2 19 (13.7%) 11 (16.9%)
## 3 38 (27.3%) 22 (33.8%)
## 4 52 (37.4%) 26 (40.0%)
## A6.1: 0.746
## 1 7 (5.04%) 4 (6.15%)
## 2 132 (95.0%) 61 (93.8%)
## A6.2: 0.625
## 1 20 (14.4%) 7 (10.8%)
## 2 119 (85.6%) 58 (89.2%)
## A6.3: 0.269
## 1 4 (2.88%) 4 (6.15%)
## 2 135 (97.1%) 61 (93.8%)
## A7.1: 0.473
## 1 51 (36.7%) 28 (43.1%)
## 2 88 (63.3%) 37 (56.9%)
## A7.2: 0.078
## 1 37 (26.6%) 26 (40.0%)
## 2 102 (73.4%) 39 (60.0%)
## A7.3: 0.058
## 1 84 (60.4%) 30 (46.2%)
## 2 4 (2.88%) 2 (3.08%)
## 3 8 (5.76%) 4 (6.15%)
## 4 3 (2.16%) 4 (6.15%)
## 5 2 (1.44%) 2 (3.08%)
## 6 1 (0.72%) 5 (7.69%)
## 7 37 (26.6%) 18 (27.7%)
## A8.1: 0.106
## 1 33 (23.7%) 24 (36.9%)
## 2 17 (12.2%) 9 (13.8%)
## 3 89 (64.0%) 32 (49.2%)
## A8.2: 0.528
## 1 19 (13.7%) 13 (20.0%)
## 2 6 (4.32%) 2 (3.08%)
## 3 114 (82.0%) 50 (76.9%)
## A8.3: 0.003
## 1 117 (84.8%) 42 (64.6%)
## 2 4 (2.90%) 2 (3.08%)
## 3 17 (12.3%) 21 (32.3%)
## A8.4: 0.007
## 1 40 (28.8%) 30 (46.2%)
## 2 17 (12.2%) 12 (18.5%)
## 3 82 (59.0%) 23 (35.4%)
## C10.1: <0.001
## 0 109 (78.4%) 0 (0.00%)
## 1 30 (21.6%) 65 (100%)
## C10.2: <0.001
## 0 139 (100%) 0 (0.00%)
## 1 0 (0.00%) 65 (100%)
## Tuoi 46.1 (9.33) 47.9 (8.50) 0.183
## Tong diem nhan thuc nhay cam 10.8 (5.63) 12.1 (5.47) 0.102
## Tong diem muc do nghiem trong 22.9 (7.15) 22.7 (7.81) 0.873
## Tong diem loi ich TKV 20.9 (9.20) 23.1 (5.96) 0.045
## Tong diem rao can 11.2 (6.01) 14.9 (6.50) <0.001
## Tong diem dong luc TKV 27.4 (6.98) 28.5 (5.13) 0.210
## Tong diem hieu qua ban than 23.5 (13.6) 31.5 (14.3) <0.001
## Tong diem raocan TKV 11.2 (6.01) 14.9 (6.50) <0.001
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯2. Factors and Group of study
t2=compareGroups(NhomNC~ A2 + A3 + A4 + A5 + A6.1 + A6.2 + A6.3 + A7.1 + A7.2 + A7.3 + A8.1 + A8.2 + A8.3 + A8.4 + C0+ C1+ C2 + C3.1 + C3.2 + C3.3 + C4 + C5 + C6 + C7 +C8 +C9.1+C9.2+C9.3+C9.4+C10.1+C10.2+C9+C10.1 + C10.2 + A1 + TB1_5 + TB6_12 + TB13_18 + TB + TB25_31 + TB32_42 + TB19_24 , data=dd2)## Warning in compare.i(X[, i], y = y, selec.i = selec[i], method.i = method[i], :
## Some levels of 'C6' are removed since no observation in that/those levels## Warning in cor(as.integer(x), as.integer(y)): the standard deviation is zero
## Warning in cor(as.integer(x), as.integer(y)): the standard deviation is zerot2##
##
## -------- Summary of results by groups of 'NhomNC'---------
##
##
## var N p.value
## 1 A2 352 0.668
## 2 A3 352 0.218
## 3 A4 352 .
## 4 A5 352 0.019**
## 5 A6.1 352 0.351
## 6 A6.2 352 0.507
## 7 A6.3 352 0.862
## 8 A7.1 352 <0.001**
## 9 A7.2 352 0.013**
## 10 A7.3 352 .
## 11 A8.1 352 0.678
## 12 A8.2 352 0.931
## 13 A8.3 351 0.004**
## 14 A8.4 352 <0.001**
## 15 Thuc hanh TKV 352 0.039**
## 16 Thoi gian th TKV 207 0.014**
## 17 Tu the coi ao boc lo nguc truoc guong hoac ke goi duoi vai 207 <0.001**
## 18 Quan sat hinh the 207 <0.001**
## 19 Quan sat da vung vu 207 0.003**
## 20 Quan sat num vu 207 0.004**
## 21 So sanh 2 vu 206 <0.001**
## 22 Danh gia thay doi huong num vu 207 0.009**
## 23 Chong tay quan sat 20 .
## 24 Kiem tra tiet dich 207 <0.001**
## 25 Kham hach 206 0.011**
## 26 ep ngon tay de len mo vu 206 <0.001**
## 27 Luc de thay doi 206 <0.001**
## 28 Dùng tay T khám vú P và tay P khám vú T 206 <0.001**
## 29 Xoay deu ko bo sot phan mo vu 206 <0.001**
## 30 C10.1 206 <0.001**
## 31 C10.2 204 0.001**
## 32 thuc hien so nan vu 11 .
## 33 C10.1 206 <0.001**
## 34 C10.2 204 0.001**
## 35 Tuoi 352 0.156
## 36 Tong diem nhan thuc nhay cam 352 <0.001**
## 37 Tong diem muc do nghiem trong 352 0.003**
## 38 Tong diem loi ich TKV 351 <0.001**
## 39 Tong diem rao can 352 <0.001**
## 40 Tong diem dong luc TKV 352 0.001**
## 41 Tong diem hieu qua ban than 352 <0.001**
## 42 Tong diem raocan TKV 352 <0.001**
## method selection
## 1 categorical ALL
## 2 categorical ALL
## 3 categorical ALL
## 4 categorical ALL
## 5 categorical ALL
## 6 categorical ALL
## 7 categorical ALL
## 8 categorical ALL
## 9 categorical ALL
## 10 categorical ALL
## 11 categorical ALL
## 12 categorical ALL
## 13 categorical ALL
## 14 categorical ALL
## 15 categorical ALL
## 16 categorical ALL
## 17 categorical ALL
## 18 categorical ALL
## 19 categorical ALL
## 20 categorical ALL
## 21 categorical ALL
## 22 categorical ALL
## 23 categorical ALL
## 24 categorical ALL
## 25 categorical ALL
## 26 categorical ALL
## 27 categorical ALL
## 28 categorical ALL
## 29 categorical ALL
## 30 categorical ALL
## 31 categorical ALL
## 32 categorical ALL
## 33 categorical ALL
## 34 categorical ALL
## 35 continuous normal ALL
## 36 continuous normal ALL
## 37 continuous normal ALL
## 38 continuous normal ALL
## 39 continuous normal ALL
## 40 continuous normal ALL
## 41 continuous normal ALL
## 42 continuous normal ALL
## -----
## Signif. codes: 0 '**' 0.05 '*' 0.1 ' ' 1t2a=createTable(t2)
print(t2a)##
## --------Summary descriptives table by 'NhomNC'---------
##
## _________________________________________________________________________________________________________
## 1 2 3 p.overall
## N=142 N=104 N=106
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯
## A2: 0.668
## 1 1 (0.70%) 0 (0.00%) 1 (0.94%)
## 2 95 (66.9%) 74 (71.2%) 80 (75.5%)
## 3 6 (4.23%) 2 (1.92%) 2 (1.89%)
## 4 40 (28.2%) 28 (26.9%) 23 (21.7%)
## A3: 0.218
## 1 5 (3.52%) 1 (0.96%) 2 (1.89%)
## 2 135 (95.1%) 98 (94.2%) 97 (91.5%)
## 3 1 (0.70%) 3 (2.88%) 5 (4.72%)
## 4 0 (0.00%) 2 (1.92%) 1 (0.94%)
## 5 1 (0.70%) 0 (0.00%) 1 (0.94%)
## A4: .
## 1 3 (2.11%) 1 (0.96%) 3 (2.83%)
## 2 43 (30.3%) 30 (28.8%) 23 (21.7%)
## 3 55 (38.7%) 52 (50.0%) 48 (45.3%)
## 4 1 (0.70%) 0 (0.00%) 0 (0.00%)
## 5 27 (19.0%) 12 (11.5%) 24 (22.6%)
## 6 8 (5.63%) 0 (0.00%) 2 (1.89%)
## 8 5 (3.52%) 9 (8.65%) 6 (5.66%)
## A5: 0.019
## 1 14 (9.86%) 25 (24.0%) 14 (13.2%)
## 2 23 (16.2%) 16 (15.4%) 9 (8.49%)
## 3 51 (35.9%) 25 (24.0%) 36 (34.0%)
## 4 54 (38.0%) 38 (36.5%) 47 (44.3%)
## A6.1: 0.351
## 1 5 (3.52%) 6 (5.77%) 2 (1.89%)
## 2 137 (96.5%) 98 (94.2%) 104 (98.1%)
## A6.2: 0.507
## 1 15 (10.6%) 14 (13.5%) 9 (8.49%)
## 2 127 (89.4%) 90 (86.5%) 97 (91.5%)
## A6.3: 0.862
## 1 4 (2.82%) 4 (3.85%) 3 (2.83%)
## 2 138 (97.2%) 100 (96.2%) 103 (97.2%)
## A7.1: <0.001
## 1 66 (46.5%) 25 (24.0%) 28 (26.4%)
## 2 76 (53.5%) 79 (76.0%) 78 (73.6%)
## A7.2: 0.013
## 1 45 (31.7%) 18 (17.3%) 20 (18.9%)
## 2 97 (68.3%) 86 (82.7%) 86 (81.1%)
## A7.3: .
## 1 70 (49.3%) 75 (72.1%) 84 (79.2%)
## 2 8 (5.63%) 1 (0.96%) 4 (3.77%)
## 3 13 (9.15%) 1 (0.96%) 2 (1.89%)
## 4 9 (6.34%) 0 (0.00%) 1 (0.94%)
## 5 3 (2.11%) 2 (1.92%) 1 (0.94%)
## 6 5 (3.52%) 2 (1.92%) 2 (1.89%)
## 7 34 (23.9%) 23 (22.1%) 12 (11.3%)
## A8.1: 0.678
## 1 37 (26.1%) 25 (24.0%) 30 (28.3%)
## 2 19 (13.4%) 9 (8.65%) 10 (9.43%)
## 3 86 (60.6%) 70 (67.3%) 66 (62.3%)
## A8.2: 0.931
## 1 18 (12.7%) 14 (13.5%) 15 (14.2%)
## 2 4 (2.82%) 4 (3.85%) 5 (4.72%)
## 3 120 (84.5%) 86 (82.7%) 86 (81.1%)
## A8.3: 0.004
## 1 99 (70.2%) 91 (87.5%) 85 (80.2%)
## 2 8 (5.67%) 0 (0.00%) 1 (0.94%)
## 3 34 (24.1%) 13 (12.5%) 20 (18.9%)
## A8.4: <0.001
## 1 80 (56.3%) 18 (17.3%) 27 (25.5%)
## 2 24 (16.9%) 16 (15.4%) 18 (17.0%)
## 3 38 (26.8%) 70 (67.3%) 61 (57.5%)
## Thuc hanh TKV: 0.039
## khong TKV 68 (47.9%) 57 (54.8%) 68 (64.2%)
## co TKV 74 (52.1%) 47 (45.2%) 38 (35.8%)
## Thoi gian th TKV: 0.014
## Không ?úng th?i ?i?m 54 (71.1%) 79 (84.9%) 24 (63.2%)
## khám ?úng th?i ?i?m 7-10 ngay 22 (28.9%) 14 (15.1%) 14 (36.8%)
## Tu the coi ao boc lo nguc truoc guong hoac ke goi duoi vai: <0.001
## khong 26 (34.2%) 64 (68.8%) 13 (34.2%)
## co 50 (65.8%) 29 (31.2%) 25 (65.8%)
## Quan sat hinh the: <0.001
## khong 28 (36.8%) 57 (61.3%) 10 (26.3%)
## co 48 (63.2%) 36 (38.7%) 28 (73.7%)
## Quan sat da vung vu: 0.003
## khong 33 (43.4%) 54 (58.1%) 10 (26.3%)
## co 43 (56.6%) 39 (41.9%) 28 (73.7%)
## Quan sat num vu: 0.004
## khong 34 (44.7%) 54 (58.1%) 10 (26.3%)
## co 42 (55.3%) 39 (41.9%) 28 (73.7%)
## So sanh 2 vu: <0.001
## Khong 18 (23.7%) 59 (63.4%) 14 (37.8%)
## Co 58 (76.3%) 34 (36.6%) 23 (62.2%)
## Danh gia thay doi huong num vu: 0.009
## Khong 40 (52.6%) 70 (75.3%) 25 (65.8%)
## Co 36 (47.4%) 23 (24.7%) 13 (34.2%)
## Chong tay quan sat: Co 11 (100%) 4 (100%) 5 (100%) .
## Kiem tra tiet dich: <0.001
## Khong 39 (51.3%) 75 (80.6%) 28 (73.7%)
## Co 37 (48.7%) 18 (19.4%) 10 (26.3%)
## Kham hach: 0.011
## khong 50 (66.7%) 78 (83.9%) 24 (63.2%)
## Co 25 (33.3%) 15 (16.1%) 14 (36.8%)
## ep ngon tay de len mo vu: <0.001
## Khong 28 (37.3%) 78 (83.9%) 18 (47.4%)
## Co 47 (62.7%) 15 (16.1%) 20 (52.6%)
## Luc de thay doi: <0.001
## khong thay doi luc de 33 (44.0%) 77 (82.8%) 22 (57.9%)
## co thay doi luc de 42 (56.0%) 16 (17.2%) 16 (42.1%)
## Dùng tay T khám vú P và tay P khám vú T: <0.001
## khong dung tay nay kham vu kia 33 (44.0%) 70 (75.3%) 12 (31.6%)
## co dung tay nay kham vu kia 42 (56.0%) 23 (24.7%) 26 (68.4%)
## Xoay deu ko bo sot phan mo vu: <0.001
## khong xoay deu cac mo vu 45 (60.0%) 80 (86.0%) 32 (84.2%)
## co xoay deu cac mo vu 30 (40.0%) 13 (14.0%) 6 (15.8%)
## C10.1: <0.001
## 0 26 (34.7%) 69 (74.2%) 16 (42.1%)
## 1 49 (65.3%) 24 (25.8%) 22 (57.9%)
## C10.2: 0.001
## 0 40 (54.8%) 76 (81.7%) 23 (60.5%)
## 1 33 (45.2%) 17 (18.3%) 15 (39.5%)
## thuc hien so nan vu: .
## khong TH 0 (.%) 9 (81.8%) 0 (.%)
## Co TH 0 (.%) 2 (18.2%) 0 (.%)
## C10.1: <0.001
## 0 26 (34.7%) 69 (74.2%) 16 (42.1%)
## 1 49 (65.3%) 24 (25.8%) 22 (57.9%)
## C10.2: 0.001
## 0 40 (54.8%) 76 (81.7%) 23 (60.5%)
## 1 33 (45.2%) 17 (18.3%) 15 (39.5%)
## Tuoi 44.8 (9.57) 46.5 (9.93) 47.0 (9.80) 0.156
## Tong diem nhan thuc nhay cam 13.1 (4.80) 9.82 (5.65) 10.0 (5.26) <0.001
## Tong diem muc do nghiem trong 24.5 (5.52) 22.0 (7.90) 21.8 (7.50) 0.003
## Tong diem loi ich TKV 23.7 (4.90) 20.0 (9.96) 20.4 (9.05) <0.001
## Tong diem rao can 14.6 (5.23) 10.6 (5.98) 12.6 (7.28) <0.001
## Tong diem dong luc TKV 29.1 (4.55) 26.8 (6.78) 26.2 (7.90) 0.001
## Tong diem hieu qua ban than 32.7 (11.5) 21.2 (13.9) 22.2 (12.4) <0.001
## Tong diem raocan TKV 14.6 (5.23) 10.6 (5.98) 12.9 (8.04) <0.001
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯export2xls(t2a, file = "C:/Users/binht/Dropbox/Hue/data/2022_data_consultant/2 dd2/R/example.xlsx")export2xls(t1a1, file = "C:/Users/binht/Dropbox/Hue/data/2022_data_consultant/2 dd2/R/t1a1.xlsx")export2xls(t1a0, file = "C:/Users/binht/Dropbox/Hue/data/2022_data_consultant/2 dd2/R/t1a0.xlsx")- Logistic model to determine the associated with BSE practice
library(epiDisplay)## Loading required package: foreign## Loading required package: survival## Loading required package: MASS## Loading required package: nnetm1= glm(C10.1 ~ NhomNC* TB1_5 + TB6_12 + TB32_42, family=binomial, data=dd2)summary(m1)##
## Call:
## glm(formula = C10.1 ~ NhomNC * TB1_5 + TB6_12 + TB32_42, family = binomial,
## data = dd2)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.7923 -0.9496 -0.4696 0.9694 2.0614
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -0.21858 0.93808 -0.233 0.8158
## NhomNC2 -1.86627 0.90863 -2.054 0.0400 *
## NhomNC3 0.28264 1.10930 0.255 0.7989
## TB1_5 0.06226 0.05226 1.191 0.2335
## TB6_12 -0.03040 0.02282 -1.333 0.1827
## TB32_42 0.02373 0.01248 1.902 0.0572 .
## NhomNC2:TB1_5 0.04898 0.06693 0.732 0.4643
## NhomNC3:TB1_5 -0.03218 0.08786 -0.366 0.7141
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 284.33 on 205 degrees of freedom
## Residual deviance: 238.71 on 198 degrees of freedom
## (146 observations deleted due to missingness)
## AIC: 254.71
##
## Number of Fisher Scoring iterations: 4logistic.display(m1)##
## Logistic regression predicting C10.1 : 1 vs 0
##
## crude OR(95%CI) adj. OR(95%CI) P(Wald's test)
## NhomNC: ref.=1
## 2 0.18 (0.09,0.36) 0.15 (0.03,0.92) 0.04
## 3 0.73 (0.33,1.62) 1.33 (0.15,11.67) 0.799
##
## TB1_5 (cont. var.) 1.12 (1.06,1.18) 1.06 (0.96,1.18) 0.234
##
## TB6_12 (cont. var.) 1 (0.96,1.04) 0.97 (0.93,1.01) 0.183
##
## TB32_42 (cont. var.) 1.05 (1.02,1.07) 1.02 (1,1.05) 0.057
##
## NhomNC:TB1_5: ref.=1:
## 2: - 1.05 (0.92,1.2) 0.464
## 3: - 0.97 (0.82,1.15) 0.714
##
## P(LR-test)
## NhomNC: ref.=1 0.029
## 2
## 3
##
## TB1_5 (cont. var.) 1
##
## TB6_12 (cont. var.) 0.179
##
## TB32_42 (cont. var.) 0.057
##
## NhomNC:TB1_5: ref.=1: 0.56
## 2:
## 3:
##
## Log-likelihood = -119.3554
## No. of observations = 206
## AIC value = 254.7107#dd2a <- na.omit(dd2) mydata1=subset(dd2, select=c("TB1_5" , "TB6_12", "TB13_18" , "TB" , "TB25_31" , "TB32_42"))library("Hmisc")## Loading required package: lattice##
## Attaching package: 'lattice'## The following object is masked from 'package:epiDisplay':
##
## dotplot## Loading required package: Formula## Loading required package: ggplot2##
## Attaching package: 'ggplot2'## The following object is masked from 'package:epiDisplay':
##
## alpha##
## Attaching package: 'Hmisc'## The following objects are masked from 'package:table1':
##
## label, label<-, units## The following objects are masked from 'package:base':
##
## format.pval, unitslibrary('PerformanceAnalytics')## Loading required package: xts## Loading required package: zoo##
## Attaching package: 'zoo'## The following objects are masked from 'package:base':
##
## as.Date, as.Date.numeric##
## Attaching package: 'PerformanceAnalytics'## The following object is masked from 'package:graphics':
##
## legendres2 <- rcorr(as.matrix(mydata1))res2## TB1_5 TB6_12 TB13_18 TB TB25_31 TB32_42
## TB1_5 1.00 0.28 0.32 0.20 0.22 0.34
## TB6_12 0.28 1.00 0.23 -0.13 0.29 0.05
## TB13_18 0.32 0.23 1.00 0.11 0.64 0.47
## TB 0.20 -0.13 0.11 1.00 0.02 0.32
## TB25_31 0.22 0.29 0.64 0.02 1.00 0.39
## TB32_42 0.34 0.05 0.47 0.32 0.39 1.00
##
## n
## TB1_5 TB6_12 TB13_18 TB TB25_31 TB32_42
## TB1_5 352 352 351 352 352 352
## TB6_12 352 352 351 352 352 352
## TB13_18 351 351 351 351 351 351
## TB 352 352 351 352 352 352
## TB25_31 352 352 351 352 352 352
## TB32_42 352 352 351 352 352 352
##
## P
## TB1_5 TB6_12 TB13_18 TB TB25_31 TB32_42
## TB1_5 0.0000 0.0000 0.0002 0.0000 0.0000
## TB6_12 0.0000 0.0000 0.0174 0.0000 0.3569
## TB13_18 0.0000 0.0000 0.0455 0.0000 0.0000
## TB 0.0002 0.0174 0.0455 0.7431 0.0000
## TB25_31 0.0000 0.0000 0.0000 0.7431 0.0000
## TB32_42 0.0000 0.3569 0.0000 0.0000 0.0000chart.Correlation(mydata1, histogram=TRUE, pch=25)
LS0tDQp0aXRsZTogJ3Byb2plY3QgMjogQlNFJw0KYXV0aG9yOiAiQmluaCBUaGFuZyINCmRhdGU6ICcyMDIyLTEwLTI1Jw0KDQpvdXRwdXQ6DQogIGh0bWxfZG9jdW1lbnQ6DQogICAgY29kZV9kb3dubG9hZDogeWVzDQogICAgY29kZV9mb2xkaW5nOiBoaWRlDQogICAgbnVtYmVyX3NlY3Rpb25zOiB5ZXMNCiAgICB0aGVtZTogam91cm5hbA0KICAgIHRvYzogeWVzDQogICAgdG9jX2Zsb2F0OiB5ZXMNCiAgd29yZF9kb2N1bWVudDoNCiAgICB0b2M6IHllcw0KLS0tDQoNCmBgYHtyfQ0KbGlicmFyeShoYXZlbikNCg0KDQpgYGANCg0KYGBge3J9DQpkZDI9cmVhZF9zYXYoIkM6L1VzZXJzL2Jpbmh0L0Ryb3Bib3gvSHVlL2RhdGEvMjAyMl9kYXRhX2NvbnN1bHRhbnQvMiBkZDIvUi9zcHNzX3Iuc2F2IikNCmBgYA0KDQpgYGB7cn0NCmNsYXNzKGRkMikNCmBgYA0KDQpgYGB7cn0NCmRhdGEubGFiZWwudGFibGUgPC0gYXR0cihkZDIsICJsYWJlbC50YWJsZSIpDQptaXNzaW5ncyA8LSBhdHRyKGRkMiwgIm1pc3NpbmdzIikNCmBgYA0KDQpgYGB7cn0NCm5hbWVzKGRkMikNCmBgYA0KDQpgYGB7cn0NCmxpYnJhcnkoY29tcGFyZUdyb3VwcykNCmxpYnJhcnkodGFibGUxKQ0KYGBgDQoNCmBgYHtyfQ0KZGQyJE5ob21OQz1hcy5mYWN0b3IoZGQyJE5ob21OQykNCmRkMiRDMTAuMT1hcy5mYWN0b3IoZGQyJEMxMC4xKQ0KZGQyJEEyPWFzLmZhY3RvcihkZDIkQTIpDQpkZDIkQTM9YXMuZmFjdG9yKGRkMiRBMykNCmRkMiRBND1hcy5mYWN0b3IoZGQyJEE0KQ0KZGQyJEE1PWFzLmZhY3RvcihkZDIkQTUpDQpkZDIkQTYuMT1hcy5mYWN0b3IoZGQyJEE2LjEpDQpkZDIkQTYuMj1hcy5mYWN0b3IoZGQyJEE2LjIpDQpkZDIkQTYuMz1hcy5mYWN0b3IoZGQyJEE2LjMpDQpkZDIkQTcuMT1hcy5mYWN0b3IoZGQyJEE3LjEpDQpkZDIkQTcuMj1hcy5mYWN0b3IoZGQyJEE3LjIpDQpkZDIkQTcuMz1hcy5mYWN0b3IoZGQyJEE3LjMpDQpkZDIkQTguMT1hcy5mYWN0b3IoZGQyJEE4LjEpDQpkZDIkQTguMj1hcy5mYWN0b3IoZGQyJEE4LjIpDQpkZDIkQTguMz1hcy5mYWN0b3IoZGQyJEE4LjMpDQpkZDIkQTguND1hcy5mYWN0b3IoZGQyJEE4LjQpDQpkZDIkQzEwLjI9YXMuZmFjdG9yKGRkMiRDMTAuMikNCg0KYGBgDQoNCioqMC4gRGVzY3JpcHRpb24qKg0KDQpgYGB7cn0NCmEwPXRhYmxlMShOaG9tTkN+QzEwLjEgKyBBMiArIEEzICsgIEE0ICArIEE1ICsgIEE2LjEgKyAgQTYuMiArICBBNi4zICsgIEE3LjEgICsgQTcuMiArICBBNy4zICsgIEE4LjEgKyAgQTguMiArICBBOC4zICsgIEE4LjQgKyAgIEMxMC4xICArIEMxMC4yICArIEExICsgVEIxXzUgKyAgVEI2XzEyICsgIFRCMTNfMTggKyAgVEIgKyAgVEIyNV8zMSAgKyBUQjMyXzQyICArIFRCMTlfMjQsIGRhdGE9ZGQyKQ0KYGBgDQoNCmBgYHtyfQ0KYTANCmBgYA0KDQoqKjEuIFJlZ3VsYXIgQlNFICYgRmFjdG9ycyoqDQoNCmBgYHtyfQ0KdDE9Y29tcGFyZUdyb3VwcyhDMTAuMX5OaG9tTkMgKyBBMiArIEEzICsgIEE0ICArIEE1ICsgIEE2LjEgKyAgQTYuMiArICBBNi4zICsgIEE3LjEgICsgQTcuMiArICBBNy4zICsgIEE4LjEgKyAgQTguMiArICBBOC4zICsgIEE4LjQgKyBDMCsgQzErIEMyICsgQzMuMSArIEMzLjIgKyBDMy4zICsgQzQgKyBDNSArIEM2ICsgQzcgK0M4ICtDOS4xK0M5LjIrQzkuMytDOS40K0MxMC4xK0MxMC4yK0M5K0MxMC4xICArIEMxMC4yICArIEExICsgVEIxXzUgKyAgVEI2XzEyICsgIFRCMTNfMTggKyAgVEIgKyAgVEIyNV8zMSAgKyBUQjMyXzQyICArIFRCMTlfMjQgLCBkYXRhPWRkMikNCmBgYA0KDQpgYGB7cn0NCnQxDQpgYGANCg0KYGBge3J9DQoNCg0KdDFhMCA8LSBjcmVhdGVUYWJsZSh0MSwgIGhpZGUubm8gPSAibm8iLCBzaG93Lm4gPSBUUlVFKQ0KDQpwcmludCh0MWEwKQ0KDQpgYGANCg0KKioxLiBEb2N0b3JzIEJTRSAmIEZhY3RvcnMqKg0KDQpgYGB7cn0NCnQxYT1jb21wYXJlR3JvdXBzKEMxMC4yfkMxMC4xK05ob21OQyArIEEyICsgQTMgKyAgQTQgICsgQTUgKyAgQTYuMSArICBBNi4yICsgIEE2LjMgKyAgQTcuMSAgKyBBNy4yICsgIEE3LjMgKyAgQTguMSArICBBOC4yICsgIEE4LjMgKyAgQTguNCArICBDMTAuMSAgKyBDMTAuMiAgKyBBMSArIFRCMV81ICsgIFRCNl8xMiArICBUQjEzXzE4ICsgIFRCICsgIFRCMjVfMzEgICsgVEIzMl80MiAgKyBUQjE5XzI0ICwgZGF0YT1kZDIpDQpgYGANCg0KYGBge3J9DQp0MWENCg0KYGBgDQoNCmBgYHtyfQ0KdDFhMT1jcmVhdGVUYWJsZSh0MWEpDQpwcmludCh0MWExKQ0KYGBgDQoNCioqMi4gRmFjdG9ycyBhbmQgR3JvdXAgb2Ygc3R1ZHkqKg0KDQpgYGB7cn0NCnQyPWNvbXBhcmVHcm91cHMoTmhvbU5DfiBBMiArIEEzICsgIEE0ICArIEE1ICsgIEE2LjEgKyAgQTYuMiArICBBNi4zICsgIEE3LjEgICsgQTcuMiArICBBNy4zICsgIEE4LjEgKyAgQTguMiArICBBOC4zICsgIEE4LjQgKyBDMCsgQzErIEMyICsgQzMuMSArIEMzLjIgKyBDMy4zICsgQzQgKyBDNSArIEM2ICsgQzcgK0M4ICtDOS4xK0M5LjIrQzkuMytDOS40K0MxMC4xK0MxMC4yK0M5K0MxMC4xICArIEMxMC4yICArIEExICsgVEIxXzUgKyAgVEI2XzEyICsgIFRCMTNfMTggKyAgVEIgKyAgVEIyNV8zMSAgKyBUQjMyXzQyICArIFRCMTlfMjQgLCBkYXRhPWRkMikNCmBgYA0KDQpgYGB7cn0NCnQyDQpgYGANCg0KYGBge3J9DQp0MmE9Y3JlYXRlVGFibGUodDIpDQpwcmludCh0MmEpDQpgYGANCg0KYGBge3J9DQpleHBvcnQyeGxzKHQyYSwgZmlsZSA9ICJDOi9Vc2Vycy9iaW5odC9Ecm9wYm94L0h1ZS9kYXRhLzIwMjJfZGF0YV9jb25zdWx0YW50LzIgZGQyL1IvZXhhbXBsZS54bHN4IikNCmBgYA0KDQpgYGB7cn0NCmV4cG9ydDJ4bHModDFhMSwgZmlsZSA9ICJDOi9Vc2Vycy9iaW5odC9Ecm9wYm94L0h1ZS9kYXRhLzIwMjJfZGF0YV9jb25zdWx0YW50LzIgZGQyL1IvdDFhMS54bHN4IikNCmBgYA0KDQpgYGB7cn0NCmV4cG9ydDJ4bHModDFhMCwgZmlsZSA9ICJDOi9Vc2Vycy9iaW5odC9Ecm9wYm94L0h1ZS9kYXRhLzIwMjJfZGF0YV9jb25zdWx0YW50LzIgZGQyL1IvdDFhMC54bHN4IikNCmBgYA0KDQo0LiAgKipMb2dpc3RpYyBtb2RlbCB0byBkZXRlcm1pbmUgdGhlIGFzc29jaWF0ZWQgd2l0aCBCU0UgcHJhY3RpY2UqKg0KDQpgYGB7cn0NCmxpYnJhcnkoZXBpRGlzcGxheSkNCmBgYA0KDQpgYGB7cn0NCm0xPSBnbG0oQzEwLjEgfiBOaG9tTkMqIFRCMV81ICsgIFRCNl8xMiAgKyAgIFRCMzJfNDIsIGZhbWlseT1iaW5vbWlhbCwgZGF0YT1kZDIpDQoNCmBgYA0KDQpgYGB7cn0NCnN1bW1hcnkobTEpDQpgYGANCg0KYGBge3J9DQpsb2dpc3RpYy5kaXNwbGF5KG0xKQ0KYGBgDQoNCmBgYHtyfQ0KI2RkMmEgPC0gbmEub21pdChkZDIpIA0KYGBgDQoNCmBgYHtyfQ0KbXlkYXRhMT1zdWJzZXQoZGQyLCBzZWxlY3Q9YygiVEIxXzUiICwgICJUQjZfMTIiLCAgIlRCMTNfMTgiICwgICJUQiIgLCAgIlRCMjVfMzEiICAsICJUQjMyXzQyIikpDQoNCmBgYA0KDQpgYGB7cn0NCmxpYnJhcnkoIkhtaXNjIikNCmxpYnJhcnkoJ1BlcmZvcm1hbmNlQW5hbHl0aWNzJykNCmBgYA0KDQpgYGB7cn0NCnJlczIgPC0gcmNvcnIoYXMubWF0cml4KG15ZGF0YTEpKQ0KYGBgDQoNCmBgYHtyfQ0KcmVzMg0KYGBgDQoNCmBgYHtyfQ0KY2hhcnQuQ29ycmVsYXRpb24obXlkYXRhMSwgaGlzdG9ncmFtPVRSVUUsIHBjaD0yNSkNCg0KYGBgDQo=