Model comparison - subcategory meta-regression

[ThisIsLorenzo]

We want to compare our current model (subgroup correlated effects model):

mod_chem <- rma.mv(ln_TMF ~ PFAS_ID - 1,
                   V = VCV_chem,
                   random = list(~ PFAS_ID | FW_ID), 
                   struct = "DIAG",
                   test = "t",
                   sparse = T,
                   data = TMF_data3,
                   #verbose = TRUE, 
                   control=list(iter.max=1000, 
                                rel.tol=1e-8))

summary(mod_chem)
Multivariate Meta-Analysis Model (k = 987; method: REML)

    logLik    Deviance         AIC         BIC        AICc   
-1892.9472   3785.8944   4069.8944   4754.3367   4122.4325  

Variance Components:

outer factor: FW_ID   (nlvls = 115)
inner factor: PFAS_ID (nlvls = 71)

             estim    sqrt  k.lvl  fixed          level 
tau^2.1     0.0079  0.0887      1     no  10:2 Cl-PFESA 
tau^2.2     0.0000  0.0000      2     no      10:2 FTSA 
tau^2.3     0.0078  0.0885      1     no      10:3 FTCA 
tau^2.4     0.0079  0.0886      1     no        3m-PFOS 
tau^2.5     0.0078  0.0885      1     no       4:2 FTSA 
tau^2.6     0.0078  0.0885      1     no        4m-PFOS 
tau^2.7     0.0078  0.0886      1     no        5m-PFOS 
tau^2.8     1.9333  1.3904      3     no      6:2 diPAP 
tau^2.9     0.9843  0.9921      4     no       6:2 FTSA 
tau^2.10    0.0000  0.0000      2     no    6:2 H-PFESA 
tau^2.11    0.0078  0.0885      1     no       7:3 FTCA 
tau^2.12    0.0079  0.0887      1     no        8:2 FTS 
tau^2.13    0.0000  0.0008      2     no       8:2 FTSA 
tau^2.14    0.0078  0.0885      1     no    8:2 H-PFESA 
tau^2.15    0.0000  0.0000      3     no          ADONA 
tau^2.16    1.9527  1.3974      3     no      br-EtFOSA 
tau^2.17    0.0078  0.0886      1     no        br-FOSA 
tau^2.18    0.0079  0.0887      1     no        br-PFDA 
tau^2.19    0.0000  0.0000      2     no      br-PFDoDA 
tau^2.20    0.0079  0.0887      1     no        br-PFOA 
tau^2.21    0.3085  0.5554     12     no        br-PFOS 
tau^2.22    0.0078  0.0886      1     no      br-PFTeDA 
tau^2.23    0.0078  0.0885      1     no     cis-PFECHS 
tau^2.24    0.0079  0.0888      1     no       Cl-PFESA 
tau^2.25    0.0000  0.0000      2     no        dm-PFOS 
tau^2.26    0.0078  0.0886      1     no        EtFOSAA 
tau^2.27    0.0537  0.2317     37     no          F-53B 
tau^2.28    0.0000  0.0000      2     no           FBSA 
tau^2.29    0.0079  0.0886      1     no          FHpPA 
tau^2.30    0.3754  0.6127     28     no           FOSA 
tau^2.31    0.0000  0.0055      2     no          FOSAA 
tau^2.32    0.6342  0.7964      2     no     H-PFMO2OSA 
tau^2.33    0.0000  0.0002      2     no        HFPO-DA 
tau^2.34    0.8065  0.8981      2     no       HFPO-TeA 
tau^2.35    0.1402  0.3745      3     no       HFPO-TrA 
tau^2.36    0.0078  0.0886      1     no       iso-PFOS 
tau^2.37    1.6174  1.2718      3     no       l-EtFOSA 
tau^2.38    0.3054  0.5526      3     no         l-FOSA 
tau^2.39    0.0078  0.0886      1     no         l-PFDA 
tau^2.40    0.0079  0.0886      1     no       l-PFDoDA 
tau^2.41    0.0078  0.0886      1     no         l-PFOA 
tau^2.42    0.2924  0.5407     20     no         l-PFOS 
tau^2.43    0.0079  0.0887      1     no       l-PFTeDA 
tau^2.44    0.0001  0.0077      3     no        MeFOSAA 
tau^2.45    0.0000  0.0000      2     no        mm-PFOS 
tau^2.46    0.5641  0.7511      6     no       NEtFOSAA 
tau^2.47    0.5011  0.7079      4     no       NMeFOSAA 
tau^2.48    0.0078  0.0886      1     no        PF4OPeA 
tau^2.49    0.0000  0.0008     15     no           PFBA 
tau^2.50    0.3373  0.5807     20     no           PFBS 
tau^2.51    0.3138  0.5602     85     no           PFDA 
tau^2.52    0.1406  0.3750     84     no         PFDoDA 
tau^2.53    0.7195  0.8483     21     no           PFDS 
tau^2.54    0.0000  0.0000      2     no         PFECHS 
tau^2.55    0.0000  0.0000     18     no          PFHpA 
tau^2.56    0.6271  0.7919      9     no          PFHpS 
tau^2.57    1.3652  1.1684     15     no          PFHxA 
tau^2.58    0.0000  0.0000      8     no         PFHxDA 
tau^2.59    0.3283  0.5730     43     no          PFHxS 
tau^2.60    0.0027  0.0523      2     no         PFMOAA 
tau^2.61    0.2921  0.5405     87     no           PFNA 
tau^2.62    0.0079  0.0887      1     no         PFNOBS 
tau^2.63    0.0079  0.0887      1     no           PFNS 
tau^2.64    0.0000  0.0000      2     no        PFO5DoA 
tau^2.65    0.5133  0.7165     71     no           PFOA 
tau^2.66    0.1766  0.4202    110     no           PFOS 
tau^2.67    0.0000  0.0000     10     no          PFPeA 
tau^2.68    0.1914  0.4375     47     no         PFTeDA 
tau^2.69    0.5311  0.7288     65     no         PFTrDA 
tau^2.70    0.2868  0.5355     93     no         PFUnDA 
tau^2.71    0.0078  0.0886      1     no   trans-PFECHS 

Test for Residual Heterogeneity:
QE(df = 916) = 13135.7882, p-val < .0001

Test of Moderators (coefficients 1:71):
F(df1 = 71, df2 = 916) = 18.5142, p-val < .0001

Model Results:

                      estimate      se     tval   df    pval    ci.lb    ci.ub      
PFAS_ID10:2 Cl-PFESA    0.5235  0.7819   0.6695  916  0.5033  -1.0110   2.0580      
PFAS_ID10:2 FTSA        1.2462  0.2490   5.0046  916  <.0001   0.7575   1.7348  *** 
PFAS_ID10:3 FTCA        0.5273  0.8847   0.5960  916  0.5513  -1.2090   2.2636      
PFAS_ID3m-PFOS          1.1100  0.4702   2.3607  916  0.0185   0.1872   2.0328    * 
PFAS_ID4:2 FTSA        -1.1700  1.8633  -0.6279  916  0.5302  -4.8269   2.4869      
PFAS_ID4m-PFOS          1.2500  0.4281   2.9195  916  0.0036   0.4097   2.0903   ** 
PFAS_ID5m-PFOS          1.4300  0.3800   3.7628  916  0.0002   0.6842   2.1758  *** 
PFAS_ID6:2 diPAP       -1.8059  0.9303  -1.9411  916  0.0526  -3.6317   0.0200    . 
PFAS_ID6:2 FTSA        -0.9780  0.6372  -1.5348  916  0.1252  -2.2286   0.2725      
PFAS_ID6:2 H-PFESA      0.9898  0.4307   2.2981  916  0.0218   0.1445   1.8350    * 
PFAS_ID7:3 FTCA        -0.0834  1.3225  -0.0630  916  0.9497  -2.6788   2.5120      
PFAS_ID8:2 FTS          0.5878  0.7297   0.8055  916  0.4207  -0.8443   2.0199      
PFAS_ID8:2 FTSA         0.8844  0.4137   2.1378  916  0.0328   0.0725   1.6964    * 
PFAS_ID8:2 H-PFESA      1.4100  0.4777   2.9519  916  0.0032   0.4726   2.3474   ** 
PFAS_IDADONA           -0.7057  0.3912  -1.8039  916  0.0716  -1.4735   0.0621    . 
PFAS_IDbr-EtFOSA       -0.0774  1.0189  -0.0760  916  0.9395  -2.0771   1.9223      
PFAS_IDbr-FOSA         -0.1863  0.1078  -1.7282  916  0.0843  -0.3979   0.0253    . 
PFAS_IDbr-PFDA          0.5766  0.8632   0.6680  916  0.5043  -1.1174   2.2706      
PFAS_IDbr-PFDoDA        1.0641  0.2368   4.4938  916  <.0001   0.5994   1.5288  *** 
PFAS_IDbr-PFOA          0.7655  0.7641   1.0017  916  0.3167  -0.7342   2.2651      
PFAS_IDbr-PFOS          0.9115  0.1780   5.1201  916  <.0001   0.5621   1.2609  *** 
PFAS_IDbr-PFTeDA        0.3507  1.0502   0.3339  916  0.7385  -1.7105   2.4118      
PFAS_IDcis-PFECHS       0.6530  0.6627   0.9853  916  0.3247  -0.6477   1.9537      
PFAS_IDCl-PFESA         1.2149  0.5586   2.1748  916  0.0299   0.1185   2.3113    * 
PFAS_IDdm-PFOS         -0.2103  0.0892  -2.3587  916  0.0185  -0.3853  -0.0353    * 
PFAS_IDEtFOSAA          0.3684  0.7922   0.4650  916  0.6420  -1.1864   1.9232      
PFAS_IDF-53B            1.1243  0.1242   9.0535  916  <.0001   0.8806   1.3680  *** 
PFAS_IDFBSA             0.2203  0.2060   1.0696  916  0.2851  -0.1839   0.6246      
PFAS_IDFHpPA            1.7600  0.3068   5.7367  916  <.0001   1.1579   2.3621  *** 
PFAS_IDFOSA             0.6388  0.1594   4.0084  916  <.0001   0.3260   0.9515  *** 
PFAS_IDFOSAA           -0.5627  1.3995  -0.4021  916  0.6877  -3.3093   2.1839      
PFAS_IDH-PFMO2OSA       1.7034  0.6281   2.7122  916  0.0068   0.4708   2.9361   ** 
PFAS_IDHFPO-DA         -0.0171  0.7161  -0.0239  916  0.9810  -1.4225   1.3883      
PFAS_IDHFPO-TeA         1.8564  0.6857   2.7076  916  0.0069   0.5108   3.2021   ** 
PFAS_IDHFPO-TrA         1.1827  0.3880   3.0481  916  0.0024   0.4212   1.9443   ** 
PFAS_IDiso-PFOS         1.4400  0.4530   3.1786  916  0.0015   0.5509   2.3291   ** 
PFAS_IDl-EtFOSA         0.0365  0.7557   0.0483  916  0.9615  -1.4465   1.5196      
PFAS_IDl-FOSA           0.1459  0.3423   0.4261  916  0.6702  -0.5260   0.8177      
PFAS_IDl-PFDA           0.9594  0.5504   1.7430  916  0.0817  -0.1209   2.0396    . 
PFAS_IDl-PFDoDA         1.5851  0.3793   4.1786  916  <.0001   0.8407   2.3296  *** 
PFAS_IDl-PFOA           0.4511  1.0247   0.4402  916  0.6599  -1.5600   2.4622      
PFAS_IDl-PFOS           0.8779  0.1444   6.0788  916  <.0001   0.5945   1.1613  *** 
PFAS_IDl-PFTeDA         1.5831  0.3438   4.6040  916  <.0001   0.9083   2.2579  *** 
PFAS_IDMeFOSAA         -0.2257  0.3431  -0.6578  916  0.5108  -0.8990   0.4476      
PFAS_IDmm-PFOS          1.0621  0.1350   7.8662  916  <.0001   0.7971   1.3271  *** 
PFAS_IDNEtFOSAA        -0.4829  0.3717  -1.2990  916  0.1943  -1.2124   0.2467      
PFAS_IDNMeFOSAA         0.1421  0.4029   0.3528  916  0.7243  -0.6485   0.9327      
PFAS_IDPF4OPeA         -1.0788  1.4419  -0.7482  916  0.4545  -3.9086   1.7509      
PFAS_IDPFBA             0.2694  0.2514   1.0715  916  0.2842  -0.2240   0.7627      
PFAS_IDPFBS            -0.0843  0.2778  -0.3036  916  0.7615  -0.6296   0.4609      
PFAS_IDPFDA             1.0288  0.0896  11.4857  916  <.0001   0.8530   1.2045  *** 
PFAS_IDPFDoDA           0.6991  0.0721   9.6994  916  <.0001   0.5576   0.8406  *** 
PFAS_IDPFDS             0.4972  0.2583   1.9246  916  0.0546  -0.0098   1.0042    . 
PFAS_IDPFECHS           0.5311  0.2027   2.6196  916  0.0089   0.1332   0.9289   ** 
PFAS_IDPFHpA           -0.0620  0.1642  -0.3777  916  0.7057  -0.3842   0.2602      
PFAS_IDPFHpS            0.6859  0.3086   2.2224  916  0.0265   0.0802   1.2916    * 
PFAS_IDPFHxA           -0.2902  0.4098  -0.7081  916  0.4791  -1.0944   0.5141      
PFAS_IDPFHxDA           0.1497  0.0837   1.7885  916  0.0740  -0.0146   0.3140    . 
PFAS_IDPFHxS            0.5640  0.1414   3.9883  916  <.0001   0.2865   0.8415  *** 
PFAS_IDPFMOAA          -1.3933  1.1180  -1.2462  916  0.2130  -3.5874   0.8009      
PFAS_IDPFNA             0.7950  0.0915   8.6880  916  <.0001   0.6154   0.9746  *** 
PFAS_IDPFNOBS           0.5953  0.6732   0.8843  916  0.3768  -0.7259   1.9165      
PFAS_IDPFNS             1.3300  0.4060   3.2757  916  0.0011   0.5332   2.1268   ** 
PFAS_IDPFO5DoA          1.7248  0.1695  10.1756  916  <.0001   1.3922   2.0575  *** 
PFAS_IDPFOA             0.2699  0.1417   1.9045  916  0.0572  -0.0082   0.5480    . 
PFAS_IDPFOS             1.1057  0.0693  15.9552  916  <.0001   0.9697   1.2417  *** 
PFAS_IDPFPeA           -0.0404  0.3065  -0.1319  916  0.8951  -0.6420   0.5611      
PFAS_IDPFTeDA           0.3505  0.1062   3.2997  916  0.0010   0.1420   0.5590   ** 
PFAS_IDPFTrDA           0.7151  0.1296   5.5178  916  <.0001   0.4608   0.9694  *** 
PFAS_IDPFUnDA           0.8816  0.0868  10.1522  916  <.0001   0.7112   1.0520  *** 
PFAS_IDtrans-PFECHS     0.8130  0.6170   1.3176  916  0.1880  -0.3979   2.0239      

with a simpler model with only one variance component:

mod_chem2 <- rma.mv(ln_TMF ~ PFAS_ID - 1,
                V = VCV_chem,
                random = list(~ 1|PFAS_ID,
                                        ~ 1|FW_ID), 
                struct = "DIAG",
                est = "t",
                sparse = T,
                data = TMF_data3,
               #verbose = TRUE, 
                   control=list(iter.max=1000, 
                         rel.tol=1e-8))

summary(mod_chem2)
Multivariate Meta-Analysis Model (k = 987; method: REML)

    logLik    Deviance         AIC         BIC        AICc   
-3286.1037   6572.2075   6718.2075   7070.0687   6731.0388   

Variance Components:

            estim    sqrt  nlvls  fixed   factor 
sigma^2.1  0.1432  0.3784     71     no  PFAS_ID 
sigma^2.2  0.3199  0.5656    115     no    FW_ID 

Test for Residual Heterogeneity:
QE(df = 916) = 13135.7882, p-val < .0001

Test of Moderators (coefficients 1:71):
QM(df = 71) = 157.6155, p-val < .0001

Model Results:

                      estimate      se     zval    pval    ci.lb    ci.ub      
PFAS_ID10:2 Cl-PFESA    0.6639  0.8888   0.7469  0.4551  -1.0782   2.4059      
PFAS_ID10:2 FTSA        0.7528  0.4652   1.6182  0.1056  -0.1590   1.6646      
PFAS_ID10:3 FTCA        0.6677  0.9805   0.6809  0.4959  -1.2541   2.5894      
PFAS_ID3m-PFOS          0.4806  0.6147   0.7819  0.4343  -0.7241   1.6853      
PFAS_ID4:2 FTSA        -1.2415  1.9083  -0.6506  0.5153  -4.9818   2.4987      
PFAS_ID4m-PFOS          0.6206  0.5831   1.0642  0.2872  -0.5223   1.7635      
PFAS_ID5m-PFOS          0.8006  0.5488   1.4589  0.1446  -0.2750   1.8761      
PFAS_ID6:2 diPAP       -1.7175  0.6007  -2.8590  0.0042  -2.8949  -0.5401   ** 
PFAS_ID6:2 FTSA        -0.3925  0.3980  -0.9862  0.3240  -1.1727   0.3876      
PFAS_ID6:2 H-PFESA      0.5497  0.5864   0.9373  0.3486  -0.5997   1.6991      
PFAS_ID7:3 FTCA         0.2460  1.4066   0.1749  0.8611  -2.5108   3.0029      
PFAS_ID8:2 FTS          0.7027  0.8457   0.8309  0.4060  -0.9548   2.3602      
PFAS_ID8:2 FTSA         0.3117  0.5740   0.5431  0.5871  -0.8132   1.4367      
PFAS_ID8:2 H-PFESA      0.7806  0.6204   1.2582  0.2083  -0.4354   1.9965      
PFAS_IDADONA           -0.6894  0.5664  -1.2172  0.2235  -1.7995   0.4207      
PFAS_IDbr-EtFOSA        1.5292  0.4129   3.7035  0.0002   0.7199   2.3385  *** 
PFAS_IDbr-FOSA          0.9102  0.3890   2.3399  0.0193   0.1478   1.6726    * 
PFAS_IDbr-PFDA          0.0119  0.9664   0.0123  0.9902  -1.8823   1.9061      
PFAS_IDbr-PFDoDA        0.4202  0.4597   0.9141  0.3607  -0.4808   1.3212      
PFAS_IDbr-PFOA          0.2007  0.8791   0.2283  0.8194  -1.5223   1.9238      
PFAS_IDbr-PFOS          1.3802  0.3844   3.5905  0.0003   0.6268   2.1335  *** 
PFAS_IDbr-PFTeDA       -0.2141  1.1366  -0.1883  0.8506  -2.4418   2.0137      
PFAS_IDcis-PFECHS       0.0236  0.7720   0.0305  0.9756  -1.4895   1.5366      
PFAS_IDCl-PFESA         1.0017  0.7697   1.3014  0.1931  -0.5068   2.5102      
PFAS_IDdm-PFOS         -0.4932  0.4080  -1.2090  0.2267  -1.2928   0.3064      
PFAS_IDEtFOSAA         -0.0147  0.8823  -0.0166  0.9868  -1.7439   1.7146      
PFAS_IDF-53B            0.9808  0.3999   2.4523  0.0142   0.1969   1.7646    * 
PFAS_IDFBSA             0.6614  0.4356   1.5185  0.1289  -0.1923   1.5151      
PFAS_IDFHpPA            1.1306  0.5008   2.2574  0.0240   0.1490   2.1122    * 
PFAS_IDFOSA             1.0168  0.3842   2.6469  0.0081   0.2639   1.7697   ** 
PFAS_IDFOSAA           -0.7705  1.4630  -0.5267  0.5984  -3.6379   2.0968      
PFAS_IDH-PFMO2OSA       1.5693  0.5571   2.8166  0.0049   0.4773   2.6613   ** 
PFAS_IDHFPO-DA          0.0374  0.8230   0.0454  0.9638  -1.5757   1.6505      
PFAS_IDHFPO-TeA         1.8131  0.5429   3.3396  0.0008   0.7490   2.8771  *** 
PFAS_IDHFPO-TrA         1.3382  0.5244   2.5517  0.0107   0.3103   2.3661    * 
PFAS_IDiso-PFOS         0.8106  0.6016   1.3473  0.1779  -0.3686   1.9897      
PFAS_IDl-EtFOSA         1.6430  0.3870   4.2455  <.0001   0.8845   2.4015  *** 
PFAS_IDl-FOSA           0.6603  0.3854   1.7134  0.0866  -0.0950   1.4156    . 
PFAS_IDl-PFDA           0.3946  0.7014   0.5626  0.5737  -0.9800   1.7693      
PFAS_IDl-PFDoDA         1.0204  0.5769   1.7687  0.0770  -0.1104   2.1512    . 
PFAS_IDl-PFOA          -0.1137  1.1131  -0.1021  0.9187  -2.2954   2.0680      
PFAS_IDl-PFOS           1.1387  0.3841   2.9645  0.0030   0.3858   1.8915   ** 
PFAS_IDl-PFTeDA         1.0184  0.5542   1.8374  0.0661  -0.0679   2.1047    . 
PFAS_IDMeFOSAA          0.1489  0.5153   0.2889  0.7726  -0.8610   1.1588      
PFAS_IDmm-PFOS          0.7767  0.4205   1.8471  0.0647  -0.0474   1.6008    . 
PFAS_IDNEtFOSAA        -0.4902  0.4000  -1.2256  0.2203  -1.2742   0.2937      
PFAS_IDNMeFOSAA         0.0580  0.3931   0.1474  0.8828  -0.7126   0.8285      
PFAS_IDPF4OPeA         -1.9793  1.4957  -1.3233  0.1857  -4.9108   0.9523      
PFAS_IDPFBA             0.3631  0.4610   0.7877  0.4309  -0.5404   1.2667      
PFAS_IDPFBS            -0.2227  0.4011  -0.5553  0.5787  -1.0090   0.5635      
PFAS_IDPFDA             0.9753  0.3838   2.5408  0.0111   0.2230   1.7276    * 
PFAS_IDPFDoDA           0.7476  0.3839   1.9474  0.0515  -0.0048   1.5001    . 
PFAS_IDPFDS             0.5971  0.3843   1.5537  0.1203  -0.1561   1.3503      
PFAS_IDPFECHS           0.9741  0.4340   2.2443  0.0248   0.1234   1.8248    * 
PFAS_IDPFHpA           -0.1135  0.4217  -0.2692  0.7878  -0.9401   0.7131      
PFAS_IDPFHpS            0.5516  0.3909   1.4111  0.1582  -0.2146   1.3179      
PFAS_IDPFHxA           -0.0463  0.4430  -0.1044  0.9168  -0.9146   0.8220      
PFAS_IDPFHxDA           0.1648  0.3949   0.4173  0.6765  -0.6092   0.9388      
PFAS_IDPFHxS            0.8260  0.3872   2.1332  0.0329   0.0671   1.5848    * 
PFAS_IDPFMOAA          -1.2846  1.2194  -1.0534  0.2921  -3.6746   1.1055      
PFAS_IDPFNA             0.8635  0.3841   2.2478  0.0246   0.1106   1.6164    * 
PFAS_IDPFNOBS           0.6971  0.7971   0.8746  0.3818  -0.8652   2.2595      
PFAS_IDPFNS             0.7006  0.5670   1.2355  0.2167  -0.4108   1.8120      
PFAS_IDPFO5DoA          1.7888  0.5157   3.4690  0.0005   0.7781   2.7995  *** 
PFAS_IDPFOA             0.6914  0.3845   1.7983  0.0721  -0.0621   1.4450    . 
PFAS_IDPFOS             0.9839  0.3838   2.5635  0.0104   0.2317   1.7362    * 
PFAS_IDPFPeA            0.2105  0.4961   0.4242  0.6714  -0.7619   1.1828      
PFAS_IDPFTeDA           0.7231  0.3848   1.8790  0.0602  -0.0312   1.4774    . 
PFAS_IDPFTrDA           0.5389  0.3840   1.4034  0.1605  -0.2137   1.2915      
PFAS_IDPFUnDA           0.8674  0.3839   2.2597  0.0238   0.1151   1.6197    * 
PFAS_IDtrans-PFECHS     0.1836  0.7331   0.2504  0.8023  -1.2533   1.6204      

The model with the lower AICc value is generally preferred because it indicates a better fit with fewer parameters. In this case: mod_chem has an AICc of 4122.4325. mod_chem2 has an AICc of 6731.0388.

@itchyshin

[itchyshin]

@ThisIsLorenzo

Thanks - actually mod_chem is a much better model, which I am a bit surprised how much better it is

I want you to run a few more models during our meeting so please be ready to run some models for me.

Thanks!

[ThisIsLorenzo]

https://github.com/ThisIsLorenzo/PFAS_Trophic_Magnification/blob/3ad70f6cf9062fbbe36c7ff25cef0a7f187ab71c/R/Analysis_code.Rmd#L898

At this point in the code, I added a section to run and compare simpler models against mod_chembased on their AICc values. If you want to run the models, simply execute the code up to this point @itchyshin