数据结构如下：
year pattern genotype y1 y2 y3 y4 y5 y6 y7 y8 y9 y10
year 有3个水平（1,2,3）
pattern有2个水平（1,2）
genotype有2个水平（1,2）
y1-y10代表一年中的十次观测值

要求：Correlation analyses are quite naïve. More sophisticated analysis, such as linear mixed effect model using year as a random factor, should be performed to reduce the number of separate analyses and type I error.
请各位帮忙指导下程序应该怎么编写啊，谢谢！