análisis repetitivo<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Hola a todos,
Quiero procesar los datos contenidos en una tabla de acuerdo con el ejemplo adjunto. La tabla consta de 100 columnas + una columna que sirve como medida de control y de 30 a 80 filas.
El objetivo es aplicar una prueba estadística de una manera iterativa comenzando con control vs columna 51, control vs columna 52, control vs columna 53, etc ... para detectar cuándo aparece una diferencia significativa, cuánto tiempo está presente en las columnas sucesivas y cuando el significado se pierde.
Por ejemplo, se puede detectar una diferencia de la columna 63 a la columna 78, y luego se pierde de c79 a c86, y luego vuelve a aparecer en c87 a c98.
No sé cómo puedo escribir un script que permita un análisis automático de este tipo de tabla, dando un informe de los momentos en que ocurren las diferencias.
Cualquier sugerencia ?
Re: análisis repetitivo¡Es una forma constructiva de intercambiar consejos! Gracias por eso Kevin. Es interesante saber que podemos ir y venir con R y JMP. Seguramente voy a echar un vistazo más de cerca a los procedimientos de cambio de punto.
Re: análisis repetitivo¡Hola Pablo!
Sin más información, no estoy haciendo mucho más que disparar en la oscuridad.
Pero si, como sospecho, está buscando datos recopilados secuencialmente en un intento por discernir un cambio en el proceso generativo de los datos, un enfoque de punto de cambio podría estar más justificado.
Hay muchas formas de detectar puntos de cambio, y una gran cantidad de referencias que se remontan a muchos años atrás.
CRAN tiene un punto de cambio El paquete que implementa varios métodos recientemente investigados en R. JMP tiene una clara integración con R, en la que puede ejecutar código R en conjuntos de datos JMP y obtener los resultados de nuevo en JMP. Pruebe el método de tiempo lineal exacto podado (PELT) al que se hace referencia en R. Killick, P. Fearnhead y IA Eckley (2012) Detección óptima de puntos de cambio con un costo computacional lineal, diario de la American Statistical Association, 107: 500, 1590-1598.
Re: análisis repetitivoGracias Kevin, cualquier sugerencia es apreciada.
Según tengo entendido, Paul no eligió ninguna prueba todavía, pero estaba más preguntando cómo ejecutar una secuencia de pruebas de dos muestras (una muestra siempre es la misma) en JMP, dado su conjunto de datos original. Mi punto era que las secuencias de comandos no serían necesarias en este caso.
Gracias de nuevo por sus contribuciones.
Re: análisis repetitivoGracias por tu mensaje Kevin,
No quiero molestar a los estadísticos preeminentes si están vivos o muertos.
yo sería satisfecho a Sabes cual es tu sugerencia.
Re: análisis repetitivoProbablemente, alguien debería mencionar que el enfoque que están tomando puede causar que un número considerable de estadísticos preeminentes muertos giren como tornos en sus respectivas tumbas.
No puedo discernir exactamente qué decisión estás intentando tomar a partir de tu descripción, pero a mí me parece más bien una detección de cambio de punto. Si es así, analizar una cantidad de valores p agregados secuencialmente de cualquier prueba de dos muestras podría considerarse, en el mejor de los casos, una forma deficiente.
Re: análisis repetitivoHola Pablo,
Gracias, eso ayudó mucho.
Parece que el scripting no es obligatorio aquí. Podrías hacer lo siguiente:
1. Tablas> Pila: obtenga una nueva tabla de datos con datos de su columna de control y todas las columnas de medición apiladas. Esto creará una etiqueta de columna con sus nombres de columna y una segunda columna con todos sus datos.
El archivo adjunto muestra mi resultado.
2. Analizar> Ajustar Y por X, con X = Etiqueta e Y = Datos> Zona activa> Comparar medias> Con Control, Dunnett's: elija su columna de control y obtendrá todas las pruebas por pares. Puede hacer clic con el botón derecho en la tabla bajo Matriz de umbral de LSD, por ejemplo, para clasificarla o convertirla en una tabla de datos (tal vez para graficar valores p en Graph Builder).
Para la prueba de Dunnett vea también Comparar medios .
Espero que ayude,
Volker
Re: análisis repetitivoHola Volker y gracias por ocuparte de mi problema.
Aquí hay una tabla de datos ya analizada en la que la comparación Ctrl / prueba ha evidenciado diferencias significativas (evaluadas con en la prueba, pero sería más apropiado un Wilcoxon).
Cada fila a continuación proporciona la comparación entre los datos de Ctrl (columna de control) y una columna de prueba dada (numerada de 100 a 200)
Col nbr p valor
...
160 0.01
161 0,0022
162 ns
163 ns
164 0,03
165 ns
166 0,008
167 0,04
168 0,04
169 ns
170 0,004
171 ns
172 ns
173 ns
174 0,04
175 ns
176 <0,0001
177 0,008
178 ns
179 0,004
180 0,02
181 0,003
182 0,004
183 0,009
184 0,04
185 0,005
186 ns
187 0,006
188 0,03
189 0,02
190 0,0002
191 ns
192 0,001
etc ...
Re: análisis repetitivoHola Pablo,
Esto es lo que aprendí de nuestro chat:
Definición de una diferencia: la prueba de Wilcoxon compara los valores en la columna Ctrl con los valores en cada una de las otras 100 columnas (repitiendo comparaciones).
La tabla de ejemplo es solo para dar una idea de la estructura de datos. El comportamiento descrito anteriormente no se aplica a esta tabla.
El objetivo sería ordenar los tiempos de inicio y desplazamiento de una diferencia para cada tabla analizada de este tipo.
Si pudiera proporcionar otra tabla junto con algunos resultados de análisis que le gustaría ver, sería de gran ayuda.
BR, Volker