Buena tarde amigos, En días pasados hice algunas consultas y ya pude salir de las dudas que tenia en ese momento, ahora requiero de su colaboración con lo siguiente: Tengo un vector dicotomico (Binario) con la siguiente información que me surgio de algunos procesos anteriores:> MuestraS[1] 1 1 1 1 0 1 1 0 1 1 Ahora necesito hallar todas las posibles combinaciones que se puedan hacer con estos elementos para luego hacer un muestro aleatorio simple con reemplazo y seleccionar algunas de las posibles muestras que se obtengan. Agradezco me puedan ayudar con este asunto Saludos, DC. [[alternative HTML version deleted]]
Hola David, ¿Seguro que buscas las combinaciones? Creo que lo que buscas es esto... #---------------- MuestraS <- c(1 ,1, 1, 1, 0, 1, 1, 0, 1, 1) library(combinat) resPer <- permn(MuestraS) matresPer <- matrix(unlist(resPer), nrow=factorial(length(MuestraS)), ncol=length(MuestraS)) head(matresPer) #---------------- Que produce esto:> matresPer <- matrix(unlist(resPer), nrow=factorial(length(MuestraS)),ncol=length(MuestraS))> head(matresPer)[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 1 1 0 1 1 1 1 0 1 1 [2,] 1 1 1 0 1 1 1 1 0 1 [3,] 1 0 1 1 1 1 0 1 1 1 [4,] 1 1 1 1 0 1 1 1 1 0 [5,] 0 1 1 1 1 0 1 1 1 1 [6,] 1 1 1 1 1 1 1 1 1 1 Cuidado que el objeto "resPer" es una lista de 3628800 elementos... 609.1Mb que como matriz sólo ocupa 276.9Mb. Saludos, Carlos Ortega www.qualityexcellence.es El 2 de marzo de 2015, 22:51, David Contreras <davidcontreras00 en gmail.com> escribió:> Buena tarde amigos, > > En días pasados hice algunas consultas y ya pude salir de las dudas que > tenia en ese momento, ahora requiero de su colaboración con lo siguiente: > > Tengo un vector dicotomico (Binario) con la siguiente información que me > surgio de algunos procesos anteriores: > > > MuestraS > [1] 1 1 1 1 0 1 1 0 1 1 > > Ahora necesito hallar todas las posibles combinaciones que se puedan hacer > con estos elementos para luego hacer un muestro aleatorio simple con > reemplazo y seleccionar algunas de las posibles muestras que se obtengan. > > Agradezco me puedan ayudar con este asunto > > Saludos, > > DC. > > [[alternative HTML version deleted]] > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]
Hola David, No sé exactamente que quieres hacer pero aparentemente con sample() podría ser suficiente. Como puedes ver en el código de Carlos, con factorial() obtendrías el numero de permutaciones, aunque no se para que lo querrías (supongo que con 1000 podría ser suficiente). Para el muestreo te recomiendo fijar la semilla con set.seed () y después llamar a sample las veces que necesites (puedes combinarlo con replicate). Por ejemplo con: nsim <- 1000 set.seed(1) res <- replicate(nsim, sample(MuestraS)) obtendrías nsim replicas aleatorias (en columnas; equivaldrían a nsim filas aleatorias de la matriz matresPer del correo de Carlos, si no me equivoco...). Un saludo, Rubén. El 02/03/2015 22:51, "David Contreras" <davidcontreras00 en gmail.com> escribió:> Buena tarde amigos, > > En días pasados hice algunas consultas y ya pude salir de las dudas que > tenia en ese momento, ahora requiero de su colaboración con lo siguiente: > > Tengo un vector dicotomico (Binario) con la siguiente información que me > surgio de algunos procesos anteriores: > > > MuestraS > [1] 1 1 1 1 0 1 1 0 1 1 > > Ahora necesito hallar todas las posibles combinaciones que se puedan hacer > con estos elementos para luego hacer un muestro aleatorio simple con > reemplazo y seleccionar algunas de las posibles muestras que se obtengan. > > Agradezco me puedan ayudar con este asunto > > Saludos, > > DC. > > [[alternative HTML version deleted]] > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es >[[alternative HTML version deleted]]
Muchas gracias por sus respuestas tan oportunas y precisas, Use partes de los dos códigos y hasta el momento todo va perfecto. Nuevamente gracias y estaré en contacto ante cualquier inconveniente que se me presente. Saludos, DC. El 2 de marzo de 2015, 17:50, Carlos Ortega <cof en qualityexcellence.es> escribió:> Hola David, > > ¿Seguro que buscas las combinaciones? > > Creo que lo que buscas es esto... > > #---------------- > MuestraS <- c(1 ,1, 1, 1, 0, 1, 1, 0, 1, 1) > library(combinat) > resPer <- permn(MuestraS) > > matresPer <- matrix(unlist(resPer), nrow=factorial(length(MuestraS)), > ncol=length(MuestraS)) > head(matresPer) > #---------------- > > Que produce esto: > > > matresPer <- matrix(unlist(resPer), nrow=factorial(length(MuestraS)), > ncol=length(MuestraS)) > > head(matresPer) > [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] > [1,] 1 1 0 1 1 1 1 0 1 1 > [2,] 1 1 1 0 1 1 1 1 0 1 > [3,] 1 0 1 1 1 1 0 1 1 1 > [4,] 1 1 1 1 0 1 1 1 1 0 > [5,] 0 1 1 1 1 0 1 1 1 1 > [6,] 1 1 1 1 1 1 1 1 1 1 > > Cuidado que el objeto "resPer" es una lista de 3628800 elementos... > 609.1Mb que como matriz sólo ocupa 276.9Mb. > > Saludos, > Carlos Ortega > www.qualityexcellence.es > > El 2 de marzo de 2015, 22:51, David Contreras <davidcontreras00 en gmail.com> > escribió: > >> Buena tarde amigos, >> >> En días pasados hice algunas consultas y ya pude salir de las dudas que >> tenia en ese momento, ahora requiero de su colaboración con lo siguiente: >> >> Tengo un vector dicotomico (Binario) con la siguiente información que me >> surgio de algunos procesos anteriores: >> >> > MuestraS >> [1] 1 1 1 1 0 1 1 0 1 1 >> >> Ahora necesito hallar todas las posibles combinaciones que se puedan hacer >> con estos elementos para luego hacer un muestro aleatorio simple con >> reemplazo y seleccionar algunas de las posibles muestras que se obtengan. >> >> Agradezco me puedan ayudar con este asunto >> >> Saludos, >> >> DC. >> >> [[alternative HTML version deleted]] >> >> _______________________________________________ >> R-help-es mailing list >> R-help-es en r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-help-es >> > > > > -- > Saludos, > Carlos Ortega > www.qualityexcellence.es >[[alternative HTML version deleted]]