search for: gsummary

Full_Name: Ben Saville Version: 2.1 OS: Windows XP Submission from: (NULL) (152.2.94.145) I'm using the gsummary function to calculate a sum of V1 (column one) from my data 'mytest' by group (V2,or column 2). If V1 (the variable of interest) is all the same value (in this case all 2's), I do not get back the correct summation. If there is at least one difference in V1 (all 2's except for one...

Could someone give me a hand with the format of the gsummary function? Basically, I have a large set of xyz coordinates generated by LiDAR data (>37 million points) and I am trying to derive various summary statistics on the z-coordinates by a grid cell. I wrote a function to do this by creating factors from the x- and y- coordinates and then using gsum...

...d 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. [From mucking around &c using getAnywhere("augPred.lme"), and setting: debug(gsummary).] Further ferreting around showed that the data structures within gsummary() are fine, but that any (numeric only?) variable that has a label attached t...

Dear R list members, I've turned up a strange discrepancy between results obtained from the lme function in the nlme package in R and results obtained with lme in S-PLUS. I'm using version 3.1-24 of nlme in R 1.4.1 under Windows 2000, and both S-PLUS 2000 and 6.0, again under Windows 2000. I've noticed discrepancies in a couple of instances. Here's one, using data from Bryk

...all, Sometime I need to change data column name,but I catn't do it. example: type<-c("max","min") test<-data.frame(class=c("a","a","a","b","b"),value=1:5) for (i in 2:length(type)) { tmp<-gsummary(test,FUN=type[i],groups=test$class) tmp1<-data.frame(class=tmp$class, type[i] = tmp$value) ... error... : I catn't assign "max" by type[i] } Thanks Ken -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

Hi: Thanks a lot to Deepayan Sarkar, author of lattice I think, who solved my first query. I am afraid that I have another one. I am plotting several mutipanels boxplots (with one conditioning variable) on one page. The x, y and conditioning variable are all continuous variables. The x and conditioning variables are transformed to shingles before being plotted. The plot looks nice but there

I've got a data frame that looks like this: subject foo bar 2 1.7 3.2 2 2.3 4.1 3 7.6 2.3 3 7.1 3.3 3 7.3 2.3 3 7.4 1.3 5 6.2 6.1 5 3.4 6.9 ... That is, I've got multiple rows per subject. I need to compute summaries within categories where the subject has the same number

...; class that extends the data.frame class to describe relationships amoung the variables and the grouping of the observations. - multiple levels of grouping in the data and nested random effects in the models. - summarization of data according to grouping of the observations via the gsummary function (don't know if this works in R yet). - trellis-based plots of data and fitted models and derived quantities (definitely doesn't work in R yet) - much more stable computational methods, selection of starting values, etc. for fitting linear mixed-effects models - maximum lik...

Hi all! I am working on a nested lme model with one fixed effect ("treatment", which 3 levels) and two random effects for "Individuals" (four of them) within "treatment" and "replicate -2 levels-" within "individual" within "treatment". For doing so, I´ve been trying to create a factor for Individual%in%Treatment, say IT by

...[[<-`(`*tmp*`, nm, value = c(1L, 1L, 1L, 1L, 5L, 5L, 5L, 5L, 9L, 9L, 9L, 9L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 12L, 12L, 12L, 12L, 3L, 3L, 3L, 3L, 7L, 7L, 7L, 7L, 11L, 11L, 11L, 11L, 1L, 1L, 1L, 1L, 5L, 5L, 5L, 5L, 9L, 9L, 9L, 9L, 2L, 2L, 2L, 2L, 6L, 6L, 6L, 6L, 10L, 10L, 10L, 10L)) 3: gsummary(data, groups = groups) 2: augPred.lme(fit1, level = 0:1) 1: augPred(fit1, level = 0:1) > version _ platform x86_64-w64-mingw32 arch x86_64 os mingw32 system x86_64, mingw32 status Under development (unstable) major 3 minor...

...hy I get the number 3206, especially for > > depD27 depD35 and depD50. > > > > Thanks > > Han > > Han-Lin.Lai at noaa.gov > > By saying that dep is a level-2 variable do you mean that it does not > change within groups? Have you checked that, say by > > gsummary(mydata, form = ~ group, invariants = TRUE) -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not...

...ue = c(1L, 1L, 1L, 1L, 5L, 5L, 5L, 5L, > 9L, 9L, 9L, 9L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 12L, 12L, 12L, > 12L, 3L, 3L, 3L, 3L, 7L, 7L, 7L, 7L, 11L, 11L, 11L, 11L, 1L, > 1L, 1L, 1L, 5L, 5L, 5L, 5L, 9L, 9L, 9L, 9L, 2L, 2L, 2L, 2L, 6L, > 6L, 6L, 6L, 10L, 10L, 10L, 10L)) > 3: gsummary(data, groups = groups) > 2: augPred.lme(fit1, level = 0:1) > 1: augPred(fit1, level = 0:1) > >> version > _ > platform x86_64-w64-mingw32 > arch x86_64 > os mingw32 > system x86_64, mingw32 > status Under de...

..., 5L, 5L, 5L, > > 9L, 9L, 9L, 9L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 12L, 12L, 12L, > > 12L, 3L, 3L, 3L, 3L, 7L, 7L, 7L, 7L, 11L, 11L, 11L, 11L, 1L, > > 1L, 1L, 1L, 5L, 5L, 5L, 5L, 9L, 9L, 9L, 9L, 2L, 2L, 2L, 2L, 6L, > > 6L, 6L, 6L, 10L, 10L, 10L, 10L)) > > 3: gsummary(data, groups = groups) > > 2: augPred.lme(fit1, level = 0:1) > > 1: augPred(fit1, level = 0:1) > > > >> version > > _ > > platform x86_64-w64-mingw32 > > arch x86_64 > > os mingw32 > > system...

...arg$main <- "ATM Experiment: Frequency of Intrusion-none Button Presses During Error Trials\nAs A Function of Hierarchy Traversal Distance To Current Expected Button Press" arg$scales <- list(rot=c(90,90)) arg$sub <- date() arg$means <- list(data=gsummary(subset(arg$data,select=c("expected.hierarchyTo","trial.bi n","Freq")),form=~expected.hierarchyTo/trial.bin,mean)) arg$par.settings <- list(background=list(col="white"),strip.background=list(col="white")) arg$panel <- functi...

...L, >>> 9L, 9L, 9L, 9L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 12L, 12L, 12L, >>> 12L, 3L, 3L, 3L, 3L, 7L, 7L, 7L, 7L, 11L, 11L, 11L, 11L, 1L, >>> 1L, 1L, 1L, 5L, 5L, 5L, 5L, 9L, 9L, 9L, 9L, 2L, 2L, 2L, 2L, 6L, >>> 6L, 6L, 6L, 10L, 10L, 10L, 10L)) >>> 3: gsummary(data, groups = groups) >>> 2: augPred.lme(fit1, level = 0:1) >>> 1: augPred(fit1, level = 0:1) >>> >>>> version >>> _ >>> platform x86_64-w64-mingw32 >>> arch x86_64 >>> os mingw32...