¿Cómo puedo encontrar las observaciones asociadas con el rango j en este código?

Answer

1. 

   davy.ai

   June 12, 2023 at 4:53 am

   Para encontrar las observaciones asociadas con el rango j en el Muestreo de Conjuntos Clasificados utilizando simulaciones, se puede modificar la función rss para que devuelva los conjuntos clasificados generados para cada simulación. Luego, se puede utilizar la función rank para encontrar el rango de cada observación dentro de cada conjunto. Finalmente, se pueden extraer las observaciones j de cada conjunto para obtener los valores Xj.

   Aquí hay un ejemplo de código que implementa este enfoque:

   # Función rss modificada que devuelve los conjuntos clasificados para cada simulación
   rss_sims <- function(X, m, r, sims) {
     sets <- matrix(0, ncol = m * r, nrow = sims)
   
     for (i in 1:sims) {
       ranked_sets <- rankedsets(X, m)
       sets[i,] <- matrix(ranked_sets, ncol = m * r, byrow = TRUE)
     }
   
     rownames(sets) <- paste0("Simulación ", 1:sims)
   
     return(sets)
   }
   
   # Ejemplo de uso
   datos <- rnorm(10000, 1, 3)
   sims <- 1000
   sets <- rss_sims(datos, m = 5, r = 3, sims = sims)
   
   # Encontrar valores Xj para el rango j = 3
   rango_j <- 3
   mtx_rango <- matrix(apply(sets, 1, rank), ncol = 5)
   Xj <- apply(mtx_rango, 1, function(x) sets[which(x == rango_j),][[1]])
   
   # Mostrar estadísticas resumidas de los valores Xj
   summary(Xj)
   

   En este código, primero definimos una versión modificada de la función rss llamada rss_sims que devuelve los conjuntos clasificados generados para cada simulación. Luego, generamos sims conjuntos clasificados usando esta función y guardamos los conjuntos resultantes en una matriz llamada sets.

   Para encontrar los valores Xj para el rango j = 3 (puedes reemplazar esto con cualquier otro rango), primero utilizamos la función rank para calcular el rango de cada observación dentro de cada conjunto. Luego, extraemos la observación j de cada conjunto utilizando una combinación de las funciones which y apply, y guardamos los valores resultantes de Xj en el vector Xj. Finalmente, mostramos las estadísticas resumidas de los valores Xj utilizando la función summary.

Comments are closed.