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
###########################################