C'est mon premier petit projet en R... Et j'ai deja des soucis...
J'ai des données sur 31 individus qui sont mesurées 7 fois dans le temps. J'ai ajusté un modèle avec la fonction "lme" à mes données, et quand je trace un graphique avec les valeurs ajustées par rapport a ma réponse, ça a l'air juste. Je cherche à retrouver la fonction correspondant a mon ajustement, et pour ça j'ai demandé un "summary" de mon modèle :
Code : Tout sélectionner
Linear mixed-effects model fit by maximum likelihood
Data: glucose2
AIC BIC logLik
-763.5905 -687.423 400.7952
Random effects:
Formula: ~TimePoint | SubjectID
Structure: General positive-definite, Log-Cholesky parametrization
StdDev Corr
(Intercept) 0.085272828 (Intr)
TimePoint 0.001960009 -0.694
Formula: ~TimePoint | Period %in% SubjectID
Structure: General positive-definite, Log-Cholesky parametrization
StdDev Corr
(Intercept) 0.116220940 (Intr)
TimePoint 0.001457469 -0.843
Residual 0.066005142
Fixed effects: log(BloodGlucose) ~ poly(TimePoint, 2) * Grade
Value Std.Error DF t-value p-value
(Intercept) 5.879948 0.0318297 341 184.73177 0.0000
poly(TimePoint, 2)1 -0.753267 0.4282114 341 -1.75910 0.0795
poly(TimePoint, 2)2 -3.949686 0.1361768 341 -29.00410 0.0000
GradeDiabGrade 1 0.260718 0.0445524 27 5.85194 0.0000
GradeDiabGrade 2 0.253268 0.0445524 27 5.68473 0.0000
GradeDiabGrade 3 0.252461 0.0456082 27 5.53544 0.0000
poly(TimePoint, 2)1:GradeDiabGrade 1 3.165915 0.6015212 341 5.26318 0.0000
poly(TimePoint, 2)2:GradeDiabGrade 1 1.368242 0.1895393 341 7.21878 0.0000
poly(TimePoint, 2)1:GradeDiabGrade 2 4.095052 0.6015212 341 6.80783 0.0000
poly(TimePoint, 2)2:GradeDiabGrade 2 1.292958 0.1895393 341 6.82158 0.0000
poly(TimePoint, 2)1:GradeDiabGrade 3 4.121715 0.6179009 341 6.67051 0.0000
poly(TimePoint, 2)2:GradeDiabGrade 3 1.472031 0.1927865 341 7.63555 0.0000
Je voudrais tracer la courbe correspondant a mon modèle, donc :
y(t)= [b0 + b01 grade_1 + b02 grade_2 + b03 grade_3]
+ [b1 + b11 grade_1 + b12 grade_2 + b13 grade_3] t
+ [b2 + b21 grade_1 + b22 grade_2 + b23 grade_3] t²
avec grade_i =0 si l'individu appartient au groupe i, 1 sinon,
b0 <- 5.889098
b1 <- -0.790891
b2 <- -3.949542
b01 <- 0.258151
b02 <- 0.243972
b03 <- 0.243311
b11 <- 3.168266
b21 <- 1.368098
b12 <- 4.130165
b22 <- 1.292814
b13 <- 4.159338
b23 <- 1.471887
je trouve quelque chose qui n'a rien a voir avec mes données. Apparement, les échelles ne sont pas bonnes du tout (il faudrait diviser par environ 2 pour l'ordonnée et environ 60 pour l'abscisse..)...
Est-ce que quelqu'un comprend mon problème et si oui, pourrait m'aider ?
merci
Gaëlle