similar to: nlme: augPred.lme for factor covariates

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"

Hi All, I sent the following message to R-help on July 14th, 2006: Let's say I have the following formula: a.formula <- x ~ y + z I want to extract the left and right-hand sides of the function so that I have two character vectors like the ones you would create using the following assignments: left.hand.side <- "x" right.hand.side <- c("y", "z") One

Dear R-help, I can't plot nls objects for some reason. The following example, taken from help(plot.nls) illustrates the problem: > data(Orthodont) > fm1 <- lme(distance ~ age, Orthodont, random = ~ age | Subject) > plot(fm1, resid(., type = "p") ~ fitted(.) | Sex, abline = 0) Error in do.call(plotFun, as.list(args)) : couldn't find function

Posted this to r-sig-mixed-models (https://stat.ethz.ch/pipermail/r-sig-mixed-models/2019q1/027620.html) but this might rather need to go to r-devel anyway, so reposting here: I came across a weird problem / bug with augPred() from nlme (nlme_3.1-137). Here is a reproducible example to illustrate the issue (tested on R 3.5.2 and R-devel 2019-03-03 r76192): library(nlme) dat <- data.frame(id

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,

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.

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),