Hi
?predict
over
?order #(ed)
data
could be one way.
HTH
Petr
On 14 Dec 2005 at 13:36, Kare Edvardsen wrote:
Date sent: Wed, 14 Dec 2005 13:36:42 +0100
From: Kare Edvardsen <ked at nilu.no>
To: R-help <r-help at stat.math.ethz.ch>
Subject: [R] Plotting model over data
> 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
> ###########################################
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
Petr Pikal
petr.pikal at precheza.cz