Dear R-list.
How would I plot model result over actual data?
Here are the lme results:
Linear mixed-effects model fit by REML
Data: D-set
AIC BIC logLik
884.3326 911.428 -432.1663
Random effects:
Formula: ~C | ID
Structure: General positive-definite, Log-Cholesky parametrization
StdDev Corr
(Intercept) 13.7293823 (Intr)
C 0.7425786 -0.544
Residual 6.7688776
Fixed effects: Y ~ A * C + B * C
Value Std.Error DF t-value p-value
(Intercept) 29.662389 7.979584 98 3.717285 0.0003
A 0.185470 0.102510 13 1.809294 0.0936
C -0.741605 0.692078 98 -1.071563 0.2865
B 0.006185 0.002265 98 2.730804 0.0075
A:C -0.003325 0.008243 98 -0.403405 0.6875
C:B -0.000679 0.000278 98 -2.444981 0.0163
Correlation:
(Intr) A C B A:C
A -0.865
C -0.540 0.437
B -0.138 -0.009 0.104
A:C 0.473 -0.559 -0.770 0.030
C:B 0.152 0.013 -0.218 -0.958 -0.042
Standardized Within-Group Residuals:
Min Q1 Med Q3 Max
-3.55799126 -0.53364258 0.01324915 0.52255745 2.68609309
Number of Observations: 117
Number of Groups: 15
########################################### Kare Edvardsen <kare.edvardsen at nilu.no> Norwegian Institute for Air Research (NILU) Polarmiljosenteret NO-9296 Tromso http://www.nilu.no Swb. +47 77 75 03 75 Dir. +47 77 75 03 90 Fax. +47 77 75 03 76 Mob. +47 90 74 60 69 ###########################################