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

Displaying 20 results from an estimated 300 matches similar to: "Weird problem / bug with augPred() from nlme"

2006 Nov 10
0
Can plot(augPred()) be used with objects that are not groupedData?
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
2006 Sep 04
1
abline and plot(augPred) help
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
2006 Nov 09
0
Help with plot(augPred()) plot(comparePred()) in nlme
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
2006 May 17
1
Fix for augPred/gsummary problem (nlme library)
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.
2007 Mar 08
2
augPred in lmer
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,
2008 Oct 08
1
Observed responses in 'augPred' data frame - Wrong order ?
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"
2004 Nov 30
1
augPred with lme(...,subset=...)
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
2006 Sep 14
1
plotting all subgroups with augPred
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)
2011 Jul 07
2
How do I overlay two trellis plots of lme fitted lines produced by plot.augPred?
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.
2010 Jul 28
0
nlme: augPred.lme for factor covariates
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 <-
2008 Mar 24
0
Error using augPred with nlme grouped data
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
2006 Jul 03
1
panel ordering in nlme and augPred plots
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
2001 Aug 08
1
NLME augPred error
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),
2017 Aug 23
0
strange nlme augpred behaviour
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
2017 Aug 24
0
strange nlme augpred behaviour
> 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
2005 Sep 29
1
plot.augPred sorted and labelled according second factor
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 --
2003 Jan 30
1
as.formula(string) and augPred in lme
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
2017 Aug 23
2
strange nlme augpred behaviour
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),
2017 Aug 23
2
strange nlme augpred behaviour
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.
2005 Apr 25
4
panel ordering in nlme and augPred plots
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