similar to: Weird problem / bug with augPred() from nlme

Dear Friends, Is there a workaround to this problem or must the object plotted by plot(augPred(object, ...)) be a grouped object? score <- c(108, 103, 96, 84, 118, 110, 129, 90, 84, 96, 105, 113, 96, 117, 107, 85, 125, 107, 128, 84, 104, 100, 114, 117, 110, 127, 106, 92, 125, 96, 123, 101, 100, 103, 105, 132, 122, 133, 107, 99, 116, 91, 128, 113, 88, 105, 112, 130) mcc

Dear all as I did not get any response on my post about abline and plot(augPred)) I try again. I hope I do not break some posting guide rules. I would try to contact package maintainer directly but there is stated to be R-core people, so I feel R-help list shall be OK. I need to draw straight lines through augPred plotted panels (vertical or horizontal) at specified point. I know I shall

Dear r-helpers, When I issue the commands plot(augPred(mcc.lme, primary = ~ mcc$age, length.out = 2, level = c (0, 1))) or plot(comparePred(mcc1.lis, mcc.lme, primary = ~ mcc$age, length.out = 2), layout = c(4, 3))) no lines are drawn. But fm1 <- lme(Orthodont) plot(augPred(fm1, level = 0:1, length.out = 2)) does plot lines. augPred() and comparePredict() seem to be working correctly, as

Dear R-users, I am a newbie to this site and a relative new-comer to S/R, so please tread lightly, for you tread... There have been several posting relating to problems with augPred() from the nlme library. Here is a "fix" for one of these problems which may lie at the root of others. In my case the problem with augPred() lay in gsummary(), which augPred() uses, causing it to fail.

I read the posts about augPred with lme, but does anyone know if there is a correlate for augPred for lmer? Specifically, I want to be able to use it to plot projections for all groups in an lmer class object using plot(augPred(lmer.object)). Thanks, Kyle *************************************** J. Kyle Roberts, Ph.D. Baylor College of Medicine Center for Educational Outreach One Baylor Plaza,

Dea-R community. I'd like to draw your attention to an issue I have recently encountered while doing my current data analysis. I've got an unexpected (to me) result from the command: > augPred(lmList(my.object)), 'my.object' being a grouped data frame of class: > class(my.object) [1] "nfnGroupedData" "nfGroupedData" "groupedData"

Hello, Is there a way to get augPred to work with lme if a subset statement is used? For example, if I modify the example from ?augPred.lme to include a subset statement, I get the following error: fm1 <- lme(Orthodont, random = ~1, subset=distance>19) augPred(fm1, length.out = 2, level = c(0,1)) Error in tapply(object[[nm]], groups, FUN[["numeric"]], ...) : arguments

All, I have a question RE plotting the prediction lines of a random effects model via augPred. I'll illustrate via the Pixel dataset within the nlme package: library(nlme) attach(Pixel) fm1Pixel = lme(pixel ~ day + I(day^2), data = Pixel, random = list(Dog = ~ 1)) plot(augPred(fm1Pixel)) ### 10 fitted lines since there are 10 dogs fm2Pixel = update(fm1Pixel, . ~ . + Side)

Hello, I want to use lme to fit two (or more) models, and then compare the fits on each individual. I know how to write my own code to do this (for each individual, plot the raw data, followed by lines() to plot each fitted curve) but I would like to use plot(augPred(... as it produces a nice trellis plot. I thought I could do this with par(new=T) but it does not seem to work.

Hi everybody, as you may be aware the function augPred.lme does not work as soon as the covariate is a factor. The problem lies in the line newprimary <- seq(from = minimum, to = maximum, length.out = length.out) which does not make sense for factors. I think augPred.lme can be useful for models with a factor covariate as well. Thus, I propose to change the code to: augPred.lme <-

Hello all, I am having a problem using augPred after I have created an nlme object. My data is 3933 rows long with no NA's, and is grouped with this function: vanfem<-groupedData(mass~age.years|name, data=av.mex, labels=list(x="Age", y="Mass"), units=list(x="(years)", y="(kg)")) Under the grouping factor "name", this is a character

Hi, I'm new at this, I'm very confused, and I think I'm missing something important here. In our pet example we have this: > fm <- lme(Orthodont) > plot(Orthodont) > plot(augPred(fm, level = 0:1)) which gives us a trellis plot with the females above the males, starting with "F03", "F04", "F11", "F06", etc. I thought the point of

Could someone explain the meaming of this error message from augPred: > augPred(area3.pen.nlme, primary=~day) Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop = FALSE], : Levels 1,2,3 not allowed for block > predict.nlme(area3.pen.nlme) does not produce an error. area3.pen.nlme was created with: > area3.pen.nlme <- nlme(area ~ SSlogis(day, Asym, xmid, scal),

Better posted on r-sig-mixed-models , no? Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Aug 23, 2017 at 5:17 AM, PIKAL Petr <petr.pikal at precheza.cz> wrote: > Dear all > > I encountered strange

> On Aug 23, 2017, at 8:08 AM, PIKAL Petr <petr.pikal at precheza.cz> wrote: > > Hi > > Well, yes I tried it about two weeks ago but my post did not get through as it still awaits moderator approval. It got through just fine. It appeared on Aug 15. It just didn't get any replies. As I read your original question in this thread, it was not clear to me that you had

Hi using this code example: library(nlme) fm1 <- lme(Orthodont, random = ~1) plot(augPred(fm1)) is there any way to have the plots in each cell labelled and ordered according to Orthodont$Sex? I.e. in addition to the bar with the label for Orthodont$Subject there is another bar labelling the Sex of the subject? thanks a lot christoph --

Using formulas constructed from strings only partially works for me in lme: library(nlme) data(Orthodont) fm2<-lme(as.formula("distance~age"),data=Orthodont,random=~1|Subject) summary(fm2) # works augPred(fm2) # fails #Error in inherits(object, "formula") : #Argument "object" is missing, with no default I assume that my use of as.formula is wrong, but

Dear all I encountered strange behaviour of augPred with virtually the same data First I made groupedData object. > mar.g<-groupedData(rutilizace~doba|int, data=mar) When I perform nlme on complete dataset I get an error with augPred > fit<-nlsList(rutilizace~SSasymp(doba, Asym, R0, lrc), data=mar.g) Warning message: c("1 error caught in nls(y ~ cbind(1 - exp(-exp(lrc) * x),

Hi Well, yes I tried it about two weeks ago but my post did not get through as it still awaits moderator approval. I could check which column is offending but actually it is only minor nuisance, I can live with selection of columns before fitting a model. What seems to me strange is that both full dataset and only selected colums gave me identical fit results but only one works within augPred.

Dear all I am trying nlme together with Pinheiro/Bates book. I constructed grouped data object with suitable plotting layout (according to some common factor, panels from bottom to top are in increasing order). When I do nlme(... some stuff...) I get fitted object which I can plot with plot(augPred(fit.nlme6, level=0:1)) but it results in completely different ordering. Is there any way