search for: nfngroupeddata

...roupedData(estimate ~ slant | subject, data = vfr) vfrg.lmL <- lmList( estimate ~ (slant + respType + visField + hand)^2 | subject, vfrg) Error in function (classes, fdef, mtable) : unable to find an inherited method for function "lmList", for signature "formula", "nfnGroupedData" The data are at http://www.people.virginia.edu/~mk9y/LongFormat.csv To get the working data.frame: vf <- read.table("http://www.people.virginia.edu/~mk9y/ LongFormat.csv", header = TRUE, sep=',') names(vf) <- c( 'slant', 'estimate', 'respTyp...

Hi! Here is some output from an str command: > str(CO2) Classes ‘nfnGroupedData’, ‘nfGroupedData’, ‘groupedData’ and 'data.frame': 84 obs. of 5 variables: $ Plant : Ord.factor w/ 12 levels "Qn1"<"Qn2"<"Qn3"<..: 1 1 1 1 1 1 1 2 2 2 ... $ Type : Factor w/ 2 levels "Quebec","Mississippi": 1 1 1 1 1...

...1 1309 1 4.1 4 1 1 544 1 5.1 5 1 1 654 1 6.1 6 1 0 NA 1 7.1 7 1 1 441 1 8.1 8 1 1 882 1 9.1 9 1 1 1097 1 10.1 10 1 1 898 1 The data are clearly grouped. Just to be sure.... str(df.group) Classes 'nfnGroupedData', 'nfGroupedData', 'groupedData' and 'data.frame': 15480 obs. of 5 variables: $ id : Ord.factor w/ 43 levels "31"<"8"<"30"<..: 41 11 33 22 4 27 5 2 37 19 ... $ trial: int 1 1 1 1 1 1 1 1 1 1 ... $ ACC : int 1 0 1 1 1 0 1 1...

...tex example missing link(s): xyplot plot.augPred text html latex example missing link(s): xyplot plot.compareFits text html latex example <snip> plot.nffGroupedData text html latex example plot.nfnGroupedData text html latex example missing link(s): xyplot plot.nmGroupedData text html latex example <snip> These missing links occur only because S(PLUS) has xyplot but R does not, correct? That's what the old mailing list messages (1999) su...

Given a data frame of a categorical variable and two continuous variables, I would like to display one continuous variable against the other for each value of the categorical variable, all superimposed on the same plot. For example: data(Indometh); str(Indometh) Classes 'nfnGroupedData', 'nfGroupedData', 'groupedData' and 'data.frame': 66 obs. of 3 variables: $ Subject: Ord.factor w/ 6 levels "1"<"4"<"2"<"5"<..: 1 1 1 1 1 1 1 1 1 1 ... $ time : num 0.25 0.5 0.75 1 1.25 2 3 4 5 6 ... $ conc...

Hallo! I want to plot multiple grouped data in a postscript file using a loop. As I use a loop no plot (or just one empty plot) is generated. Here an example: library(nlme) data(Loblolly) # example data from nlme postscript("PSFile.ps") for (i in 1:1) # just as example { plot(Loblolly) } dev.off() Result: Just an empty PSFile.ps. (Withoput the loop it works. May anybody help? Karl

...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" "data.frame" The problem is more easily explained by showing these two to Trellis plot: > plot(my.object,some_options) > plot(augPred(lmList(my.object)),some_options) [link to output: http://www.palug.net/Members/jabba ]...

...lot.gls* [16] plot.hclust* plot.histogram* plot.HoltWinters* [19] plot.intervals.lmList* plot.isoreg* plot.lm [22] plot.lme* plot.lmList* plot.medpolish* [25] plot.mlm plot.nffGroupedData* plot.nfnGroupedData* [28] plot.nls* plot.nmGroupedData* plot.pdMat* [31] plot.ppr* plot.prcomp* plot.princomp* [34] plot.profile.nls* plot.ranef.lme* plot.ranef.lmList* [37] plot.shingle* plot.simulate.lme* plot.spec...

Dear All, I would be grateful if you can tell how can I suppress the legend (automatically created) when I plot a groupedData. suppose that I have data farme df_0 which contains the following columns: ID represents the clusters, t represents time of observation, Y: the variable of interest and X a binary covariate. Then: df_1 <- groupedData(Y ~ t | ID) plot(df_1,

Hi, R users, I'm using the "lme" function in R to estimate a 2 level mixed effects model, in which the size of the subject groups are different. It turned out that It takes forever for R to converge. I also tried the same thing in SPSS and SPSS can give the results out within 20 minutes. Anyone can give me some advice on the lme function in R, especially why R does not converge?

Hi, I would be grateful for help with a problem which is irritating me. I am quite sure that I am doing something stupid, but I can't see what it is. I am running R 1.7 on Windows 2000. The graphics device is the PC screen. The graphics from the nlme demonstration in Bates an Pinheiro's manual work just as advertised. The CO2 data and the Orthodont data dsiplay

...77L, 95L, 96L, 97L, 98L, 99L, 100L, 101L, 102L, 103L, 104L, 105L, 106L, 107L, 108L, 109L, 110L, 127L, 128L, 129L, 130L, 131L, 132L, 133L, 134L, 135L, 136L, 137L, 138L, 139L, 140L, 141L, 159L, 160L, 161L, 162L, 163L, 164L, 165L, 166L, 167L, 168L, 169L, 170L, 171L, 172L, 173L, 174L ), class = c("nfnGroupedData", "nfGroupedData", "groupedData", "data.frame"), formula = Prop ~ end | temp, FUN = function (x) max(x, na.rm = TRUE), order.groups = TRUE) > > sessionInfo() R Under development (unstable) (2017-07-31 r73003) Platform: x86_64-w64-mingw32/x64 (64-bit) Runnin...

> 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 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.