Bonjour,
Je travaille sur la densité de diatomées en fonction du temps (date) et du traitement. J'aimerai démontrer que l'on a une différence de densité selon la date et selon le traitement
Pour cela, j'utilise le modèle linéaire à effet fixe.
ou EU: traitement
DW: densité
date=date
> fm=lme(DW~EU+date,random=~1|EU)
> summary(fm)
je ne comprends la sortie de R (notament les NA), comment dois je conclure???
Linear mixed-effects model fit by REML
Data: NULL
AIC BIC logLik
224.3678 233.1622 -106.1839
Random effects:
Formula: ~1 | EU
(Intercept) Residual
StdDev: 4.241099 5.509349
Fixed effects: DW ~ EU + date
Value Std.Error DF t-value p-value
(Intercept) 3.496153 4.788676 32 0.730088 0.4706
EUEU2 4.144664 6.405674 0 0.647030 NaN
EUEU3 -8.630701 6.405674 0 -1.347353 NaN
date 4.449244 0.478171 32 9.304718 0.0000
Correlation:
(Intr) EUEU2 EUEU3
EUEU2 -0.669
EUEU3 -0.669 0.500
date -0.325 0.000 0.000
Standardized Within-Group Residuals:
Min Q1 Med Q3 Max
-2.3430222 -0.5752665 0.1968327 0.5759900 2.0478112
Number of Observations: 36
Number of Groups: 3
Warning message:
production de NaN in: pt(q, df, lower.tail, log.p)
quelles sont les tests à faire pour valider le modèle?
Merci
Seb