Wildlife on stamps reveal taxonomic bias in biodiversity representativeness in China

Code and data for "Wildlife on stamps reveal taxonomic bias in biodiversity representativeness in China" This repository contains three R scripts developed by Wendong Xie for analyzing taxonomic representation bias on Chinese wildlife stamps. Script 1: Chi_square_tests_and_correlation_of_taxonomic_bias.R With data inside the script. This script performs chi-square tests comparing stamp representation of taxonomic groups against three biodiversity baselines: global described species, China's described species, and IUCN-assessed species. Outputs standardized residual bar plots (Figure S2) and scatter plots with Pearson correlations between stamp and IUCN residuals (Figure S4). Script 2: GLM_of_repeated_species_issuance.R With data in Ne_of_species_on_stamp.csv. This script analyzes factors influencing the number of stamp issuance events per species using Poisson GLM selection, model averaging, and robustness checks. Produces Table 1 (model-averaged coefficients) and Table 2 (top candidate models) in the manuscript. Script 3: annual_trend_analysis_of_taxonomic_bias.R With data in stamp_grouped_by_year.csv. This script conducts year-by-year chi-square tests comparing stamp counts per taxonomic group against global and China described species proportions. Computes Pearson residuals and generates data for Figure 2c and Figure S3, tracking temporal changes in representation bias.