search for: oxboy

Am I missing something obvious on how to draw multi-line plots in base graphics? In ggplot2, I can do: data(Oxboys, package = "nlme") library(ggplot2) qplot(age, height, data = Oxboys, geom = "line", group = Subject) But in base graphics, the best I can come up with is this: with(Oxboys, plot(age, height, type = "n")) lapply(split(Oxboys[c("age", "height")],...

Dear all, I've got a problem in including square variables in lme and lmlist functions. I've tried to work on Oxboys data of Pinheiro and Bates'book, which consist of the heights of 26 boys, each mesured on nine different occasions : > Oxboys Grouped Data: height ~ age | Subject Subject age height Occasion 1 1 -1.0000 140.50 1 2 1 -0.7479 143.40 2 3 1 -0.4...

...math.ethz.ch > Subject: [R] [R ] Query : problems with the arithmetic operator "^" > with function "lme" and "lmList" > > Dear all, > > I've got a problem in including square variables in lme and lmlist > functions. I've tried to work on Oxboys data of Pinheiro and > Bates'book, which consist of the heights of 26 boys, each mesured on > nine different occasions : > > > Oxboys > Grouped Data: height ~ age | Subject > Subject age height Occasion > 1 1 -1.0000 140.50 1 > 2 1...

I have a number of x, y observations (Time, Conc) for a number of Subjects (with subject number Subj) and Doses. I can plot the individual points with xyplot fine: xyplot(Conc ~ Time | Subj, Groups=Dose, data=myData, panel = function(x,y) { panel.xyplot(x, y) panel.superpose(???) # Needs more here } ) I also like to plot on

Full_Name: Vaidotas Zemlys Version: 2.9.0 OS: Ubuntu 8.10 Submission from: (NULL) (213.197.173.50) Steps to reproduce the bug: library(nlme) data(Oxboys) qm=lmList(height~age|Subject,data=Oxboys) grid=with(Oxboys,expand.grid(age=seq(min(age),max(age),length=50),Subject=levels(Subject))) res <- predict(qm,grid,se=TRUE) Erreur dans if (pool) { : l'argument est de longueur nulle res <- predict(qm,grid,se=TRUE,pool=TRUE) Erreur dans sprint...

Hi all, I am trying to plot the fitted trajectories for each individual from an individual growth model (fit with a linear mixed effects model in lme). How can I plot each person's trajectory in the *same* panel, along with the mean-level trajectory? Below is an image of a plot similar to what I'm trying to create (from: http://jpepsy.oxfordjournals.org/content/31/10/1002/F6.large.jpg):

Dear R-users Coming from a proc mixed (SAS) background I am trying to get into the use of (n)lme. In this connection, I have some (presumably stupid) questions which I am sure someone out there can answer: 1) With proc mixed it is easy to get a hold on the estimated variance parameters as they can be put out into a SAS data set. How do I do the same with lme-objects? For example, I can see the

Dear R users, I was wondering if anyone knows how to obtain(subset) the first and/or the last record of a subject in a longitudinal setup. Normally in SAS one uses first.variable1 and last.variable1. So my question is that is there an R way of doing this. Regards, -- Luwis Diya, Phd student (Biostatistics), Biostatistical Center, School Of Public Health, Catholic University of Leuven, U.Z. St

I wish to draw a graph representing multiple series (sets of x,y points). Each series has its own label and points within each series are joined by a line ordered by their X cooridnate. I would also like a legend automatically showing which each series is. Which high-level command can serve this purpose? I looked at my book but can't find such a command. Thanks! Mark

...], type="l") =========================================================================== I thought that I should be able to do it following the example from pg.50 of Hadley's ggplot pdf book. ===========================Hadleys example=========================== library(nlme) # to get Oxboy dataset p <- ggplot(Oxboys, aes(age, height, group = Subject)) + geom_line() p ========================NON-WORKING CODE============================= bmm <- melt(bmi, id=c("pct")) names(bmm) <- c("age", "pp", "bm") bmm[,2] <- as.ordered(bmm[,2])...

I'm looking at Pinheiro & Bates "Mixed-Effects Models in S and S-PLUS" at the moment. Several datasets are used, one of which is called "PBIB" (a partially balanced incomplete block design). All the other datasets can be found somewhere or other in R. However, I cannot locate PBIB, and it does not seem to be mentioned in the latest edition of the R Full Reference

...e, we have about 800 children with measurements between 0 and 2.5 years. Unfortunatly, the number of the measurements as well as the timing differs a lot, it is only guranteed that each child is measured at least once in all of the 5 first quarters. Chosing a model like height ~ age, like in the Oxboys example from "Mixed Effext Models in S and S-PLUS" does not make sense since growth velocity changes rapidly in the first two years. So, I sticked to the spline approach, since I thought, that 3 knots would suffice this time. Unfortunatly, the computation would not stop after more tha...