Raw counts

Table 1: Data table of artifact class by site

Figure S1: Stacked bar plots of artifact class in systematic collections Figure S2: Stacked bar plots of artifact class in piece plot collections

Ratios

Figure S3: Artifact density Figure S4: Ratio of Levallois products to Nubian Levallois cores Figure S5: Ratio of retouched tools to cores Figure S6: Ratio of debitage to cores

Summary for all sites

Weighted patination scores for each site

# A tibble: 15 × 3
   Site   total_count weighted_score
   <chr>        <int>          <dbl>
 1 TH.069        1833           NA  
 2 TH.076          89            5.1
 3 TH.123          57            4.2
 4 TH.143          97            4.7
 5 TH.187          54            4.7
 6 TH.188          33            4.4
 7 TH.191         929           NA  
 8 TH.205         189            3.5
 9 TH.268         662            2.9
10 TH.323         356            5.4
11 TH.363         202            3.4
12 TH.418         868            3.2
13 TH.419         287            3.1
14 TH.571          88           NA  
15 TH.591         797           NA  

Figure 2: Histograms of patination scores per site Figure S7: Stacked bar plot core type % by patination score per site Figure S8: Boxplots of geo mean organized by patina score

Cluster Analysis

Figure 3: Violin plot of patination scores by site with three-centroid cluster analysis Kruskal-Wallis and Dunn’s Tests

    Kruskal-Wallis rank sum test

data:  Patination by Cluster_3
Kruskal-Wallis chi-squared = 1991, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 1990.962, df = 2, p-value = 0

                           Comparison of x by group                            
                                (No adjustment)                                
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   4.067598
         |    0.0000*
         |
     III |   33.52105   34.35357
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

$Adjusted_P_Value
numeric(0)

Violin plot of bi cores versus Nub cores by patina

Raw counts

Data table of core types by site

Data table of core types by cluster

Table 2: Data table of core types by site and cluster

Figure S9: Stacked bar plot of total core counts per site Figure S10: Stacked bar plot % of total core types per cluster

Reduction strategy summary

Figure 5: Ternary plot of main reduction strategies by site Chi Square Test of differences in main reduction strategies

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 1401.1, df = 42, p-value < 2.2e-16

Fisher’s Exact Test

  Cluster1 Cluster2   P_Value Adjusted_P_Value
1        I       II 9.999e-05       0.00029997
2        I      III 9.999e-04       0.00299970
3       II      III 9.999e-05       0.00029997

Descriptive Statistics

Table S1: Nub core descriptive statistics by site

Table S2: Nub core descriptive statistics by cluster

Morphology

Figure S11: MDR angle box plots by site Figure 7: MDR angle box plots by cluster MDR boxplot tests

    Kruskal-Wallis rank sum test

data:  Nub_MDR_Angle by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 10.918, df = 2, p-value = 0.004258

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 10.9181, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |  -0.975660
         |     0.4938
         |
     III |  -2.218928  -2.959322
         |     0.0397    0.0046*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S12: MDR prep pie charts by cluster Chi Square Test of Nub core MDR prep between clusters

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 3.6126, df = 4, p-value = 0.461

Figure S13: Ternary plot of MDR prep by cluster Figure S14: Nubian platform prep pie charts for all sites Figure S15: Nubian platform prep pie charts by cluster Figure S16: Nubian core shape pie charts for all sites Figure S17: Nubian core shape pie charts by cluster Chi Square Test of Nub core shape between clusters

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 15.537, df = 6, p-value = 0.01647

Fisher’s Exact Test

# A tibble: 3 × 4
  Cluster_1 Cluster_2 Fisher_p_value p_value_adjusted
  <chr>     <chr>              <dbl>            <dbl>
1 I         II               0.576            1      
2 I         III              0.0750           0.225  
3 II        III              0.00150          0.00450

Figure S18: Nubian core final scar for all sites Figure S19: Nubian core final scar pie charts by cluster Chi Square Test of Nub core final scar between clusters

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 3.8127, df = 3, p-value = 0.2824

Metric Analysis

Figure S20: Boxplot of Nub core length by site Figure S21: Boxplot of Nub core length by cluster Kruskal Wallis Nub Core Length Test

    Kruskal-Wallis rank sum test

data:  Max_Thick by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 235.01, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 180.9646, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   3.702671
         |    0.0003*
         |
     III |   9.689317   11.93788
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S22: Boxplot of Nub core width by site Figure S23: Boxplot of Nub core width by cluster Kruskal Wallis Nub Core Width Test

    Kruskal-Wallis rank sum test

data:  Max_Width by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 157.25, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 157.2549, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   2.573247
         |    0.0151*
         |
     III |   8.357723   11.52055
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S24: Boxplot of Nubian core thickness by site Figure S25: Boxplot of Nub core thickness by cluster Kruskal Wallis Nub Core Thick Test

    Kruskal-Wallis rank sum test

data:  Max_Thick by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 239.23, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 239.2257, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   1.326070
         |     0.2772
         |
     III |   8.847278   14.83982
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S26: Boxplot of Nubian core Index of Flattening by site Figure S27: Boxplot of Nubian core Index of Flattening by cluster Kruskal Wallis Test of Nub Core IF

    Kruskal-Wallis rank sum test

data:  Index_F by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 42.007, df = 2, p-value = 7.558e-10

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 42.0065, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |  -1.101186
         |     0.4062
         |
     III |   2.236916   6.478489
         |     0.0379    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S28: Boxplot of Nubian core Index of Elongation by site Figure S29: Boxplot of Nubian core Index of Elongation by cluster Kruskal Wallis and Dunn’s Tests of Nub Core IE

    Kruskal-Wallis rank sum test

data:  Index_E by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 16.605, df = 2, p-value = 0.0002479

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 16.6048, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   2.121828
         |     0.0508
         |
     III |   3.580412   2.987359
         |    0.0005*    0.0042*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S30: Nub core L x W scatterplot by cluster Figure 8: Boxplot of Nubian core geometric mean by site Figure S31: Boxplot of Nubian core geometric mean by cluster Kruskal Wallis & Dunns Tests of Nub core geo mean

    Kruskal-Wallis rank sum test

data:  Geo_Mean by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 240.74, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 240.7385, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   2.699195
         |    0.0104*
         |
     III |   10.00168   14.43400
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Descriptive Statistics

Table S3: Bi core descriptive statistics by site

Table S4: Bi core descriptive statistics by cluster

Morphology

Figure S32: Bi core dorsal scar pie charts for all sites Figure 11: Bi core dorsal scar pie charts by cluster Chi Square Test of Bi core dorsal scars between clusters

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 36.011, df = 6, p-value = 2.743e-06

Fisher’s Exact Test

# A tibble: 3 × 4
  Cluster_1 Cluster_2 Fisher_p_value p_value_adjusted
  <chr>     <chr>              <dbl>            <dbl>
1 I         II              0.518             1      
2 I         III             0.225             0.676  
3 II        III             0.000500          0.00150

Figure S33: Bi core platform prep pie charts for all sites Figure S34: Bi core platform prep pie charts by cluster Chi Square Test of Platform differences between clusters

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 7.5627, df = 4, p-value = 0.109

Figure S35: Bi core shape pie charts for all sites Figure S36: Bi core shape pie charts by cluster Chi Square test of shape differences between clusters

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 38.562, df = 8, p-value = 5.926e-06

Fisher’s Exact Test

# A tibble: 3 × 4
  Cluster_1 Cluster_2 Fisher_p_value p_value_adjusted
  <chr>     <chr>              <dbl>            <dbl>
1 I         II              0.0951           0.0951  
2 I         III             0.00198          0.00198 
3 II        III             0.000310         0.000310

Figure S37: Bi core final scar pie charts for all sites Figure S38: Bi core final scar pie charts by cluster Chi Square Test of bi core final scar between clusters

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 3.8127, df = 3, p-value = 0.2824

Metric Analysis

Figure S39: Boxplot of bi core length by site Figure S40: Boxplot of Nub core length by cluster Kruskal Wallis and Dunn’s Test of bi core length

    Kruskal-Wallis rank sum test

data:  Max_Thick by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 94.684, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 133.6285, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   0.238863
         |     1.0000
         |
     III |   5.267489   10.78207
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S41: Boxplot of bi core width by site Figure S42: Boxplot of bi core width by cluster Kruskal Wallis and Dunn’s tests of bi core width

    Kruskal-Wallis rank sum test

data:  Max_Width by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 126.32, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 126.3235, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   0.121488
         |     1.0000
         |
     III |   5.022640   10.52265
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S43: Boxplot of bi core thickness by site Figure S44: Boxplot of bi core thickness by cluster Kruskal Wallis and Dunn’s tests of bi core thickness

    Kruskal-Wallis rank sum test

data:  Max_Thick by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 94.684, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 94.6845, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |  -0.802867
         |     0.6331
         |
     III |   3.517468   9.389299
         |    0.0007*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S45: Boxplot of bi core Index of Flattening by site Figure S46: Boxplot of bi core Index of Flattening by cluster Kruskal Wallis and Dunn’s tests of bi core IF

    Kruskal-Wallis rank sum test

data:  Index_F by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 13.031, df = 2, p-value = 0.00148

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 13.0309, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |  -1.091507
         |     0.4126
         |
     III |  -2.485955  -2.861845
         |    0.0194*    0.0063*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S47: Boxplot of bi core Index of Elongation by site Figure S48: Boxplot of bi core Index of Elongation by cluster Kruskal Wallis test of bi core IE

    Kruskal-Wallis rank sum test

data:  Index_E by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 4.123, df = 2, p-value = 0.1273

Figure S49: bi core L x W scatterplot by cluster Figure 12: Boxplot of bi core geometric mean by site Figure S50: Boxplot of bi core geometric mean by cluster Kruskal Wallis & Dunn’s Tests of bi core geo mean

    Kruskal-Wallis rank sum test

data:  Geo_Mean by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 151.8, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 151.7967, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   0.037183
         |     1.0000
         |
     III |   5.419765   11.56829
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Descriptive Statistics

Table S5: Uni core descriptive statistics by site

Table S6: Uni core descriptive statistics by cluster

Morphology

Figure S51: Uni core dorsal scar pie charts for all sites Figure 14: Uni core dorsal scar pie charts by cluster Chi Square Test of uni core dorsal scars between clusters

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 31.274, df = 8, p-value = 0.0001256

Fisher’s Exact Test

# A tibble: 3 × 4
  Cluster_1 Cluster_2 Fisher_p_value p_value_adjusted
  <chr>     <chr>              <dbl>            <dbl>
1 I         II              0.389             1      
2 I         III             0.00850           0.0255 
3 II        III             0.000500          0.00150

Figure S52: Uni core platform prep pie charts for all sites Figure S53: Uni core platform prep pie charts by cluster Chi Square Test of Platform differences between clusters

    cortical&flat dihedral faceted
I              22        0       1
II             29        2       2
III           269        6      42

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 5.7461, df = 4, p-value = 0.2189

Figure S54: Uni core shape pie charts for all sites Figure 15: uni core shape pie charts by cluster Chi Square Test of shape differences between clusters

    cortical&flat dihedral faceted
I              22        0       1
II             29        2       2
III           269        6      42

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 49.266, df = 8, p-value = 5.654e-08

Fisher’s Exact Test

# A tibble: 3 × 4
  Cluster_1 Cluster_2 Fisher_p_value p_value_adjusted
  <chr>     <chr>              <dbl>            <dbl>
1 I         II            0.216             0.648    
2 I         III           0.0000200         0.0000600
3 II        III           0.00001000        0.0000300

Metric Analysis

Figure S55: Boxplot of uni core length by site Figure S56: Boxplot of uni core length by cluster Kruskal Wallis and Dunn’s Test of uni core length

    Kruskal-Wallis rank sum test

data:  Max_Length by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 101.21, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 101.2057, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   0.550643
         |     0.8728
         |
     III |   8.026905   6.841831
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S57: Boxplot of uni core width by site Figure S58: Boxplot of uni core width by cluster Kruskal Wallis and Dunn’s tests of uni core width

    Kruskal-Wallis rank sum test

data:  Max_Width by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 45.989, df = 2, p-value = 1.032e-10

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 45.9893, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |  -1.387422
         |     0.2480
         |
     III |   4.172433   5.737139
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S59: Boxplot of uni core thickness by site Figure S60: Boxplot of uni core thickness by cluster Kruskal Wallis and Dunn’s tests of uni core thickness

    Kruskal-Wallis rank sum test

data:  Max_Thick by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 95.855, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 95.855, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   0.936877
         |     0.5232
         |
     III |   8.046092   6.356117
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S61: Boxplot of uni core Index of Flattening by site Figure S62: Boxplot of uni core Index of Flattening by cluster Kruskal Wallis and Dunn’s tests of uni core IF

    Kruskal-Wallis rank sum test

data:  Index_F by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 15.682, df = 2, p-value = 0.0003932

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 15.6823, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   1.981956
         |     0.0712
         |
     III |   3.897608   1.087643
         |    0.0001*     0.4151

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S63: Boxplot of uni core Index of Elongation by site Figure S64: Boxplot of uni core Index of Elongation by cluster Kruskal Wallis test of uni core IE

    Kruskal-Wallis rank sum test

data:  Index_E by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 9.8925, df = 2, p-value = 0.00711

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 9.8925, df = 2, p-value = 0.01

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   2.496398
         |    0.0188*
         |
     III |   3.072735  -0.358490
         |    0.0032*     1.0000

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S65: uni core L x W scatterplot by cluster Figure 16: Boxplot of uni core geometric mean by site Figure S66: Boxplot of uni core geometric mean by cluster Kruskal Wallis & Dunn’s Tests of uni core geo mean

    Kruskal-Wallis rank sum test

data:  Geo_Mean by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 160.86, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 160.8607, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   0.207413
         |     1.0000
         |
     III |   9.818731   8.976391
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Metric Analysis

Figure S67: debitage L x W scatterplot by site Figure S68: debitage L x W scatterplot by cluster Kruskal Wallis test of debitage IE

    Kruskal-Wallis rank sum test

data:  Index_E by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 63.416, df = 1, p-value = 1.674e-15

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 63.4157, df = 1, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I
---------+-----------
     III |   7.963397
         |    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Descriptive Statistics

Table S7: Levallois point descriptive statistics by site

Table S8: Levallois point descriptive statistics by cluster

Metric Analysis

Figure S69: Boxplot of Levallois point length by site Figure S70: Boxplot of Levallois point length by cluster Kruskal Wallis and Dunn’s test of Levallois point length

    Kruskal-Wallis rank sum test

data:  Max_Length by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 13.348, df = 2, p-value = 0.001263

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 13.3479, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   0.300930
         |     1.0000
         |
     III |   2.584102   3.407263
         |    0.0146*    0.0010*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S71: Boxplot of Levallois point width by site Figure S72: Boxplot of Levallois point width by cluster Kruskal Wallis test of Levallois point width

    Kruskal-Wallis rank sum test

data:  Max_Width by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 3.0931, df = 2, p-value = 0.213

Figure S73: Boxplot of Levallois point thickness by site Figure S74: Boxplot of Levallois point thickness by cluster Kruskal Wallis test of Levallois point thickness

    Kruskal-Wallis rank sum test

data:  Max_Thick by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 2.1718, df = 2, p-value = 0.3376

Figure S75: Boxplot of Levallois point Index of Flattening by site Figure S76: Boxplot of Levallois point Index of Flattening by cluster Kruskal Wallis test of Levallois point IF

    Kruskal-Wallis rank sum test

data:  Index_F by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 3.6636, df = 2, p-value = 0.1601

Figure S77: Boxplot of Levallois point Index of Elongation by site Figure S78: Boxplot of Levallois point Index of Elongation by cluster Kruskal Wallis and Dunn’s Tests of Levallois point IE

    Kruskal-Wallis rank sum test

data:  Index_E by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 8.0129, df = 2, p-value = 0.0182

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 8.0129, df = 2, p-value = 0.02

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |  -0.345512
         |     1.0000
         |
     III |   1.553295   2.792537
         |     0.1805    0.0078*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S79: Levallois point L x W scatterplot by cluster Figure S80: Boxplot of Levallois point geometric mean by site Figure 18: Boxplot of Levallois point geometric mean by cluster Kruskal Wallis & Dunn’s tests of Levallois point geo mean

    Kruskal-Wallis rank sum test

data:  Geo_Mean by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 9.6275, df = 2, p-value = 0.008117

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 9.6275, df = 2, p-value = 0.01

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   0.315790
         |     1.0000
         |
     III |   2.236957   2.871348
         |     0.0379    0.0061*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S81: Boxplot of Levallois point TCSA by site Figure S82: Boxplot of Levallois point TCSA by cluster Kruskal Wallis test of Levallois point TCSA

    Kruskal-Wallis rank sum test

data:  TCSA by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 2.8092, df = 2, p-value = 0.2455

Figure S83: Boxplot of Levallois point TCSP by site Figure S84: Boxplot of Levallois point TCSP by cluster Kruskal Wallis test of Levallois point TCSP

    Kruskal-Wallis rank sum test

data:  TCSP by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 2.9824, df = 2, p-value = 0.2251

Figure S85: Boxplot of preferential scar length by site Figure S86: Boxplot of preferential scar length by cluster Kruskal Wallis and Dunn’s test of preferential scar length

    Kruskal-Wallis rank sum test

data:  Lev_Scar_Length by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 78.904, df = 2, p-value < 2.2e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 78.9043, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   2.106848
         |     0.0527
         |
     III |   5.618736   8.202080
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S87: Boxplot of preferential scar width by site Figure S88: Boxplot of preferential scar width by cluster Kruskal Wallis test of preferential scar width

    Kruskal-Wallis rank sum test

data:  Lev_Scar_Width by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 49.313, df = 2, p-value = 1.958e-11

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 49.3133, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |  -0.252625
         |     1.0000
         |
     III |   2.791355   6.973859
         |    0.0079*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S89: Boxplot of preferential scar Index of Elongation by site Figure S90: Boxplot of preferential scar Index of Elongation by cluster Kruskal Wallis and Dunn’s Tests of preferential scar IE

    Kruskal-Wallis rank sum test

data:  Lev_Scar_Index_E by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 7.6151, df = 2, p-value = 0.0222

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 7.6151, df = 2, p-value = 0.02

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   2.756268
         |    0.0088*
         |
     III |   2.545415  -0.291302
         |    0.0164*     1.0000

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S91: Preferential scar L x W scatterplot by cluster Figure S92: Boxplot of preferential scar size by site Figure S93: Boxplot of preferential scar size by cluster Kruskal Wallis and Dunn’s Tests of preferential scar size

    Kruskal-Wallis rank sum test

data:  Lev_Scar_Size by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 71.46, df = 2, p-value = 3.038e-16

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 71.4602, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   0.860724
         |     0.5841
         |
     III |   4.399454   8.176696
         |    0.0000*    0.0000*

alpha = 0.05
Reject Ho if p <= alpha/2

Levallois points versus final scars

Figure S94: Boxplots of point vs. scar length by cluster Kruskal Wallis test of Levallois point versus scar length

[[1]]

    Kruskal-Wallis rank sum test

data:  Value by Variable
Kruskal-Wallis chi-squared = 0.66516, df = 1, p-value = 0.4147


[[2]]

    Kruskal-Wallis rank sum test

data:  Value by Variable
Kruskal-Wallis chi-squared = 22.295, df = 1, p-value = 2.338e-06


[[3]]

    Kruskal-Wallis rank sum test

data:  Value by Variable
Kruskal-Wallis chi-squared = 17.166, df = 1, p-value = 3.426e-05

Figure S95: Boxplots of point vs. scar width by cluster Kruskal Wallis tests of point versus scar width

[[1]]

    Kruskal-Wallis rank sum test

data:  Value by Variable
Kruskal-Wallis chi-squared = 5.9864, df = 1, p-value = 0.01442


[[2]]

    Kruskal-Wallis rank sum test

data:  Value by Variable
Kruskal-Wallis chi-squared = 39.006, df = 1, p-value = 4.224e-10


[[3]]

    Kruskal-Wallis rank sum test

data:  Value by Variable
Kruskal-Wallis chi-squared = 61.615, df = 1, p-value = 4.176e-15

Figure 20: Boxplots of discarded point vs. scar size by cluster Kruskal Wallis tests of point versus scar size

[[1]]

    Kruskal-Wallis rank sum test

data:  Value by Variable
Kruskal-Wallis chi-squared = 2.474, df = 1, p-value = 0.1157


[[2]]

    Kruskal-Wallis rank sum test

data:  Value by Variable
Kruskal-Wallis chi-squared = 43.256, df = 1, p-value = 4.802e-11


[[3]]

    Kruskal-Wallis rank sum test

data:  Value by Variable
Kruskal-Wallis chi-squared = 45.101, df = 1, p-value = 1.871e-11

Figure 20b: Boxplots of point vs. scar size by cluster w/o TH.69 Figure 20c: Boxplots of point vs. scar size at TH.69

Types

Table 4: Data table of restricted tool types

Table S9: Data table of unrestricted tool types

Figure S96: Bar plots of retouched tool counts per site (no Levallois products) Figure S97: Bar plot of retouched tool counts per cluster

Blank Selection

Figure S98: Bar plots of retouched tool blank selection count per site Figure 19: Pie charts of retouched tool blank selection by cluster Chi Square Test of differences between clusters

    blade cortical flake Lev_product
I      10        1     5           6
II      6        8    17          15
III    10       20    29           7

    Pearson's Chi-squared test

data:  contingency_table
X-squared = 23.678, df = 6, p-value = 0.0005986

Fisher’s Exact Test

# A tibble: 3 × 4
  Cluster_1 Cluster_2 Fisher_p_value p_value_adjusted
  <chr>     <chr>              <dbl>            <dbl>
1 I         II              0.0299           0.0299  
2 I         III             0.000724         0.000724
3 II        III             0.0335           0.0335  

Metric Analysis

Figure S99: Boxplot of retouched tool length by site Figure S100: Boxplot of retouched tool length by cluster Kruskal Wallis test of retouched tool length

    Kruskal-Wallis rank sum test

data:  Max_Thick by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 0.42104, df = 2, p-value = 0.8102

Figure S101: Boxplot of retouched tool width by site Figure S102: Boxplot of retouched tool width by cluster Kruskal Wallis test of retouched tool width

    Kruskal-Wallis rank sum test

data:  Max_Width by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 10.221, df = 2, p-value = 0.006034

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 10.2208, df = 2, p-value = 0.01

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   0.386974
         |     1.0000
         |
     III |   2.503890   2.655513
         |    0.0184*    0.0119*

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S103: Boxplot of retouched tool thickness by site Figure S104: Boxplot of retouched tool thickness by cluster Kruskal Wallis test of retouched tool thickness

    Kruskal-Wallis rank sum test

data:  Max_Thick by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 0.42143, df = 2, p-value = 0.81

Figure S105: Boxplot of retouched tool Index of Flattening by site Figure S106: Boxplot of retouched tool Index of Flattening by cluster Kruskal Wallis test of retouched tool IF

    Kruskal-Wallis rank sum test

data:  Index_F by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 5.1134, df = 2, p-value = 0.07756

Figure S107: Boxplot of retouched tool Index of Elongation by site Figure S108: Boxplot of retouched tool Index of Elongation by cluster Kruskal Wallis test of uni core IE

    Kruskal-Wallis rank sum test

data:  Index_E by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 8.8316, df = 2, p-value = 0.01208

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 8.8316, df = 2, p-value = 0.01

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   2.200984
         |     0.0416
         |
     III |   2.965897   0.700145
         |    0.0045*     0.7258

alpha = 0.05
Reject Ho if p <= alpha/2

Figure S109: retouched tool L x W scatterplot by cluster Figure S110: Boxplot of retouched tool geometric mean by site Figure S111: Boxplot of retouched tool geometric mean by cluster Kruskal Wallis & Dunn’s Tests of retouched tool geo mean

    Kruskal-Wallis rank sum test

data:  Geo_Mean by as.factor(Cluster_3)
Kruskal-Wallis chi-squared = 16.802, df = 2, p-value = 0.0002246

  Kruskal-Wallis rank sum test

data: x and group
Kruskal-Wallis chi-squared = 16.8024, df = 2, p-value = 0

                           Comparison of x by group                            
                                 (Bonferroni)                                  
Col Mean-|
Row Mean |          I         II
---------+----------------------
      II |   1.779196
         |     0.1128
         |
     III |   3.888477   2.391189
         |    0.0002*     0.0252

alpha = 0.05
Reject Ho if p <= alpha/2