求助：混合模型的双因素方差分析

在做产量丰产性分析时，想用到混合模型的双因素方差分析，其中sites为试验点，为random因子；variety为品种，固定因子；rep为每个sites内的3次重复；yield为产量。求助：在R里面如何实现，用aov函数的话如何编写公式？谢谢！

sites	variety	yield	rep
94	粳籼89	11.25	1
94	盐稻3号	11.4	1
94	中籼糯909	11.18	1
94	澳青占	9.85	1
94	七秀占3号	10.6	1
94	青珍8号	10.05	1
94	豫籼3号	11.18	1
94	特三矮2号(ck)	12.6	1
100	粳籼89	10.82	1
100	盐稻3号	12.08	1
100	中籼糯909	10.82	1
100	澳青占	9.76	1
100	七秀占3号	11.4	1
100	青珍8号	11.1	1
100	豫籼3号	11.4	1
100	特三矮2号(ck)	12.84	1

不懂，也想学习