0

我观察了十个不同地点和不同年龄的鱼(同一物种)的大小。对于我在 1、2、3、4 岁等时知道的每个人的大小 -> 我对某个年龄的一个人有一个值(没有针对特定年龄的多个测量值!)。[根据数据,我知道至少有一个人口在所有年龄组中都变大了 -> 可能是一种不依赖于其他固定因素的人口效应]。

> 'data.frame': 688 obs. of  24 variables: 
>$ FishID       : int  99191 99192 99193 99203 99206 99207 99208 99213 99215 99217 ...  
>$ Sex: Factor w/ 3 levels "F","M","U": 1 2 1 2 2 1 2 2 2 2 ...  
>$ Elevation: int  829 829 829 829 829 829 829 829 829 829 ...  
>$ Mean_Sed: num  14992 14992 14992 14992 14992 ...  
>$ Age : int  1 1 1 1 1 1 1 1 1 1 ...  
>$ Length: num  113 111 117 106 111 ...  
>$ Location     : Factor w/ 10 levels " 1 Tavanasa",..: 1 1 1 1 1 1 1 1 1
>$ Catchment_km2: num  598 598 598 598 598 ...      
>$ Slope : num  1.08 1.08 1.08 1.08 1.08 ...      
>$ No_Trout: int  36 36 36 36 36 36 36 36 36 36 ...      
>$ Qt: Factor w/ 197 Levels "1005109","1011605",..: 122 122 122 122 122 122 122 122 122 122 ...

为简化起见,我在年龄组内进行比较。

当然,鱼不仅地点不同,独立种群也不同。我尝试了这种嵌套方法:Size~(1|Site/FishID)+.... 但得到一个错误:每个分组因子的级别数必须小于观察数。

我认为,错误的原因是我没有一个单独的 FishID 的多个鱼大小读数。显然,在我的情况下,不能使用这种方法/嵌套设计(School-Teacher-Scores)。

=> 我如何最好地解决这个问题?是否有另一种方法可以告诉 lmer 该位置内的鱼的大小可能更相似,因为它们属于同一种群。

提前谢谢了

4

0 回答 0