Visualisation of raw data for a smoothing-splines mixed-effects model

Share:

Description

A visualisation of the raw data used to fit the smoothing-splines mixed-effects model is obtained. This is presented as a trellis plot where each panel shows the observations for one of the subjects as a scatter plot (with optional lines connecting the observations if showRawLines is TRUE). If showModelFits is TRUE, then these observations are overlaid with the fitted individual curve for that subject.

Usage

1
plotSmeRaw(x,xlab="Time",ylab="Y",mainTitle="",showModelFits=TRUE,showRawLines=FALSE)

Arguments

x

an object inheriting from class sme, representing a fitted smoothing-splines mixed-effects model

xlab

the label for the x-axis, defaulting to "Time"

ylab

the label for the y-axis, defaulting to "Y"

mainTitle

the main title for the plot, defaulting to none

showModelFits

if TRUE, the fitted individual curve for each subject will be overlaid on their observations in their panel

showRawLines

if TRUE, the raw observations will be connected by lines

Value

A trellis plot as described above.

Author(s)

Maurice Berk maurice@mauriceberk.com

References

Berk, M. (2012). Smoothing-splines Mixed-effects Models in R. Preprint

See Also

plot.sme, plotSmeModel, plotSmeDiagnostic

Examples

1
2
3
4
5
  data(MTB)
  fit <- sme(MTB[MTB$variable==6031,c("y","tme","ind")])
  plotSmeRaw(fit,showModelFits=FALSE)
  plotSmeRaw(fit,showModelFits=FALSE,showRawLines=TRUE)
  plotSmeRaw(fit)