Plant–insect interactions from the mid-Cretaceous at Puy-Puy (Aquitaine Basin, western France) indicates preferential herbivory for angiosperms amid a forest of ferns, gymnosperms, and angiosperms

The nine in-text figures and table below (Appendices S1–S10), and the additional text and excel files attached, provide the raw data, summaries of the raw data, rarefaction analyses, and nonmetric multidimensional scale analyses (NMDS) that support the discussions of the main text. The raw data and their summaries of provide for each plant species or morphotype values important for assessment of their herbivory: percentage of specimens herbivorized, damage type (DT) richness, DT frequency, DT host-plant specificity, herbivorized surface area as a proportion of total surface area, and feeding event occurrences. The rarefaction analyses furnished evaluations of whether the number of samples was sufficient, given the surface area covered by those samples. For comparison, the number of samples was rarified to the number of DTs in those samples. Lastly, two NMDS analyses produced the relationships between the plant orders present in the plant assemblage and their interactive functional feeding groups (FFGs). A separate NMDS analysis shows the association between the three most herbivorized species and their FFGs.