icbiplot function

Interval-censored biplot

Principal Component Analysis for interval-censored data as described in Cecere, Groenen and Lesaffre (2013).

icbiplot(L, R, p = 2, MaxIter = 10000, tol = 1e-06, plotit = TRUE, seed = NULL, ...)

Arguments

• L: Matrix of dimension number of individuals/samples by number of variables with left endpoints of observed intervals.
• R: Matrix of dimension number of individuals/samples by number of variables with right endpoints of observed intervals.
• p: Dimension of the solution. Default value is $p = 2$.
• MaxIter: Maximum number of iterations in the iterative minimazation algorithm
• tol: Tolerance when convergence is declared
• plotit: Logical value. Default equals TRUE. A biplot in dimension 2 is plotted.
• seed: The seed for the random number generator. If NULL, current R system seed is used.
• ...: further arguments to be passed.

Returns

Returns a list with the following components - X: matrix of number of individuals times 2 (p) with coordinates representing the individuals

• Y: matrix of number of variables times 2 (p) with coordinates representing the variables

• H: matrix of number of individuals times number of variables with approximated events

• DAF: Disperssion accounted for (DAF) index

• FpV: matrix showing the fit per variable

• iter: number of iterations performed

References

Cecere, S., Groenen, P. J. F., and Lesaffre, E. (2013). The interval-censored biplot. Journal of Computational and Graphical Statistics, 22 (1), 123-134.

Author(s)

Silvia Cecere, port into icensBKL by Arnošt Komárek arnost.komarek@mff.cuni.cz

Examples

data("tandmob", package = "icensBKL")
Boys <- subset(tandmob, fGENDER=="boy")
L <- cbind(Boys$L14, Boys$L24, Boys$L34, Boys$L44)
R <- cbind(Boys$R14, Boys$R24, Boys$R34, Boys$R44)
L[is.na(L)] <- 0
R[is.na(R)] <- 20    ## 20 = infinity in this case

icb <- icbiplot(L, R, p = 2, MaxIter = 10000, tol = 1e-6,
                plotit = TRUE, seed = 12345)