rs <- lmer(y ~ time + (time| id), data=dataset, REML=FALSE)

rs
Linear mixed model fit by maximum likelihood  ['lmerMod']
Formula: y ~ time + (time | id)
   Data: dataset
      AIC       BIC    logLik  deviance  df.resid
  66.2413   70.0756   -27.1206  54.2413         8
Random effects:
 Groups   Name       Variance  Std.Dev.   Corr
 id      (Intercept) 8.846024. 2.97423
          time       0.006426  0.08016   -1.00
 Residual            1.473614. 1.21393
Number of obs: 14, groups:  id, 4

Fixed Effects:
              Estimate Std.Error t value
(Intercept)   99.04580  1.57192   63.01
time          -0.16086  0.04457  -3.61