[R-sig-ME] nlme issue since 3.40

[Jochen Wirsing]

Dear All,


thanks again for your nice feedback. The Session Info etc was mentioned
in the previous email, but one of the folks here on the list asked me to
use dput() so I could paste my data in textform here.

So here you go:

structure(list(Place = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label =
c("Cambridge Bay",
"Iqaluit", "Rankin Inlet"), class = "factor"), year_c = c(-2,
-1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, -2,
-1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, -2,
-1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15), SR_Total_ctd
= c(NA,
NA, 5.48961420000001, 7.24637679999999, 3.76569040000001, 6.9014085,
9.3088858, 5.1771117, 7.51366120000002, 7.18157180000001, 3.7275064,
5.7397959, 6.87500000000001, 6.25, 5.84577109999999, 6.2421973,
4.70446319999999, NA, NA, NA, 11.9192688, 9.6291013, 8.54640980000001,
8.7893297, 7.42343539999999, 5, 5.3674121, 6.1322261, 6.4182843,
5.75692960000001, 6.3861534, 6.6706021, 7.76642339999999, 6.91983119999999,
5.3634631, NA, NA, NA, 0.0870321999999959, -1.19047620000001,
-1.8723404, -0.4244482, 1.18443319999999, 2.19594590000001,
1.57545610000001,
3.080766, 1.8744906, 5.51053480000002, 6.31412790000002, 7.4451411,
8.1226054, 9.52743899999999, 9.4339623, NA)), .Names = c("Place",
"year_c", "SR_Total_ctd"), row.names = c(NA, -54L), class = "data.frame")
> 


I was also asked, whether using lme4 would be an option, but I don't
think that this will help fix the issue that I believe occurred through
updating to R 3.40. Also, nlme provides some information in a nicer way
(covariances etc, if I remember correctly - all this is still very new
to me....), so I'd really rather see nlme fixed.


Thanks again!


	[[alternative HTML version deleted]]