[R] add spline to longitudinal data - preferably similar to SAS's 'I=SM50S' routine

Dylan Beaudette dylan.beaudette at gmail.com
Tue Jan 19 00:49:41 CET 2010


On Mon, Jan 18, 2010 at 3:25 PM, Eric Fail <e at it.dk> wrote:
> Hi Ruser
>
> I'm trying to replicate some SAS code. I have to add a spline to my
> longitudinal spaghetti plot.
>
> I have the plot, but I can't add the spline, a overall trend line. In the
> SAS code they use the command   'I=SM50S' and I would prefer something
> similar. I’m using R 2.10.1 on windows XP…
>
> I have made this working example.
>
> tolerance.pp <-
> read.table("http://www.ats.ucla.edu/stat/R/examples/alda/tolerance1_pp.txt",
> sep=",", header=T)
> # install.packages("lattice", dep = T)
> library(lattice)
> xyplot(tolerance ~ age, groups = id, data=tolerance.pp, type = "l")

Hi,

How about something like:

# your data
tolerance.pp <-
read.table("http://www.ats.ucla.edu/stat/R/examples/alda/tolerance1_pp.txt",
sep=",", header=T)

# almost works...
xyplot(tolerance ~ age, groups = id, data=tolerance.pp, type = c('l','smooth'))

see ?panel.loess for more ideas

Cheers,
Dylan

>
> This is where I want to add a overall spline.
>
> Hope someone out there can figure this out.
>
> Thanks
>
> Eric
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>



More information about the R-help mailing list