Mixed Model Approach in Archaeology

Dear Stackexchange Community, I need some advice regarding a statistical analysis I want to conduct as part of my PhD. The general question of this is whether we can detect statistically significant (and relevant) differences in vessel morphology (shape) and morphometry (dimensions) of specific types (i got ten) between specific sites (five in total). Note, these types are similar across sites, as they follow an inter-regionally established vessel shape or form, produced by—supposedly—many producers across a wider geographic area. However, based on several theoretical concepts, it is assumed that each producer created their own respective variant due to different traditions of learning as well as environmental and social influences. So the idea behind this is to verify production patterns. Basically, the question is whether we can differentiate between places of production based on the morphology and morphometry of standardized vessels. The problem is manifold, and the main problem is the quantity and availability of the data, which is not good. However, there is nothing I can do about this, as this is simply the way it is due to excavation/publication bias, etc. I have one main site with the majority of the data, and I now wanted to, exploratively, use the available data to analyze tendencies. I have attached a crosstab which shows my data and the respective availability problems. Type Site 1 Site 2 Site 3 Site 4 Site 5 Type 1 73 28 32 0 3 Type 2 92 60 2 0 7 Type 3 74 78 2…