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 =
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
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]]