#################################
#DNE calculated on whole teeth
> summary(DNE_mod, WAIC=T)
 Family: categorical 
  Links: muFL = logit; muIN = logit 
Formula: diet ~ mesh_08 + (1 | Tooth) + (1 | genus) 
   Data: data (Number of observations: 2668) 
Samples: 4 chains, each with iter = 3000; warmup = 1000; thin = 1;
         total post-warmup samples = 8000

Group-Level Effects: 
~genus (Number of levels: 22) 
                   Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sd(muFL_Intercept)     5.06      0.56     4.05     6.22 1.00     8013     6135
sd(muIN_Intercept)     4.14      0.55     3.14     5.28 1.00     7973     5981

~Tooth (Number of levels: 223) 
                   Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sd(muFL_Intercept)     0.37      0.28     0.02     1.04 1.00     7305     3849
sd(muIN_Intercept)     0.45      0.35     0.02     1.29 1.00     6109     4202

Population-Level Effects: 
               Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
muFL_Intercept    -2.55      1.23    -4.99    -0.19 1.00     3252     4800
muIN_Intercept    -3.96      1.20    -6.37    -1.67 1.00     4220     5633
muFL_mesh_08       1.15      0.82    -0.46     2.80 1.00     8887     6750
muIN_mesh_08       3.99      1.19     1.84     6.49 1.00     6275     5669


#################################
#DNE calculated on segments
> summary(seg_DNE_mod, WAIC=T)
 Family: categorical 
  Links: muFL = logit; muIN = logit 
Formula: diet ~ ep_08 + (1 | Tooth) + (ep_08 | segment) + (1 | taxon) + (1 | genus) 
   Data: data (Number of observations: 870) 
Samples: 4 chains, each with iter = 3000; warmup = 1000; thin = 1;
         total post-warmup samples = 8000

Group-Level Effects: 
~genus (Number of levels: 20) 
                   Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sd(muFL_Intercept)     0.81      0.63     0.03     2.35 1.00     5524     3223
sd(muIN_Intercept)     0.84      0.62     0.03     2.34 1.00     5829     2706

~segment (Number of levels: 15) 
                               Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sd(muFL_Intercept)                 0.62      0.48     0.02     1.78 1.00     6314     3311
sd(muFL_ep_08)                     0.59      0.46     0.02     1.70 1.00     6448     4004
sd(muIN_Intercept)                 0.60      0.47     0.03     1.76 1.00     6100     3461
sd(muIN_ep_08)                     0.60      0.45     0.02     1.69 1.00     7242     3649
cor(muFL_Intercept,muFL_ep_08)     0.00      0.57    -0.95     0.94 1.00     8521     4961
cor(muIN_Intercept,muIN_ep_08)    -0.00      0.58    -0.95     0.95 1.00    10428     5150

~taxon (Number of levels: 20) 
                   Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sd(muFL_Intercept)     2.49      0.57     1.50     3.71 1.00     7032     5861
sd(muIN_Intercept)     2.59      0.58     1.58     3.85 1.00     6122     5562

~Tooth (Number of levels: 58) 
                   Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
sd(muFL_Intercept)     0.67      0.51     0.02     1.92 1.00     7335     3709
sd(muIN_Intercept)     0.65      0.50     0.03     1.85 1.00     5297     3179

Population-Level Effects: 
               Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
muFL_Intercept    -1.05      2.89    -6.68     4.61 1.00     7795     6104
muIN_Intercept    -0.96      2.87    -6.62     4.61 1.00     6995     6037
muFL_ep_08         0.10      1.14    -2.19     2.37 1.00     6847     5330
muIN_ep_08         0.38      1.16    -1.85     2.73 1.00     6874     5328