t = "C:\\Users\\Admin\\Desktop\\tonghop.csv"
duy = read.csv(t)
head(duy)## id sex age weigh time_day time_month c0 dose_mg dose_kg mmf glucose
## 1 1 1 37 54 1 0 3.2 6.0 0.1111 2.0 6.8
## 2 3 0 24 39 1 0 7.2 5.0 0.1282 1.5 11.3
## 3 4 1 30 52 1 0 18.1 6.0 0.1154 2.0 6.6
## 4 5 1 44 68 1 0 16.4 5.5 0.0809 2.0 8.7
## 5 6 1 41 62 1 0 10.8 7.0 0.1129 2.0 8.3
## 6 7 1 36 53 1 0 6.3 6.0 0.1132 2.0 7.4
## creatinin egfr alb cho hdl ldl triglycerid uric bil_toltal
## 1 453 13.5593 31.4 NA NA NA NA NA 8.0
## 2 128 47.2302 34.9 NA NA NA NA NA 8.6
## 3 606 10.1132 28.9 3.4 NA NA NA NA 8.7
## 4 546 10.5530 33.2 4.1 NA NA NA NA 8.5
## 5 438 13.8058 35.5 NA NA NA NA NA 6.6
## 6 205 34.0428 41.3 NA NA NA NA NA 7.6
## bil_directed ast alt amylase ggt hct hb
## 1 1.2 24 15 74 46 0.281 92
## 2 2.5 19 14 142 24 0.295 100
## 3 1.5 22 13 217 18 0.246 80
## 4 2.7 13 15 41 59 0.292 96
## 5 0.1 15 13 87 20 0.351 112
## 6 2.0 12 11 105 NA 0.248 81library("BMA")## Loading required package: survival## Loading required package: leaps## Loading required package: robustbase##
## Attaching package: 'robustbase'## The following object is masked from 'package:survival':
##
## heart## Loading required package: inline## Loading required package: rrcov## Scalable Robust Estimators with High Breakdown Point (version 1.4-7)t1 = subset(duy, time_day == 1)
t3 = subset(duy, time_day == 3)
t7 = subset(duy, time_day == 7)
t14 = subset(duy, time_day == 14)
t30 = subset(duy, time_day == 30)
t90 = subset(duy, time_day == 90)
t180 = subset(duy, time_day == 180)
t360 = subset(duy, time_day == 360)Phân tích tại thời điểm T1
y1 = t1[,7]
xvars1= t1[,c(2,3,4,8:27)]
mh1 = bicreg(xvars1, y1 , strict = F, OR = 20)## Warning in bicreg(xvars1, y1, strict = F, OR = 20): There were 92 records
## deleted due to NA's## Warning in leaps.setup(x, y, wt = weights, nbest = nbest, nvmax = nvmax, :
## 6 linear dependencies found## Warning in leaps.setup(x, y, wt = weights, nbest = nbest, nvmax = nvmax, :
## nvmax reduced to 17## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs producedsummary(mh1)##
## Call:
## bicreg(x = xvars1, y = y1, strict = F, OR = 20)
##
##
## 309 models were selected
## Best 5 models (cumulative posterior probability = 0.1068 ):
##
## p!=0 EV SD model 1 model 2
## Intercept 100.0 -2.485e+01 132.75256 -22.09475 -25.28753
## sex 62.7 -2.185e+00 14.44582 . .
## age 93.5 1.073e+00 0.55337 1.32231 1.31366
## weigh 41.8 -1.059e-01 0.65831 . .
## dose_mg 30.5 1.822e-01 7.35310 0.89484 .
## dose_kg 33.6 -1.518e+02 448.82784 . .
## mmf 41.8 3.480e-02 40.09703 . .
## glucose 72.1 1.121e+00 0.86950 1.46059 1.55221
## creatinin 80.4 -6.215e-02 0.09224 -0.10419 -0.10400
## egfr 77.7 -1.607e+00 1.51226 -2.15426 -2.22956
## alb 70.2 2.682e-01 1.34811 0.72687 0.57984
## cho 48.4 2.648e+00 5.61658 . 1.79294
## hdl 27.8 2.852e+00 12.91649 . .
## ldl 54.6 -2.917e+00 11.00065 . -1.55030
## triglycerid 46.3 -5.233e+00 16.68896 . .
## uric 53.2 1.161e-02 0.02948 . .
## bil_toltal 70.4 -9.019e-01 1.88412 -2.17052 -1.83493
## bil_directed 80.1 -1.550e-02 3.66925 2.52964 1.84675
## ast 73.8 4.229e-01 1.14304 0.67022 0.62702
## alt 43.2 3.791e-03 0.45121 . .
## amylase 78.0 -1.331e-01 0.13265 -0.16672 -0.17199
## ggt 86.3 -1.192e-01 0.26542 -0.21774 -0.20457
## hct 91.6 1.698e+02 142.32547 210.10660 192.29421
## hb 46.5 1.418e-01 0.87431 -0.15772 .
##
## nVar 13 13
## r2 0.999 0.999
## BIC -86.76476 -86.76476
## post prob 0.021 0.021
## model 3 model 4 model 5
## Intercept -288.33401 -86.33111 169.91959
## sex -26.25025 -16.55037 21.19928
## age 2.29058 0.72344 0.86384
## weigh . 0.22585 -0.60934
## dose_mg . . .
## dose_kg -893.88654 . .
## mmf 102.11262 . -49.52776
## glucose . . 2.26023
## creatinin . 0.02978 -0.18262
## egfr . . -3.68463
## alb 3.56866 . -0.88922
## cho . 9.07336 .
## hdl . 16.98412 .
## ldl . -22.47585 7.72421
## triglycerid -16.20476 . .
## uric -0.03210 0.02942 .
## bil_toltal -3.06051 . .
## bil_directed 4.85382 -0.81182 -2.02002
## ast 2.72839 . 0.39745
## alt . -0.57537 .
## amylase -0.46993 -0.11378 .
## ggt -0.61007 0.20324 -0.02757
## hct 180.23480 212.71458 187.38037
## hb . . .
##
## nVar 13 13 13
## r2 0.999 0.999 0.999
## BIC -86.76476 -86.76476 -86.76476
## post prob 0.021 0.021 0.021Xây dựng mô hình tại T1
m1 = lm(c0 ~ age + dose_mg + glucose + creatinin + egfr + alb + bil_toltal + bil_directed + ast + amylase + ggt + hct + hb, data = t1)
summary(m1)##
## Call:
## lm(formula = c0 ~ age + dose_mg + glucose + creatinin + egfr +
## alb + bil_toltal + bil_directed + ast + amylase + ggt + hct +
## hb, data = t1)
##
## Residuals:
## Min 1Q Median 3Q Max
## -8.749 -3.503 -1.005 2.334 10.447
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -9.561358 8.283047 -1.154 0.25176
## age 0.013387 0.062212 0.215 0.83016
## dose_mg -0.475848 0.745732 -0.638 0.52521
## glucose -0.087854 0.237434 -0.370 0.71234
## creatinin 0.014419 0.005278 2.732 0.00773 **
## egfr 0.204806 0.093469 2.191 0.03131 *
## alb 0.002322 0.142747 0.016 0.98706
## bil_toltal -0.371006 0.259013 -1.432 0.15588
## bil_directed 1.307167 0.715407 1.827 0.07136 .
## ast -0.005193 0.028317 -0.183 0.85496
## amylase -0.001902 0.002617 -0.727 0.46949
## ggt -0.013971 0.018584 -0.752 0.45436
## hct 67.462697 27.383066 2.464 0.01587 *
## hb -0.060446 0.095281 -0.634 0.52761
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 5.033 on 81 degrees of freedom
## (15 observations deleted due to missingness)
## Multiple R-squared: 0.2994, Adjusted R-squared: 0.1869
## F-statistic: 2.662 on 13 and 81 DF, p-value: 0.003742Do có các yếu tố không ảnh hưởng tới C0 (p< 0.05) nên xây dựng lại mô hình với các yếu tố có p<0.05
m1 = lm(c0 ~ creatinin + egfr + hct, data = t1)
summary(m1)##
## Call:
## lm(formula = c0 ~ creatinin + egfr + hct, data = t1)
##
## Residuals:
## Min 1Q Median 3Q Max
## -9.8205 -3.6507 -0.8899 2.3432 12.5159
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -12.85651 4.37418 -2.939 0.004047 **
## creatinin 0.01589 0.00466 3.411 0.000921 ***
## egfr 0.24847 0.08169 3.042 0.002970 **
## hct 37.61112 10.57973 3.555 0.000568 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 5.037 on 105 degrees of freedom
## (1 observation deleted due to missingness)
## Multiple R-squared: 0.2114, Adjusted R-squared: 0.1888
## F-statistic: 9.38 on 3 and 105 DF, p-value: 1.509e-05Phân tích tại thời điểm T3
y3 = t3[,7]
xvars3= t3[,c(2,3,4,8:27)]
mh3 = bicreg(xvars3, y3 , strict = F, OR = 20)## Warning in bicreg(xvars3, y3, strict = F, OR = 20): There were 91 records
## deleted due to NA's## Warning in leaps.setup(x, y, wt = weights, nbest = nbest, nvmax = nvmax, :
## 4 linear dependencies found## Warning in leaps.setup(x, y, wt = weights, nbest = nbest, nvmax = nvmax, :
## nvmax reduced to 19## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs producedsummary(mh3)##
## Call:
## bicreg(x = xvars3, y = y3, strict = F, OR = 20)
##
##
## 217 models were selected
## Best 5 models (cumulative posterior probability = 0.0921 ):
##
## p!=0 EV SD model 1 model 2
## Intercept 100.0 -5.930e+00 1.917e+01 -9.427e+00 -3.520e+01
## sex 68.4 4.131e+00 3.852e+00 8.608e+00 .
## age 84.9 4.174e-01 3.794e-01 8.351e-01 -4.297e-01
## weigh 77.5 2.145e-01 2.047e-01 3.843e-01 .
## dose_mg 46.8 1.954e-01 1.945e+00 . .
## dose_kg 63.5 7.916e-01 1.114e+02 6.540e+01 -2.108e+02
## mmf 32.4 1.663e+00 3.390e+00 . 1.011e+01
## glucose 75.9 -1.241e+00 1.048e+00 -2.182e+00 -1.488e+00
## creatinin 69.4 -6.154e-02 6.533e-02 -1.308e-01 .
## egfr 48.3 -7.954e-03 5.533e-02 . .
## alb 59.4 -4.667e-01 5.512e-01 . -3.634e-01
## cho 65.2 5.329e+00 5.668e+00 . 2.006e+01
## hdl 61.6 2.504e+00 3.025e+00 4.814e+00 -2.980e+00
## ldl 99.5 -1.642e+01 5.256e+00 -1.379e+01 -3.101e+01
## triglycerid 100.0 7.143e+00 1.485e+00 8.721e+00 6.188e+00
## uric 88.4 2.763e-02 3.915e-02 6.928e-02 .
## bil_toltal 65.5 1.598e+00 1.548e+00 . 5.187e+00
## bil_directed 87.8 -2.425e+00 2.672e+00 -6.950e-01 -9.690e+00
## ast 59.9 -3.684e-01 4.259e-01 -8.717e-01 4.901e-01
## alt 83.4 8.005e-02 2.857e-01 4.060e-01 -3.141e-01
## amylase 62.6 4.712e-02 6.471e-02 1.245e-01 -9.597e-02
## ggt 11.8 9.054e-04 4.335e-03 . .
## hct 53.2 -2.582e+01 3.809e+01 . .
## hb 64.7 2.217e-01 2.093e-01 . 5.073e-01
##
## nVar 14 15
## r2 0.999 0.999
## BIC -9.621e+01 -9.322e+01
## post prob 0.049 0.011
## model 3 model 4 model 5
## Intercept 1.975e+01 1.511e+01 6.957e+00
## sex . . .
## age 1.364e-01 1.486e-01 1.430e-01
## weigh . . 1.830e-01
## dose_mg 1.829e+00 1.904e+00 .
## dose_kg -4.519e+01 -5.353e+01 4.062e+01
## mmf . . .
## glucose . . .
## creatinin -7.623e-03 . .
## egfr -7.868e-02 -3.922e-02 -5.158e-02
## alb -1.227e+00 -1.151e+00 -1.198e+00
## cho 9.072e+00 7.566e+00 7.968e+00
## hdl . 1.739e+00 1.468e+00
## ldl -1.718e+01 -1.552e+01 -1.581e+01
## triglycerid 5.686e+00 5.772e+00 5.700e+00
## uric -2.127e-02 -1.844e-02 -2.002e-02
## bil_toltal 2.640e+00 2.594e+00 2.663e+00
## bil_directed -2.336e+00 -2.344e+00 -2.441e+00
## ast . . .
## alt -2.432e-01 -2.338e-01 -2.260e-01
## amylase . . .
## ggt . . .
## hct -8.934e+01 -7.475e+01 -7.309e+01
## hb 4.890e-01 4.223e-01 4.257e-01
##
## nVar 15 15 15
## r2 0.999 0.999 0.999
## BIC -9.322e+01 -9.322e+01 -9.322e+01
## post prob 0.011 0.011 0.011Xây dựng mô hình tại T3
m3 = lm(c0 ~ sex + age + dose_mg + creatinin + cho + hdl + ldl + uric + bil_directed + ast + alt + amylase , data = t3)
summary(m3)##
## Call:
## lm(formula = c0 ~ sex + age + dose_mg + creatinin + cho + hdl +
## ldl + uric + bil_directed + ast + alt + amylase, data = t3)
##
## Residuals:
## Min 1Q Median 3Q Max
## -5.2757 -1.8422 -0.3094 0.9070 6.5261
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 29.398564 20.168405 1.458 0.1789
## sex 1.922053 3.235818 0.594 0.5671
## age 0.485463 0.192105 2.527 0.0324 *
## dose_mg -4.041691 2.315494 -1.745 0.1149
## creatinin -0.019867 0.046478 -0.427 0.6791
## cho -2.490504 4.521924 -0.551 0.5952
## hdl -4.532999 6.386442 -0.710 0.4958
## ldl -0.436873 5.035662 -0.087 0.9328
## uric 0.021962 0.019464 1.128 0.2884
## bil_directed 0.571191 1.172160 0.487 0.6377
## ast -0.262322 0.326413 -0.804 0.4423
## alt 0.239060 0.178528 1.339 0.2134
## amylase 0.001579 0.035787 0.044 0.9658
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 4.035 on 9 degrees of freedom
## (89 observations deleted due to missingness)
## Multiple R-squared: 0.7144, Adjusted R-squared: 0.3337
## F-statistic: 1.876 on 12 and 9 DF, p-value: 0.1752Xây dựng lại mô hình tại T3
m3 = lm(c0 ~ age , data = t3)
summary(m3)##
## Call:
## lm(formula = c0 ~ age, data = t3)
##
## Residuals:
## Min 1Q Median 3Q Max
## -7.619 -3.896 -1.767 2.399 23.028
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 6.04118 2.20203 2.743 0.00711 **
## age 0.07893 0.06072 1.300 0.19632
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 5.524 on 109 degrees of freedom
## Multiple R-squared: 0.01527, Adjusted R-squared: 0.006235
## F-statistic: 1.69 on 1 and 109 DF, p-value: 0.1963#Không có ảnh hưởng của tuổi tới C0Phân tích dữ liệu tại T7
y7 = t7[,7]
xvars7= t7[,c(2,3,4,8:27)]
mh7 = bicreg(xvars7, y7 , strict = F, OR = 20)## Warning in bicreg(xvars7, y7, strict = F, OR = 20): There were 72 records
## deleted due to NA's## Warning in leaps.setup(x, y, wt = weights, nbest = nbest, nvmax = nvmax, :
## 19 linear dependencies found## Reordering variables and trying again:## Warning in log(vr): NaNs produced## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs producedsummary(mh7)##
## Call:
## bicreg(x = xvars7, y = y7, strict = F, OR = 20)
##
##
## 303 models were selected
## Best 5 models (cumulative posterior probability = 0.0264 ):
##
## p!=0 EV SD model 1 model 2
## Intercept 100.0 1.425e+00 NaN -10.683972 -10.966852
## sex 43.0 9.594e-01 NaN . .
## age 26.1 6.995e-02 NaN 0.236646 0.236752
## weigh 21.5 -4.701e-02 NaN . .
## dose_mg 7.5 4.656e-03 0.6224 . .
## dose_kg 18.4 -7.280e+00 NaN . .
## mmf 6.0 4.949e-01 NaN . .
## glucose 10.6 -1.654e-02 NaN . .
## creatinin 11.3 1.882e-03 NaN . .
## egfr 8.2 -4.529e-04 NaN . .
## alb 15.2 4.469e-03 NaN 0.171419 0.168233
## cho 9.0 -7.123e-02 NaN . .
## hdl 8.6 5.575e-01 NaN . .
## ldl 10.7 -1.011e-01 NaN . .
## triglycerid 40.2 3.722e-01 NaN . .
## uric 15.1 -3.037e-07 NaN . .
## bil_toltal 9.1 4.389e-02 NaN . .
## bil_directed 8.3 -7.727e-04 NaN . .
## ast 12.2 1.690e-02 NaN . .
## alt 11.6 9.906e-03 NaN . .
## amylase 6.6 -6.065e-03 NaN -0.004419 .
## ggt 11.4 1.647e-03 NaN . 0.003936
## hct 7.1 3.185e+01 NaN . .
## hb 17.7 -9.443e-02 NaN . .
##
## nVar 3 3
## r2 0.999 0.999
## BIC -29.710463 -29.710463
## post prob 0.005 0.005
## model 3 model 4 model 5
## Intercept -6.272960 -2.628954 -3.217650
## sex . . .
## age 0.142084 0.168839 0.327735
## weigh . . .
## dose_mg . . .
## dose_kg . . .
## mmf . . .
## glucose . . .
## creatinin . . -0.022968
## egfr . . .
## alb 0.106978 . .
## cho . . .
## hdl . . .
## ldl . . .
## triglycerid 0.485686 . .
## uric . . .
## bil_toltal . . .
## bil_directed . . -1.238080
## ast . . .
## alt . 0.120087 .
## amylase . -0.018885 .
## ggt . . .
## hct . . .
## hb . . .
##
## nVar 3 3 3
## r2 0.999 0.999 0.999
## BIC -29.710463 -29.710463 -29.710463
## post prob 0.005 0.005 0.005Xây dựng mô hình tại T7
m7 = lm(c0 ~ age +egfr + amylase, data = t7)
summary(m7)##
## Call:
## lm(formula = c0 ~ age + egfr + amylase, data = t7)
##
## Residuals:
## Min 1Q Median 3Q Max
## -5.5541 -2.9142 -0.8849 3.0862 6.8586
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 5.273255 3.233260 1.631 0.1124
## age 0.151008 0.078056 1.935 0.0616 .
## egfr -0.044555 0.023543 -1.892 0.0672 .
## amylase -0.006461 0.017073 -0.378 0.7075
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 3.745 on 33 degrees of freedom
## (40 observations deleted due to missingness)
## Multiple R-squared: 0.1718, Adjusted R-squared: 0.09652
## F-statistic: 2.282 on 3 and 33 DF, p-value: 0.09735#Không có yếu tố ảnh hưởngPhân tích dữ liệu tại T30
y30 = t30[,7]
xvars30= t30[,c(2,3,4,8:27)]
mh30 = bicreg(xvars30, y30 , strict = F, OR = 20)## Warning in bicreg(xvars30, y30, strict = F, OR = 20): There were 73 records
## deleted due to NA's## Warning in leaps.setup(x, y, wt = weights, nbest = nbest, nvmax = nvmax, :
## 18 linear dependencies found## Warning in leaps.setup(x, y, wt = weights, nbest = nbest, nvmax = nvmax, :
## nvmax reduced to 5## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs producedsummary(mh30)##
## Call:
## bicreg(x = xvars30, y = y30, strict = F, OR = 20)
##
##
## 306 models were selected
## Best 5 models (cumulative posterior probability = 0.0235 ):
##
## p!=0 EV SD model 1 model 2 model 3
## Intercept 100.0 12.3119129 NaN 16.89951 45.92082 47.48075
## sex 50.4 -3.2000834 NaN . . .
## age 37.2 -0.0037299 NaN -0.38162 -0.39990 -0.21968
## weigh 27.7 0.0435369 NaN . . .
## dose_mg 24.2 -0.0092179 NaN . . .
## dose_kg 17.3 -9.0438031 NaN . . .
## mmf 19.1 1.0989448 NaN . . .
## glucose 10.3 -0.0386237 NaN . . .
## creatinin 10.6 -0.0007234 NaN . . -0.11027
## egfr 11.9 -0.0086790 NaN . . -0.33347
## alb 16.1 -0.0640626 NaN . . .
## cho 18.0 -0.0092642 NaN . -6.18066 .
## hdl 16.9 -1.0560713 NaN . . .
## ldl 11.2 0.0095851 NaN . 3.07908 .
## triglycerid 10.7 0.0909010 NaN . . .
## uric 17.2 0.0016439 NaN . . .
## bil_toltal 20.5 -0.0402806 NaN . -0.24875 .
## bil_directed 24.9 -0.1982019 NaN . . .
## ast 13.8 0.0045963 NaN 0.42632 . 0.08587
## alt 12.9 -0.0019898 NaN -0.10824 . .
## amylase 12.8 -0.0051824 NaN . . .
## ggt 14.2 -0.0015633 NaN -0.03985 . .
## hct 13.9 -0.2899715 NaN . . .
## hb 16.2 0.0019355 NaN . . .
##
## nVar 4 4 4
## r2 0.999 0.999 0.999
## BIC -34.27949 -34.27949 -34.27949
## post prob 0.005 0.005 0.005
## model 4 model 5
## Intercept 42.86685 20.82166
## sex . .
## age -0.47364 -0.46534
## weigh . .
## dose_mg . .
## dose_kg -86.83042 -76.59576
## mmf . .
## glucose . .
## creatinin . .
## egfr . -0.07705
## alb . .
## cho -1.66825 .
## hdl . .
## ldl . .
## triglycerid . .
## uric . 0.04820
## bil_toltal . .
## bil_directed . .
## ast . .
## alt . .
## amylase . .
## ggt -0.01649 .
## hct . .
## hb . .
##
## nVar 4 4
## r2 0.999 0.999
## BIC -34.27949 -34.27949
## post prob 0.005 0.005Xây dựng mô hình T30
m30 = lm(c0 ~ age + ast + alt + amylase, data = t30)
summary(m30)##
## Call:
## lm(formula = c0 ~ age + ast + alt + amylase, data = t30)
##
## Residuals:
## Min 1Q Median 3Q Max
## -3.7169 -2.4803 -1.6459 0.8533 10.1653
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.41286 7.66692 -0.184 0.857
## age 0.23947 0.22835 1.049 0.317
## ast 0.37251 0.24361 1.529 0.154
## alt -0.13427 0.13383 -1.003 0.337
## amylase -0.03035 0.03990 -0.761 0.463
##
## Residual standard error: 4.662 on 11 degrees of freedom
## (63 observations deleted due to missingness)
## Multiple R-squared: 0.3902, Adjusted R-squared: 0.1685
## F-statistic: 1.76 on 4 and 11 DF, p-value: 0.2071#Không có yếu tố ảnh hưởngPhân tích dữ liệu tại T90
y90 = t90[,7]
xvars90= t90[,c(2,3,4,8:27)]
mh90 = bicreg(xvars90, y90 , strict = F, OR = 20)## Warning in bicreg(xvars90, y90, strict = F, OR = 20): There were 59 records
## deleted due to NA's## Warning in leaps.setup(x, y, wt = weights, nbest = nbest, nvmax = nvmax, :
## 21 linear dependencies found## Warning in leaps.setup(x, y, wt = weights, nbest = nbest, nvmax = nvmax, :
## nvmax reduced to 2## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs produced
## Warning in log(vr): NaNs producedsummary(mh90)##
## Call:
## bicreg(x = xvars90, y = y90, strict = F, OR = 20)
##
##
## 147 models were selected
## Best 5 models (cumulative posterior probability = 0.0342 ):
##
## p!=0 EV SD model 1 model 2 model 3
## Intercept 100.0 66.053604 NaN 165.7079 25.5473 13.2507
## sex 1.4 -0.094923 NaN . . .
## age 14.4 -0.118720 NaN 0.1950 -0.4397 0.2992
## weigh 8.2 -0.046164 NaN . . .
## dose_mg 10.3 -0.742080 NaN . . .
## dose_kg 10.3 154.596280 NaN . . .
## mmf 6.2 3.132492 NaN . . .
## glucose 6.8 -1.124936 NaN . . .
## creatinin 9.6 0.026278 NaN . . .
## egfr 6.2 -0.007835 NaN . . .
## alb 13.7 -0.523329 NaN -3.8529 . .
## cho 13.7 0.537027 NaN . . .
## hdl 8.2 -9.327963 NaN . . .
## ldl 12.3 -0.499129 NaN . . .
## triglycerid 8.2 -0.788113 NaN . . .
## uric 8.9 -0.057279 NaN . . .
## bil_toltal 8.2 -0.040439 NaN . . .
## bil_directed 6.2 -0.264778 NaN . . .
## ast 10.3 -0.037497 NaN . . .
## alt 6.8 -0.011301 NaN . -0.1013 .
## amylase 9.6 -0.006902 NaN . . -0.1460
## ggt 5.5 -0.005308 NaN . . .
## hct 7.6 -9.723192 NaN . . .
## hb 7.5 -0.040698 NaN . . .
##
## nVar 2 2 2
## r2 0.999 0.999 0.999
## BIC -18.5260 -18.5260 -18.5260
## post prob 0.007 0.007 0.007
## model 4 model 5
## Intercept 29.7193 -0.2154
## sex . .
## age -0.4837 0.4182
## weigh . .
## dose_mg . .
## dose_kg . .
## mmf . .
## glucose . .
## creatinin . .
## egfr . .
## alb . .
## cho . .
## hdl . .
## ldl . .
## triglycerid . .
## uric . .
## bil_toltal . .
## bil_directed . -2.9274
## ast -0.2368 .
## alt . .
## amylase . .
## ggt . .
## hct . .
## hb . .
##
## nVar 2 2
## r2 0.999 0.999
## BIC -18.5260 -18.5260
## post prob 0.007 0.007Xây dựng mô hình tại T90
m90 = lm(c0 ~ age + alb, data = t90)
summary(m90)##
## Call:
## lm(formula = c0 ~ age + alb, data = t90)
##
## Residuals:
## Min 1Q Median 3Q Max
## -5.415 -2.296 -1.073 1.910 11.655
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.61503 6.66177 0.393 0.6961
## age 0.10778 0.05547 1.943 0.0569 .
## alb 0.02888 0.13199 0.219 0.8275
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 3.702 on 58 degrees of freedom
## (1 observation deleted due to missingness)
## Multiple R-squared: 0.06161, Adjusted R-squared: 0.02926
## F-statistic: 1.904 on 2 and 58 DF, p-value: 0.1581#Không có yếu tố ảnh hưởng