==============================================================================
A1. PRESENT-NATURAL MAMMALS: logit(IUCN threatened) ~ log10(body mass)
==============================================================================
N species with LC/NT/VU/EN/CR (terrestrial): 3700
  threatened (VU/EN/CR): 1001  non-threatened (LC/NT): 2699
mass range (g): 1.6 - 4400000.0
Pooled terrestrial: intercept = -2.267 +/- 0.096
  slope per log10 g = 0.5011 +/- 0.0331  z = 15.13  p = 0
  OR per 10x mass = exp(beta) = 1.651  AIC = 4084.0  N = 3700
Stepwise (4 mass bins): AIC = 4105.9  (continuous AIC = 4084.0, dAIC = +21.9; smaller = better)
  bin threshold log10g: ['1.57', '2.13', '3.18'] -> ['38', '136', '1,500'] g
  bin coefficients: 0.355, 0.790, 1.490 (vs bin1)
------------------------------------------------------------------------------
Per-region threatened slopes (terrestrial):
  NorthAmerica   slope = +0.1921 +/- 0.1045  z = +1.84  N = 557
  SouthAmerica   slope = +0.6978 +/- 0.1022  z = +6.83  N = 696
  Sahul          slope = +0.6078 +/- 0.1388  z = +4.38  N = 325
  Eurasia        slope = +0.7457 +/- 0.0598  z = +12.47  N = 1129
  Africa         slope = +0.1814 +/- 0.0796  z = +2.28  N = 539

==============================================================================
A2. LATE-QUATERNARY EXTINCT MAMMALS: logit(extinct) ~ log10(body mass)
==============================================================================
N terrestrial species in 5 land regions: 3928  (extinct EP+EX+EW: 222)
  NorthAmerica   N= 633  extinct=  43  extant=590
  SouthAmerica   N= 924  extinct=  59  extant=865
  Sahul          N= 413  extinct=  62  extant=351
  Eurasia        N=1345  extinct=  47  extant=1298
  Africa         N= 613  extinct=  11  extant=602
statuses excluded from this pool: DD, and region=Islands/Marine/Ocean/Unassigned (see report)
Pooled 5-region: intercept = -8.297 +/- 0.354
  slope per log10 g = 1.5196 +/- 0.0778  z = 19.54  p = 0
  OR per 10x mass = 4.570  AIC = 999.7
Stepwise (4 mass bins): AIC = 1297.2  (continuous AIC = 999.7, dAIC = +297.5)
  bin coefficients: 1.511, 2.269, 4.780 (vs bin1)
------------------------------------------------------------------------------
Joint model: logit(extinct) ~ region + log10mass + region:log10mass
  intercept            -11.6914 +/- 1.1702  z=-9.99  p=0.0000
  log10mass            +2.0127 +/- 0.2285  z=+8.81  p=0.0000
  Americas             +0.8860 +/- 1.4100  z=+0.63  p=0.5298
  Sahul                +6.0857 +/- 1.3055  z=+4.66  p=0.0000
  Africa               +1.5274 +/- 2.0984  z=+0.73  p=0.4667
  log10mass:Americas   +0.2121 +/- 0.2882  z=+0.74  p=0.4619
  log10mass:Sahul      -0.7894 +/- 0.2771  z=-2.85  p=0.0044
  log10mass:Africa     -0.4324 +/- 0.4150  z=-1.04  p=0.2974

Per-region extinction slopes from joint model:
  Americas   slope = +2.2247
  Sahul      slope = +1.2233
  Eurasia    slope = +2.0127 +/- 0.2285  (SE from cov)
  Africa     slope = +1.5802
    (Americas: slope = +2.2247 +/- 0.1757, z = +12.66, p = 0.0000)
    (Sahul: slope = +1.2233 +/- 0.1567, z = +7.81, p = 0.0000)
    (Africa: slope = +1.5802 +/- 0.3464, z = +4.56, p = 0.0000)

Slope-difference tests vs Eurasia (climate reference):
  slope(Americas) - slope(Eurasia) = +0.2121 +/- 0.2882  z = +0.74  p = 0.4619
  slope(Sahul) - slope(Eurasia) = -0.7894 +/- 0.2771  z = -2.85  p = 0.0044
  slope(Africa) - slope(Eurasia) = -0.4324 +/- 0.4150  z = -1.04  p = 0.2974
------------------------------------------------------------------------------
Hunting vs climate: logit(extinct) ~ hunt + log10mass + log10mass:hunt
  intercept            -11.6914 +/- 1.1702  z=-9.99  p=0.0000
  log10mass            +2.0127 +/- 0.2285  z=+8.81  p=0.0000
  hunt                 +3.3709 +/- 1.2497  z=+2.70  p=0.0070
  log10mass:hunt       -0.2668 +/- 0.2521  z=-1.06  p=0.2899
  N(hunting Americas+Sahul) = 1970, N(climate Eurasia) = 1345, extinct: 164 vs 47

Figure saved: /Volumes/Storage/DeepSeek/workspace/desktop-publish-2026-10-30/tests/12-extinction-size/figures/size_selectivity.png