Workflow_WTA (original) (raw)
Running the DE model with default values.
NBthDEmod <- fitNBthDE(form = ~region,
split = FALSE,
object = kidney)
#> Iteration = 1, squared error = 1.545817e-05
#> Iteration = 2, squared error = 1.218597e-05
str(NBthDEmod)
#> List of 12
#> $ X : num [1:16, 1:2] 1 1 1 1 1 1 1 1 1 1 ...
#> ..- attr(*, "dimnames")=List of 2
#> .. ..$ : chr [1:16] "DSP-1001250007851-H-A02.dcc" "DSP-1001250007851-H-A03.dcc" "DSP-1001250007851-H-A04.dcc" "DSP-1001250007851-H-A05.dcc" ...
#> .. ..$ : chr [1:2] "(Intercept)" "regiontubule"
#> ..- attr(*, "assign")= int [1:2] 0 1
#> ..- attr(*, "contrasts")=List of 1
#> .. ..$ region: chr "contr.treatment"
#> $ para0 : num [1:4, 1:1500] 5.75 2.35 23.13 181 7.41 ...
#> ..- attr(*, "dimnames")=List of 2
#> .. ..$ : chr [1:4] "(Intercept)" "regiontubule" "r" "threshold"
#> .. ..$ : chr [1:1500] "DDT" "FOSL1" "FKBP11" "PIM1" ...
#> $ para : num [1:4, 1:18439] 12.38 -1.58 3.12 181 6.61 ...
#> ..- attr(*, "dimnames")=List of 2
#> .. ..$ : chr [1:4] "(Intercept)" "regiontubule" "r" "threshold"
#> .. ..$ : chr [1:18439] "A2M" "NAT2" "ACADM" "ACADS" ...
#> $ sizefact : Named num [1:16] 0.0683 0.0457 0.0426 0.0549 0.0729 ...
#> ..- attr(*, "names")= chr [1:16] "DSP-1001250007851-H-A02.dcc" "DSP-1001250007851-H-A03.dcc" "DSP-1001250007851-H-A04.dcc" "DSP-1001250007851-H-A05.dcc" ...
#> $ sizefact0 : Named num [1:16] 0.0671 0.0452 0.0427 0.0548 0.073 ...
#> ..- attr(*, "names")= chr [1:16] "DSP-1001250007851-H-A02.dcc" "DSP-1001250007851-H-A03.dcc" "DSP-1001250007851-H-A04.dcc" "DSP-1001250007851-H-A05.dcc" ...
#> $ preci1 : num [1:2, 1:2] 0.04 0.02 0.02 1.4
#> ..- attr(*, "dimnames")=List of 2
#> .. ..$ : chr [1:2] "(Intercept)" "regiontubule"
#> .. ..$ : chr [1:2] "(Intercept)" "regiontubule"
#> $ Im0 :List of 1500
#> ..$ DDT : num[0 , 0 ]
#> ..$ FOSL1 : num[0 , 0 ]
#> ..$ FKBP11 : num[0 , 0 ]
#> ..$ PIM1 : num[0 , 0 ]
#> ..$ RABEP2 : num[0 , 0 ]
#> ..$ ASAP2 : num[0 , 0 ]
#> ..$ IRX2 : num[0 , 0 ]
#> ..$ CENPX : num[0 , 0 ]
#> ..$ NDUFS2 : num[0 , 0 ]
#> ..$ MED11 : num[0 , 0 ]
#> ..$ PDP1 : num[0 , 0 ]
#> ..$ DCAF11 : num[0 , 0 ]
#> ..$ GLCE : num[0 , 0 ]
#> ..$ ELAVL1 : num[0 , 0 ]
#> ..$ SSH2 : num[0 , 0 ]
#> ..$ CFHR1 : num[0 , 0 ]
#> ..$ SEC23B : num[0 , 0 ]
#> ..$ DGCR8 : num[0 , 0 ]
#> ..$ PVALEF : num[0 , 0 ]
#> ..$ BEGAIN : num[0 , 0 ]
#> ..$ FBXO21 : num[0 , 0 ]
#> ..$ ERP29 : num[0 , 0 ]
#> ..$ VPS13A : num[0 , 0 ]
#> ..$ CYBA : num[0 , 0 ]
#> ..$ S1PR3 : num[0 , 0 ]
#> ..$ RDH14 : num[0 , 0 ]
#> ..$ DIS3L2 : num[0 , 0 ]
#> ..$ RNF185 : num[0 , 0 ]
#> ..$ SLC50A1 : num[0 , 0 ]
#> ..$ UPF3B : num[0 , 0 ]
#> ..$ ART1 : num[0 , 0 ]
#> ..$ FAHD2B : num[0 , 0 ]
#> ..$ ZNF684 : num[0 , 0 ]
#> ..$ PTGER4 : num[0 , 0 ]
#> ..$ DMPK : num[0 , 0 ]
#> ..$ ID3 : num[0 , 0 ]
#> ..$ STIM1 : num[0 , 0 ]
#> ..$ PLRG1 : num[0 , 0 ]
#> ..$ MPC2 : num[0 , 0 ]
#> ..$ ENY2 : num[0 , 0 ]
#> ..$ FUCA1 : num[0 , 0 ]
#> ..$ TCAF1 : num[0 , 0 ]
#> ..$ MTSS1 : num[0 , 0 ]
#> ..$ STAU2 : num[0 , 0 ]
#> ..$ MIS12 : num[0 , 0 ]
#> ..$ ZBTB3 : num[0 , 0 ]
#> ..$ GLRX3 : num[0 , 0 ]
#> ..$ PLEKHB2 : num[0 , 0 ]
#> ..$ BAD : num[0 , 0 ]
#> ..$ HIPK1 : num[0 , 0 ]
#> ..$ ATP6V1D : num[0 , 0 ]
#> ..$ TMEM140 : num[0 , 0 ]
#> ..$ WRNIP1 : num[0 , 0 ]
#> ..$ OXCT2 : num[0 , 0 ]
#> ..$ PAQR5 : num[0 , 0 ]
#> ..$ DHDDS : num[0 , 0 ]
#> ..$ ARSL : num[0 , 0 ]
#> ..$ HSD3B2 : num[0 , 0 ]
#> ..$ GAS8 : num[0 , 0 ]
#> ..$ CLUAP1 : num[0 , 0 ]
#> ..$ INO80D : num[0 , 0 ]
#> ..$ ROMO1 : num[0 , 0 ]
#> ..$ ZNF785 : num[0 , 0 ]
#> ..$ AFM : num[0 , 0 ]
#> ..$ FAM180A : num[0 , 0 ]
#> ..$ MAP3K20 : num[0 , 0 ]
#> ..$ ATF6B : num[0 , 0 ]
#> ..$ EYA3 : num[0 , 0 ]
#> ..$ BAIAP2 : num[0 , 0 ]
#> ..$ USP8 : num[0 , 0 ]
#> ..$ C1RL : num[0 , 0 ]
#> ..$ CERS5 : num[0 , 0 ]
#> ..$ ZIC3 : num[0 , 0 ]
#> ..$ WDR83OS : num[0 , 0 ]
#> ..$ CHRNB1 : num[0 , 0 ]
#> ..$ TSEN15 : num[0 , 0 ]
#> ..$ ERP27 : num[0 , 0 ]
#> ..$ RIT1 : num[0 , 0 ]
#> ..$ PSMG3 : num[0 , 0 ]
#> ..$ GIPC2 : num[0 , 0 ]
#> ..$ LMBRD1 : num[0 , 0 ]
#> ..$ KCNJ15 : num[0 , 0 ]
#> ..$ CRTC1 : num[0 , 0 ]
#> ..$ ZNF300 : num[0 , 0 ]
#> ..$ EP300 : num[0 , 0 ]
#> ..$ PAOX : num[0 , 0 ]
#> ..$ PRPF6 : num[0 , 0 ]
#> ..$ YES1 : num[0 , 0 ]
#> ..$ FAM107B : num[0 , 0 ]
#> ..$ AFAP1L1 : num[0 , 0 ]
#> ..$ TNRC18 : num[0 , 0 ]
#> ..$ MYO18A : num[0 , 0 ]
#> ..$ TOPORS : num[0 , 0 ]
#> ..$ LATS2 : num[0 , 0 ]
#> ..$ METRN : num[0 , 0 ]
#> ..$ ADRB1 : num[0 , 0 ]
#> ..$ PCGF2 : num[0 , 0 ]
#> ..$ JCAD : num[0 , 0 ]
#> ..$ THNSL2 : num[0 , 0 ]
#> .. [list output truncated]
#> $ Im :List of 18439
#> ..$ A2M : num [1:4, 1:4] 21.1303 8.5032 -0.1702 0.0115 8.5032 ...
#> ..$ NAT2 : num [1:4, 1:4] 13.976574 8.887219 -0.000738 0.204521 8.887219 ...
#> ..$ ACADM : num [1:4, 1:4] 29.2401 19.2994 -0.0804 0.0802 19.2994 ...
#> ..$ ACADS : num [1:4, 1:4] 80.4057 58.1105 -0.0033 0.3098 58.1105 ...
#> ..$ ACAT1 : num [1:4, 1:4] 20.1191 12.8542 -0.0879 0.0348 12.8542 ...
#> ..$ ACVRL1 : num [1:4, 1:4] 50.7264 17.4294 -0.0244 0.2293 17.4294 ...
#> ..$ PSEN1 : num [1:4, 1:4] 7.97e+01 5.61e+01 5.33e-05 3.96e-01 5.61e+01 ...
#> ..$ ADA : num [1:4, 1:4] 13.372464 8.412482 -0.000587 0.179358 8.412482 ...
#> ..$ SGCA : num [1:4, 1:4] 23.2053 14.9894 -0.0154 0.1293 14.9894 ...
#> ..$ ADRB2 : num [1:4, 1:4] 2.67e+01 1.73e+01 1.11e-05 2.80e-01 1.73e+01 ...
#> ..$ ADRB3 : num [1:4, 1:4] 3.89 2.24 7.35e-06 1.23e-01 2.24 ...
#> ..$ AGL : num [1:4, 1:4] 3.21e+01 2.39e+01 4.28e-05 2.94e-01 2.39e+01 ...
#> ..$ AGXT : num [1:4, 1:4] 13.2306 8.44311 -0.00958 0.16776 8.44311 ...
#> ..$ ABCD1 : num [1:4, 1:4] 2.31e+01 1.12e+01 2.63e-04 2.47e-01 1.12e+01 ...
#> ..$ ALDOB : num [1:4, 1:4] 4.21742 2.36829 -0.89891 0.00155 2.36829 ...
#> ..$ AMPD1 : num [1:4, 1:4] 1.34e+01 8.33 8.22e-05 2.11e-01 8.33 ...
#> ..$ APC : num [1:4, 1:4] 3.60e+01 1.83e+01 -6.30e-06 3.10e-01 1.83e+01 ...
#> ..$ APOA1 : num [1:4, 1:4] 29.737555 17.393144 -0.000214 0.294107 17.393144 ...
#> ..$ APOC3 : num [1:4, 1:4] 44.97451 28.34827 -0.00416 0.30889 28.34827 ...
#> ..$ APOE : num [1:4, 1:4] 8.6294 4.8548 -0.2787 0.0135 4.8548 ...
#> ..$ APOH : num [1:4, 1:4] 3.48 2.33 9.79e-06 1.22e-01 2.33 ...
#> ..$ ARG1 : num [1:4, 1:4] 2.22e+01 1.24e+01 3.96e-05 2.59e-01 1.24e+01 ...
#> ..$ ARSB : num [1:4, 1:4] 55.08493 32.15418 -0.00107 0.35459 32.15418 ...
#> ..$ ASL : num [1:4, 1:4] 110.25378 82.71718 -0.00753 0.39028 82.71718 ...
#> ..$ ASPA : num [1:4, 1:4] 39.45687 31.0848 -0.00809 0.27177 31.0848 ...
#> ..$ ASS1 : num [1:4, 1:4] 20.8068 11.2059 -0.1244 0.0192 11.2059 ...
#> ..$ ATM : num [1:4, 1:4] 8.79e+01 5.15e+01 -7.71e-06 4.10e-01 5.15e+01 ...
#> ..$ ATP7A : num [1:4, 1:4] 3.03e+01 1.78e+01 -9.99e-05 3.01e-01 1.78e+01 ...
#> ..$ ATP7B : num [1:4, 1:4] 18.9938 12.1616 -0.0223 0.1923 12.1616 ...
#> ..$ AVPR2 : num [1:4, 1:4] 5.3372 3.8251 0.0119 0.0785 3.8251 ...
#> ..$ BLM : num [1:4, 1:4] 11.0912 6.8233 -0.0199 0.153 6.8233 ...
#> ..$ BRCA2 : num [1:4, 1:4] 12.86985 8.58233 -0.00286 0.16091 8.58233 ...
#> ..$ BTD : num [1:4, 1:4] 6.12e+01 3.71e+01 -4.94e-05 3.76e-01 3.71e+01 ...
#> ..$ BTK : num [1:4, 1:4] 9.01 5.79 2.38e-05 1.80e-01 5.79 ...
#> ..$ SERPING1 : num [1:4, 1:4] 49.4303 25.9894 -0.0531 0.026 25.9894 ...
#> ..$ C2 : num [1:4, 1:4] 5.1053 2.8807 -0.0282 0.0938 2.8807 ...
#> ..$ C3 : num [1:4, 1:4] 5.5626 4.09169 -0.4753 0.00522 4.09169 ...
#> ..$ C6 : num [1:4, 1:4] 1.71e+01 1.04e+01 -1.87e-05 2.39e-01 1.04e+01 ...
#> ..$ C8B : num [1:4, 1:4] 10.09948 6.10957 -0.00615 0.15522 6.10957 ...
#> ..$ CA2 : num [1:4, 1:4] 13.8934 9.074 -0.1683 0.0134 9.074 ...
#> ..$ CACNA1S : num [1:4, 1:4] 1.36e+01 7.37 -2.42e-05 2.17e-01 7.37 ...
#> ..$ CBS : num [1:4, 1:4] 53.8122 32.9343 -0.0282 0.2443 32.9343 ...
#> ..$ CD36 : num [1:4, 1:4] 2.01e+01 9.88 -7.53e-05 2.51e-01 9.88 ...
#> ..$ CD3G : num [1:4, 1:4] 8.57 5.47 8.15e-07 1.79e-01 5.47 ...
#> ..$ CD40LG : num [1:4, 1:4] 7.45 5.07 -8.81e-06 1.70e-01 5.07 ...
#> ..$ CDK4 : num [1:4, 1:4] 1.13e+02 7.64e+01 -2.24e-04 4.43e-01 7.64e+01 ...
#> ..$ CDKN1C : num [1:4, 1:4] 20.0781 5.4559 -0.1557 0.0181 5.4559 ...
#> ..$ CETP : num [1:4, 1:4] 7.99 5.21 -5.77e-05 1.73e-01 5.21 ...
#> ..$ CHRNE : num [1:4, 1:4] 20.9623 14.4304 -0.0016 0.2291 14.4304 ...
#> ..$ LYST : num [1:4, 1:4] 7.79e+01 4.83e+01 -7.13e-06 3.96e-01 4.83e+01 ...
#> ..$ ERCC8 : num [1:4, 1:4] 2.59e+01 2.00e+01 -3.06e-05 2.76e-01 2.00e+01 ...
#> ..$ CLCN1 : num [1:4, 1:4] 6.86468 4.755111 0.000854 0.145611 4.755111 ...
#> ..$ CLCN5 : num [1:4, 1:4] 12.3042 8.8935 -0.0995 0.1002 8.8935 ...
#> ..$ CLN3 : num [1:4, 1:4] 7.92e+01 5.70e+01 9.41e-05 3.92e-01 5.70e+01 ...
#> ..$ CNGA1 : num [1:4, 1:4] 2.21e+01 1.39e+01 6.66e-05 2.56e-01 1.39e+01 ...
#> ..$ COL1A1 : num [1:4, 1:4] 14.2016 9.3211 -0.0846 0.0986 9.3211 ...
#> ..$ COL1A2 : num [1:4, 1:4] 15.9234 8.5407 -0.1707 0.0259 8.5407 ...
#> ..$ COL4A3 : num [1:4, 1:4] 44.8727 14.5751 -0.0866 0.0889 14.5751 ...
#> ..$ COL5A1 : num [1:4, 1:4] 9.4344 5.697 -0.1002 0.0801 5.697 ...
#> ..$ COMP : num [1:4, 1:4] 2.56e+01 1.57e+01 -5.74e-05 2.82e-01 1.57e+01 ...
#> ..$ CP : num [1:4, 1:4] 15.15634 9.35931 0.00319 0.1581 9.35931 ...
#> ..$ CPOX : num [1:4, 1:4] 4.513219 3.030791 0.000223 0.127448 3.030791 ...
#> ..$ CPT2 : num [1:4, 1:4] 46.5891 33.1228 -0.0215 0.24 33.1228 ...
#> ..$ CST3 : num [1:4, 1:4] 126.713 71.5443 -0.0789 0.1423 71.5443 ...
#> ..$ CSTB : num [1:4, 1:4] 219.079 130.367 0.039 0.23 130.367 ...
#> ..$ CYP17A1 : num [1:4, 1:4] 9.6723 5.8395 -0.1018 0.0993 5.8395 ...
#> ..$ CYP19A1 : num [1:4, 1:4] 1.05e+01 6.49 2.92e-05 1.92e-01 6.49 ...
#> ..$ CYP1B1 : num [1:4, 1:4] 17.768 8.54 -0.1127 0.0577 8.54 ...
#> ..$ CYP2D6 : num [1:4, 1:4] 24.6712 14.2669 -0.0171 0.2062 14.2669 ...
#> ..$ DDB2 : num [1:4, 1:4] 1.79e+01 1.12e+01 -9.89e-05 2.48e-01 1.12e+01 ...
#> ..$ DLD : num [1:4, 1:4] 31.1109 21.3943 -0.0193 0.1057 21.3943 ...
#> ..$ SLC26A3 : num [1:4, 1:4] 6.27 4.29 -2.19e-05 1.58e-01 4.29 ...
#> ..$ SLC26A2 : num [1:4, 1:4] 2.58e+01 1.44e+01 3.20e-05 2.76e-01 1.44e+01 ...
#> ..$ TOR1A : num [1:4, 1:4] 24.58101 14.65125 0.00944 0.22075 14.65125 ...
#> ..$ EDN3 : num [1:4, 1:4] 6.2312 4.4511 -0.0122 0.1259 4.4511 ...
#> ..$ EMD : num [1:4, 1:4] 1.53e+02 1.02e+02 3.06e-04 4.58e-01 1.02e+02 ...
#> ..$ EPB42 : num [1:4, 1:4] 4.98e+01 2.70e+01 -4.45e-05 3.49e-01 2.70e+01 ...
#> ..$ EPHX1 : num [1:4, 1:4] 89.4385 50.2844 -0.0287 0.2761 50.2844 ...
#> ..$ EPOR : num [1:4, 1:4] 4.73e+01 3.83e+01 1.07e-04 3.27e-01 3.83e+01 ...
#> ..$ ERCC5 : num [1:4, 1:4] 55.2361 31.7118 -0.0106 0.2411 31.7118 ...
#> ..$ ESR1 : num [1:4, 1:4] 2.89e+01 1.75e+01 -9.72e-05 2.93e-01 1.75e+01 ...
#> ..$ EXT1 : num [1:4, 1:4] 6.06e+01 3.97e+01 -2.75e-06 3.72e-01 3.97e+01 ...
#> ..$ F11 : num [1:4, 1:4] 9.2722 6.9846 -0.0479 0.1282 6.9846 ...
#> ..$ F13A1 : num [1:4, 1:4] 1.97e+01 1.52e+01 -5.97e-05 2.49e-01 1.52e+01 ...
#> ..$ F5 : num [1:4, 1:4] 8.9738 2.5771 -0.0193 0.0899 2.5771 ...
#> ..$ F9 : num [1:4, 1:4] 5.30e-02 -1.27e-02 -8.84e-07 1.97e-02 -1.27e-02 ...
#> ..$ FABP2 : num [1:4, 1:4] 1.6058 1.040901 0.000356 0.081299 1.040901 ...
#> ..$ FANCA : num [1:4, 1:4] 1.24e+01 8.58 -1.56e-05 2.10e-01 8.58 ...
#> ..$ FAH : num [1:4, 1:4] 33.6237 20.5502 -0.0247 0.1941 20.5502 ...
#> ..$ MS4A2 : num [1:4, 1:4] 2.36 1.51 -6.48e-06 1.02e-01 1.51 ...
#> ..$ FECH : num [1:4, 1:4] 19.8438 14.0706 -0.109 0.0737 14.0706 ...
#> ..$ FGFR2 : num [1:4, 1:4] 2.73e+01 2.22e+01 9.89e-05 2.70e-01 2.22e+01 ...
#> ..$ FH : num [1:4, 1:4] 30.8169 20.7292 -0.0877 0.1206 20.7292 ...
#> ..$ FSHR : num [1:4, 1:4] 5.00 3.48 -1.62e-05 1.42e-01 3.48 ...
#> ..$ FUCA1 : num [1:4, 1:4] 47.0674 35.6609 0.0474 0.1502 35.6609 ...
#> ..$ FUT1 : num [1:4, 1:4] 1.05e+01 5.59 2.34e-05 1.90e-01 5.59 ...
#> ..$ FUT6 : num [1:4, 1:4] 9.4824 6.308 -0.1181 0.0734 6.308 ...
#> ..$ G6PC : num [1:4, 1:4] 4.9649 3.043 -0.161 0.0424 3.043 ...
#> ..$ GALC : num [1:4, 1:4] 47.60447 23.17732 -0.00838 0.2488 23.17732 ...
#> .. [list output truncated]
#> $ conv0 : num [1:1500, 1] 0 0 0 0 0 0 0 0 0 0 ...
#> ..- attr(*, "names")= chr [1:1500] "DDT" "FOSL1" "FKBP11" "PIM1" ...
#> $ conv : num [1:18439, 1] 0 0 0 0 0 0 0 0 0 0 ...
#> ..- attr(*, "names")= chr [1:18439] "A2M" "NAT2" "ACADM" "ACADS" ...
#> $ features_high: chr [1:1500] "DDT" "FOSL1" "FKBP11" "PIM1" ...
#> $ features_all : chr [1:18439] "A2M" "NAT2" "ACADM" "ACADS" ...