[R-sig-ME] Fwd: Re: R-core post from jw1085 at wildcats.unh.edu requires approval

[Mitchell Maltenfort]

Can you use lme4 instead?

On Mon, May 8, 2017 at 1:35 PM Jochen Wirsing <jw1085 at wildcats.unh.edu>
wrote:

> Dear R Mixed Models SIG,
>
> below you will find an email that I originally sent to the R core group,
> who technically is the maintainer of the nlme package, with which I've
> had some serious issues since the last update. Since this seems to keep
> preventing me from using R for my project, I am reaching out to you,
> hoping this problem can be fixed in a timely manner, as it probably
> affects a lot of other people as well.
>
>
> Thank you very much,
>
> Jochen Wirsing
>
> -------- Forwarded Message --------
> Subject:        Re: R-core post from jw1085 at wildcats.unh.edu requires
> approval
> Date:   Wed, 3 May 2017 15:11:09 +0530
> From:   Deepayan Sarkar <deepayan.sarkar at gmail.com>
> To:     r-core-owner at r-project.org, Jochen Wirsing <
> jw1085 at wildcats.unh.edu>
>
>
>
> Dear Jochen Wirsing,
>
> Although r-core is technically the maintainer of nlme, you should
> write to one of the public R mailing lists get help regarding such
> problems: see
>
> https://www.r-project.org/mail.html
>
> Either R-help or R-SIG-mixed-models should be appropriate.
>
> Best,
> -Deepayan
>
> > ---------- Forwarded message ----------
> > From: Jochen Wirsing <jw1085 at wildcats.unh.edu>
> > To: R-core at r-project.org
> > Cc:
> > Bcc:
> > Date: Mon, 1 May 2017 11:38:09 -0400
> > Subject: Problem with nlme after update to R 3.4.0
> >
> > Dear Sir or Madam,
> >
> > I am a grad student at UNH, teaching myself how to use R. Not having a
> background in programming, it is a bit hard to figure out what's wrong, so
> I thought I should let you know.
> >
> > Last night, I ran an MLM, and it worked just fine. This morning though,
> I saw that there is a new Version of RStudio as well as R itself, so I
> updated both. After doing so, the very identical code that worked last
> night, doesn't work anymore and throws a rather obscure (at least to me)
> error:
> >
> >
> > 15.
> > pdFactor.pdLogChol(X[[i]], ...)
> > 14.
> > FUN(X[[i]], ...)
> > 13.
> > lapply(object, pdFactor)
> > 12.
> > unlist(lapply(object, pdFactor))
> > 11.
> > pdFactor.reStruct(object)
> > 10.
> > pdFactor(object)
> > 9.
> > unlist(pdFactor(object))
> > 8.
> > MEEM(object, conLin, control$niterEM)
> > 7.
> > Initialize.reStruct(X[[i]], ...)
> > 6.
> > FUN(X[[i]], ...)
> > 5.
> > lapply(object, Initialize, data, conLin, control)
> > 4.
> > Initialize.lmeStruct(lmeSt, dataMix, grps, control = controlvals)
> > 3.
> > Initialize(lmeSt, dataMix, grps, control = controlvals)
> > 2.
> > lme.formula(fixed = year_c ~ 1, random = ~1 | ID, data = data, method =
> "ML")
> > 1.
> > lme(fixed = year_c ~ 1, random = ~1 | ID, data = data, method = "ML")
> >
> >
> >
> > The code I used was:
> >
> > ### Unconditional Model
> > library(nlme)
> > library(broom)
> >
> > uncmod <- lme(fixed = year_c~1, random = ~1|ID, data = data, method
> ="ML")
> > tidy(uncmod)
> > summary(uncmod)
> > varCorr(uncmod)
> >
> >
> >
> > I hope this helps and the problem can be fixed soon, because right now,
> I am stopped dead in my work, unable using the nlme package, on which my
> analysis depends.
> >
> >
> > Thank you very much for your consideration, help, and good work!
> >
> >
> > Best,
> >
> >
> > Jochen Wirsing
>
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
-- 
Sent from Gmail Mobile

	[[alternative HTML version deleted]]