我看看十个不同地点和不同年龄的鱼(同一物种)的大小。对于我知道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 ...
为了简化我在年龄组内的比较。
当然,鱼不仅在网站方面有所不同,在独立人群方面也有所不同。我尝试了这种嵌套方法:大小〜(1 | Site / FishID)+ .... 但是得到一个错误:每个分组因子的级别数必须是<观察次数。
我认为,错误的原因在于我没有一个单独的FishID的鱼尺寸读数。 显然,在我的情况下,这个apporach /嵌套设计(School-Teacher-Scores)不能使用。
=>我如何才能最好地解决这个问题?还有另一种方法告诉lmer,该地区内的鱼可能更大,因为它们属于同一种群。
非常感谢提前