Re: [R-es] DUDA LLENAR MATRIZ CREADA
Hola David,
Puedes hacer como ya comentaron en las otras respuestas. Pero por
si sirve de utilidad, pego al final parte de un código que empleo para
ilustrar la repetición de contrastes en simulación.
Está con un bucle y guardando en vectores, si se quisiera modificar
para obtener una matriz y emplear lapply (o sapply), tendrías que crear
una función del tipo:
fun.test - function(muestra)
unlist(test(muestra)[c('statistic','p.value')])
reemplazando test por el contraste que te interesa (y suponiendo que
devuelve el objeto estándar de clase htest).
sapply(muestras, fun.test) podría ser lo que buscas (suponiendo que
muestras es una lista con los datos en cada componente). Para que
conste, para ciertas cosas yo no soy reacio a emplear bucles en lugar de
applies...
Un saludo, Rubén.
# --
# Repetición de contrastes
# --
# --
# Valores iniciales
set.seed(1) # Fijar semilla para reproductibilidad
n - 500
nsim - 1000
estadistico - numeric(nsim)
pvalor - numeric(nsim)
# --
# Realizar contrastes
for(isim in 1:nsim) {
u - runif(n)
tmp - ks.test(u, punif, min=0, max=1)
estadistico[isim] - tmp$statistic
pvalor[isim] - tmp$p.value
}
#Probar a cambiar por u - rnorm(n) y tmp - ks.test(u, pnorm, mean =
mean(u), sd = sd(u))
# Mantengo algo del resto del código por si resulta de interés...
# --
# Proporción de rechazos
cat(\nProporción de rechazos al 1% =, mean(pvalor 0.01), \n)
cat(Proporción de rechazos al 5% =, mean(pvalor 0.01), \n)
cat(Proporción de rechazos al 10% =, mean(pvalor 0.01), \n)
# --
# Análisis del estadístico contraste
# Histograma
hist(estadistico, freq=FALSE)
#...
# --
# Análisis de los p-valores
# (si todo 'correcto' con distribución uniforme)
# Histograma
hist(pvalor, freq=FALSE)
curve(dunif(x,0,1), add=TRUE) #abline(h=1)
# Test de Kolmogorov-Smirnov
ks.test(pvalor, punif, min=0, max=1)
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
[R-es] DUDA LLENAR MATRIZ CREADA
Buena tarde, Estoy llevando a cabo un trabajo y no encuentro la forma de llenar una matriz con el p_value y un estadístico calculado. Un poco mas detallado, tengo muestras aleatorias, calculo por ejemplo la prueba t, wilcoxon, etc y requiero llevar por una parte los p_value de k muestras a una matriz creada anteriormente y por otro lado el valor del estadístico de las mismas k muestras a otra matriz creada anteriormente. Las k muetras estarían determinadas por un ciclo. Agradezco me puedan ayudar con el tema que es de mucha urgencia. Saludos, David C. [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
Re: [R-es] DUDA LLENAR MATRIZ CREADA
Estimado David Contreras Mi respuesta es genérica, primero no se si usted se refiere a matriz o un data.frame, no por confusión pero si por como se trabaja en R (donde prefiere almacenar los resultados de). La segunda parte es sobre los resultados que desea almacenar en lo que usted dice matriz, en forma genérica puede usar $ y el nombre correspondiente, coloco un pequeño ejemplo según yo pienso que podría serle útil. x - c(1,2,3,4,5,6) y - c(3,4,5,3,2,5) datos - data.frame(x,y) modelo -lm(datos) # para ver que puedo pedir str(modelo) # del modelo pido coefficients modelo$coefficients El 25 de febrero de 2015, 16:14, David Contreras davidcontrera...@gmail.com escribió: Buena tarde, Estoy llevando a cabo un trabajo y no encuentro la forma de llenar una matriz con el p_value y un estadístico calculado. Un poco mas detallado, tengo muestras aleatorias, calculo por ejemplo la prueba t, wilcoxon, etc y requiero llevar por una parte los p_value de k muestras a una matriz creada anteriormente y por otro lado el valor del estadístico de las mismas k muestras a otra matriz creada anteriormente. Las k muetras estarían determinadas por un ciclo. Agradezco me puedan ayudar con el tema que es de mucha urgencia. Saludos, David C. [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
Re: [R-es] DUDA LLENAR MATRIZ CREADA
Hola,
A falta de que puedas dar más detalles con un ejemplo con el que podríamos
entender mejor el caso, si te he entendido bien el pseudo-código sería algo
así como esto.
#-
# Entiendo que el ciclo que comentas es con un bucle for.
# Vamos a guardar en cada iteración un p-value y el valor de Wilcoxon-test
matrizConestadisticos - 0
for(i in ...) {
# Guardo la iteración
matrizConestadisticos[ i,1] - i
#... a) aquí estaría tu cálculo del que extraes el p-value que guardamos en
matrizConestadisticos
matrizConestadisticos[ i, 2] - valor.p.value
#... b) aquí estaría tu cálculo del wilcoxon-test... que también
guardamos...
matrizConestadisticos[ i, 3] - varlos.wilcoxon
} #fin del bucle for
#---
Saludos,
Carlos Ortega
www.qualityexcellence.es
El 25 de febrero de 2015, 20:14, David Contreras davidcontrera...@gmail.com
escribió:
Buena tarde,
Estoy llevando a cabo un trabajo y no encuentro la forma de llenar una
matriz con el p_value y un estadístico calculado.
Un poco mas detallado, tengo muestras aleatorias, calculo por ejemplo la
prueba t, wilcoxon, etc y requiero llevar por una parte los p_value de k
muestras a una matriz creada anteriormente y por otro lado el valor del
estadístico de las mismas k muestras a otra matriz creada anteriormente.
Las k muetras estarían determinadas por un ciclo.
Agradezco me puedan ayudar con el tema que es de mucha urgencia.
Saludos,
David C.
[[alternative HTML version deleted]]
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
[[alternative HTML version deleted]]
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
Re: [R-es] DUDA LLENAR MATRIZ CREADA
David, No hay mucho detalle, pero aqui van algunas ideas: 1. Construye una funcion que tome como argumentos los parametros con los que generas tus muestras y construye una lista con los resultados. Cada componente de la lista contendra las *B *muestras sobre las cuales vas a realizar los analisis. 2. Utiliza lapply() para tomar cada entrada de la lista creada en 1. y aplica las pruebas que necesitas. Esta funcion debe tomar como argumentos (posiblemente) un data.frame() y realizar las pruebas. Por supuesto dentro de esa funcion solo extraerias los estadisticos de prueba y los valores p asociados a t.test(), wilcox.test, etc... 3. Usa do.call(rbind, L) para organizar los resultados. 4. Usa write.table() para exportarlos. Espero sea de ayuda. Saludos, Jorge.- 2015-02-26 6:14 GMT+11:00 David Contreras davidcontrera...@gmail.com: Buena tarde, Estoy llevando a cabo un trabajo y no encuentro la forma de llenar una matriz con el p_value y un estad�stico calculado. Un poco mas detallado, tengo muestras aleatorias, calculo por ejemplo la prueba t, wilcoxon, etc y requiero llevar por una parte los p_value de k muestras a una matriz creada anteriormente y por otro lado el valor del estad�stico de las mismas k muestras a otra matriz creada anteriormente. Las k muetras estar�an determinadas por un ciclo. Agradezco me puedan ayudar con el tema que es de mucha urgencia. Saludos, David C. [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
