| bigarx_boot | R Documentation |
Calcula intervalos de confianza bootstrap para los parámetros de un
biplot, siguiendo la metodología de Nieto, Galindo, Leiva y
Vicente-Galindo (2014), implementada originalmente en biplotbootGUI.
bigarx_boot(
X,
B = 1000,
method = c("hj", "gh", "jk", "sqrt"),
transform = c("standardize", "center", "double", "grand", "none"),
alpha = NULL,
axes = c(1, 2),
conf = 0.95,
ci = c("percentile", "normal"),
seed = NULL,
groups = NULL,
progress = interactive()
)
X |
Matriz o data.frame numérico. |
B |
Número de réplicas bootstrap (por defecto 1000). |
method, transform, alpha |
Como en |
axes |
Plano en el que se calculan coordenadas, calidades y ángulos. |
conf |
Nivel de confianza (por defecto 0,95). |
ci |
|
seed |
Semilla opcional para reproducir los resultados. |
groups |
Factor opcional para colorear el biplot. |
progress |
Mostrar una barra de progreso. |
Algoritmo. Se extraen B muestras de individuos (filas) con
reemplazo. En cada una se repite el preprocesamiento (por ejemplo, se
vuelve a estandarizar) y el biplot, y se guardan los parámetros:
inercia: porcentaje de inercia de cada eje;
bondad: bondad de ajuste de la matriz en el plano;
coordenada: coordenadas de cada variable en cada eje del plano;
calidad: calidad de representación de cada variable en el plano;
contribucion: contribución de cada variable a cada eje;
angulo_variables: ángulo entre cada par de variables (en grados);
angulo_eje: ángulo entre cada variable y cada eje (en grados).
Alineación de ejes. La SVD de cada réplica puede devolver un eje con
el signo invertido, lo que inflaría artificialmente los intervalos de las
coordenadas y de los ángulos con los ejes. Por eso, en cada réplica se
invierte el signo del eje si sus cargas apuntan en sentido contrario a las
del análisis original. Los ángulos entre variables no dependen del signo.
Si dos valores propios son muy parecidos, los ejes pueden intercambiarse
entre réplicas; conviene revisar inercia antes de interpretar.
Intervalos (ci):
"percentile" (boot-q): cuantiles \alpha/2 y 1-\alpha/2 de
las réplicas (Algoritmo 3 del artículo);
"normal" (boot-t): estimación \pm z_{1-\alpha/2} por el error
estándar bootstrap, con t_{n-1} si n < 25 (Algoritmo 2).
Un objeto "bigarx_boot" con
fitEl biplot original (bigarx()).
tabledata.frame con una fila por parámetro: tipo, elemento, eje, estimación, media y error estándar bootstrap, sesgo, límites del intervalo y método.
replicatesMatriz de réplicas (una fila por réplica, una columna por parámetro).
B, conf, ci, axesConfiguración.
Nieto, A. B., Galindo, M. P., Leiva, V. y Vicente-Galindo, P. (2014). A methodology for biplots based on bootstrapping with R. Revista Colombiana de Estadística, 37(2), 367-397. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.15446/rce.v37n2spe.47944")}
Efron, B. y Tibshirani, R. J. (1993). An Introduction to the Bootstrap. Chapman & Hall.
bt <- bigarx_boot(USArrests, B = 200, seed = 1)
bt
summary(bt, type = "angulo_variables")
plot(bt)
plot(bt, type = "angulo_variables")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.