在lmer中连接两个因素(人口和站点)而没有多次测量?

时间:2017-06-01 10:29:51

标签: r nested lme4 replicate population

我看看十个不同地点和不同年龄的鱼(同一物种)的大小。对于我知道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,该地区内的鱼可能更大,因为它们属于同一种群。

非常感谢提前

0 个答案:

没有答案