plotMean3d: ~ Function: plotMeans3d for ClusterLongData3d ~
In kml3d: K-Means for Joint Longitudinal Data

plotMeans3d,ClusterLongData3d

R Documentation

~ Function: plotMeans3d for ClusterLongData3d ~

Description

Plot the means of two variables of a ClusterLongData3d object in 3D relatively to a Partition.

Usage



## S4 method for signature 'ClusterLongData3d,numeric'
plotMeans3d(x,y,varY=1,varZ=2,
   parTraj=parTRAJ(type="n"),parMean=parMEAN(),...)

Arguments

`x`	`[ClusterLongData3d]`: Object containing the trajectories to plot.
`y`	`[numeric]` or `[vector2(numeric)]`: Define the `Partition` P that will be use to plot the object. `P` is a `Partition` hold in the field `c2`, `c3`, ... `c26`. If `y=c(a,b)`, then `P` is the `Partition` number `b` with `a` clusters. If `y=a`, then `P` is the partition number `1` with `a` clusters.
`varY`	`[numeric]` or `[character]`: either the number or the name of the first variable to display. 1 by default.
`varZ`	`[numeric]` or `[character]`: either the number or the name of the second variable to display. 2 by default.
`parTraj`	`[ParLongData]`: Set the graphical parameters used to plot the trajectories of the `ClusterLongData3d`. See `ParLongData` for details.
`parMean`	`[ParLongData]`: Set the graphical parameters used to plot the mean trajectories of each clusters `ClusterLongData3d` (only when `y` is non missing). See `ParLongData` for details.

...

Arguments to be passed to methods, such as graphical parameters.

Details

Plot two variables of a ClusterLongData3d object in 3D. It use the rgl library. The user can make the graphical representation turn using its mouse.

Examples

### Move to tempdir
wd <- getwd()
setwd(tempdir()); getwd()

##################
### Real example on array

time=c(1,2,3,4,8,12,16,20)
id2=1:120
f <- function(id,t)((id-1)%%3-1) * t
g <- function(id,t)(id%%2+1)*t
h <- function(id,t)(id%%4-0.5)*(20-t)
myCld <- clusterLongData3d(array(cbind(outer(id2,time,f),outer(id2,time,g),
   outer(id2,time,h))+rnorm(120*8*3,0,3),dim=c(120,8,3)))
kml3d(myCld,3:4,2)

### Basic plot
plotMeans3d(myCld,3)

### plotMeans3d, variable 1 and 3
plotMeans3d(myCld,4,varZ=3)
plotMeans3d(myCld,3,parTraj=parTRAJ(col="red"))

### Go back to current dir
setwd(wd)

kml3d documentation built on Feb. 16, 2023, 9:44 p.m.