similar to: change fitted line colour in lme() trellis plot?

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.

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"

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

> 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

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.

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

All, I'm trying to create an augPred plot in the nlme library, similar to the plot on p.43 of Pinheiro & Bates (Mixed Effects Models in S and S-Plus) for their Pixel data. My data structure is the same as the example but I still get the error msg below. > comp.adj.UKV <- groupedData(adj.UKV ~ Time | Patient_no/Lisinopril, data = comp.adj.UKV.frm, order.groups = F) >

Hello everyone I've got another question on Trellis plots. I plotted my fitted model with plot(augPred(fitted model)). According to the description of this method, one can use any optional argument which is passed down to the xyplot function for Trellis plots. Additionally to the trivial plot i would like to insert a horizontal line. I tried to do it with the argument panel= function (...)

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

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

If I make a scatterplot and several (e.g. 5) points lie on top of each other at a given x,y location I would like the plot symbol to be the number of superimposed points (e.g. "5"). Could someone please tell me how to do this in R? Thanks! Bill Simpson

Dear List: I have a data set with over 7000 students with about 4 observations over time per student. I want to examine the within-group fits of a random sample of this group as it takes forever to compute and draw all 7000 regressions. Here is the code I have used so far. >group<-groupedData(math~year|childid, data=scores) >group.list<-lmList(group)

Using R v. 1.7.0 on Windows 2000 I would like to plot the fitted values of a model as a function of a continuous covariate, augmented with data (e.g., augPred) grouping by combinations of fixed effects. I have not been able to use augPred effectively, and am wondering if it does not handle unbalanced data (3 out of 192 missing). I include below the model and an xyplot that almost does the

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

Dear alltogether, how can a delay like possible with par(ask=TRUE) be attained while using trellis-plots within a loop or something like that? the following draws each plot without waiting for a signal (mouse-klick), so par() does not work for that: library(nlme) for(i in 1:3) { fitlme <- lme(Orthodont) par(ask=TRUE) # does not work with trellis.... print(

Dear R-supporters, I have results of lme similar to those shown in Fig. 1.21, p.51 of Pinheiro/Bates. However, In my data set, one of the panels is missing, leading to an ugly frame shift of the following panels. How can I replace one of the panels by an empty one to restore the raster? Dieter Menne --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel

Dear all I try to refine my nlme models and with partial success. The model is refined and fitted (using Pinheiro/Bates book as a tutorial) but when I try to plot plot(augPred(fit4)) I obtain Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop = FALSE], : Levels (0,3.5],(3.5,5],(5,7],(7,Inf] not allowed for vykon.fac > Is it due to the fact that I have unbalanced

I generate a groupedData object library(nlme) obj <- groupedData(mg10 ~ time | gp, data = common, outer = ~pct) gp has 101 levels, and pct has 3. There are 38, 25, 38 gps in each of the levels of pct respectively. I fit my model fit.rtg <- lme(mg10 ~ time * group, data = obj, random = ~time * group | gp) Now I try to plot the results. I would like to print 40 panels on each

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