Usage
convexLogisticPCA(x, k = 2, m = 4, quiet = TRUE, partial_decomp = FALSE, max_iters = 1000, conv_criteria = 1e-06, random_start = FALSE, start_H, mu, main_effects = TRUE, ss_factor = 4, weights, M)
Arguments
x
matrix with all binary entries
k
number of principal components to return
m
value to approximate the saturated model
quiet
logical; whether the calculation should give feedback
partial_decomp
logical; if TRUE
, the function uses the rARPACK package
to quickly initialize H
when ncol(x)
is large and k
is small
max_iters
number of maximum iterations
conv_criteria
convergence criteria. The difference between average deviance
in successive iterations
random_start
logical; whether to randomly inititalize the parameters. If FALSE
,
function will use an eigen-decomposition as starting value
start_H
starting value for the Fantope matrix
mu
main effects vector. Only used if main_effects = TRUE
main_effects
logical; whether to include main effects in the model
ss_factor
step size multiplier. Amount by which to multiply the step size. Quadratic
convergence rate can be proven for ss_factor = 1
, but I have found higher values
sometimes work better. The default is ss_factor = 4
.
If it is not converging, try ss_factor = 1
.
weights
an optional matrix of the same size as the x
with non-negative weights
M
depricated. Use m
instead