reduceNbTimes: ~ Function: reduceNbTimes ~

In kmlShape: K-Means for Longitudinal Data using Shape-Respecting Distance

Description

reduceNbTimes simplify some trajectories (in long format) by reducing their number of points.

Usage

  reduceNbTimes(trajLong, nbPoints, spar=NA)

Arguments

trajLong

[data.frame]: data.frame that hold the trajectories in long format. The data.frame has to be (no choice!) in the following format: the first column should be the individual indentifiant. The second should be the times at which the measurement are made. The third one should be the measurement.

nbPoints	[numeric]: fixe the number of that the simplified trajectories should have.
spar	[numeric]: smoothing parameter that is used if the trajectories shall be smoothed before being simplified.

Details

reduceNbTimes simplify some trajectories by reducing their number of points. The trajectories should be in long format. If a value is given to spar (different from NA), trajectories are smoothed using smooth.spline.

The reduction of the number of point is done using a variation of Douglas-Peucker algorithme based on the number of points instead of an epsilon.

Value

A data.frame holding the simplified trajectories, in long format.

Author(s)

Christophe Genolini

See Also

reduceNbTimes, DouglasPeuckerEpsilon, DouglasPeuckerNbPoints

Examples

  require(lattice)

  ### Some artificial data
  g <- function(x)dnorm(x,3)+dnorm(x,7)+x/10
  dn <- data.frame(id=rep(1:20,each=101),
    times=rep((0:100)/10,times=20),
    traj=rep(g((0:100)/10),20)+rep(runif(20),each=101)+rnorm(20*101,,0.1))

  xyplot(traj ~ times, data=dn, groups=id,type="l")

  ### Reduction to 50 points
  dn2 <- reduceNbTimes(trajLong=dn,nbPoints=50)
  xyplot(traj ~ times, data=dn2, groups=id,type="l")

  ### Reduction to 20 points
  dn3 <- reduceNbTimes(trajLong=dn,nbPoints=20)
  xyplot(traj ~ times, data=dn3, groups=id,type="l")

  ### Smoothing then reduction to 20 points
  dn4 <- reduceNbTimes(trajLong=dn,nbPoints=20,spar=0.5)
  xyplot(traj ~ times, data=dn4, groups=id,type="l")

Example output

Loading required package: class
Loading required package: longitudinalData
Loading required package: clv
Loading required package: cluster
Loading required package: rgl
Loading required package: misc3d
Loading required package: kml
Loading required package: lattice

Attaching package: 'kmlShape'

The following object is masked from 'package:longitudinalData':

    distFrechet

Warning messages:
1: In rgl.init(initValue, onlyNULL) : RGL: unable to open X11 display
2: 'rgl_init' failed, running with rgl.useNULL = TRUE 
3: .onUnload failed in unloadNamespace() for 'rgl', details:
  call: fun(...)
  error: object 'rgl_quit' not found 

kmlShape documentation built on May 1, 2019, 7:50 p.m.