icbiplot: Interval-censored biplot

In icensBKL: Accompanion to the Book on Interval Censoring by Bogaerts, Komarek, and Lesaffre

Interval-censored biplot

Description

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

Usage

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.

Value

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

Author(s)

Silvia Cecere, port into icensBKL by Arnošt Komárek [e-mail]

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.

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) 

icensBKL documentation built on Sept. 19, 2022, 5:06 p.m.