[R-sig-ME] lme4 at flexLambda Error

[Steve Walker]

I'm sorry but I couldn't reproduce the error.  Could you provide 
sessionInfo()?  Mine is below.

Steve

 > sessionInfo()
R version 3.1.1 (2014-07-10)
Platform: x86_64-apple-darwin13.1.0 (64-bit)

locale:
[1] C

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] lme4_2.0      Rcpp_0.11.2.1 Matrix_1.1-4  setup_0.0-1

loaded via a namespace (and not attached):
[1] MASS_7.3-33     compiler_3.1.1  grid_3.1.1      lattice_0.20-29
[5] minqa_1.2.3     nlme_3.1-117    nloptr_1.0.0    splines_3.1.1
[9] tools_3.1.1
 >



On 2015-02-01 5:03 PM, bbonit at tin.it wrote:
> set.seed(2352)
> nsubj <- 20
> nreps <- 2
> length <- 20
> n <- nsubj*nreps*length
> fr <- expand.grid(
>      time=1:length,
>      id=factor(1:nsubj),
>      rep=factor(1:nreps))
>
> rho <- .8
> sd.b <- 3
> b <- sd.b * as.vector(replicate(nsubj, arima.sim(n=length, model=list(ar=c(rho)))))
> sd.eps <- 1
> beta <- .5
> fr <- within(fr,{
>      mu <- beta*time
>      y <- mu + model.matrix(~0+as.factor(time):id)%*%b +
>          rnorm(n, sd=sd.eps)
> })
>
> m <-  flexLmer(y ~ time + ar1d(~(time|id)), data=fr)