线性混合效应模型教程

shrimp.lm.sex.tank <- lm(M2BW~SexID+TankID,shrimp)

shrimp.lm.sex.tank.sum <- summary(shrimp.lm.sex.tank)

shrimp.lm.sex.tank.sum

shrimp[,mean(M2BW,na.rm=TRUE),by=.(SexID,TankID)]

>shrimp.lm.sum

Call:

lm(formula = M2BW ~ SexID + TankID, data = shrimp)

Residuals:
     Min       1Q   Median       3Q      Max 
-21.4916  -3.0478  -0.1697   2.8166  20.4084 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  35.8560     0.1276  281.10   <2e-16 ***
SexIDMale    -6.2082     0.1493  -41.58   <2e-16 ***
TankIDT2     -2.9644     0.1492  -19.87   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1<