4 votos

Es posible (o conveniente) para hacer el Cambio de Punto de Análisis en la secuencia de los grupos con R?

Estoy familiarizado con post-hoc de las pruebas con ANOVA para explorar las diferencias entre una secuencia de grupos, pero últimamente he estado leyendo sobre el Cambio de Punto de Análisis (especialmente la R paquetes bcp, changepoint y strucchange).

Parece que los paquetes sólo se ocupan de datos donde hay un punto de datos por unidad de tiempo. Tengo curiosidad por ver si pueden ser utilizados con los datos que existen múltiples puntos de datos por unidad de tiempo. He aquí un ejemplo de datos que representan la medición de una sola variable continua en un número de ejemplares que se han fechado a momentos específicos en el tiempo (no hay mediciones repetidas):

a<-data.frame(time="1000",x=rnorm(10,12,3))
b<-data.frame(time="2000",x=rnorm(50,13,4))
c<-data.frame(time="3500",x=rnorm(50,12,4))
d<-data.frame(time="5000",x=rnorm(7,14,5))
e<-data.frame(time="7000",x=rnorm(20,10,3))
f<-data.frame(time="7500",x=rnorm(15,11,3))
g<-data.frame(time="9000",x=rnorm(15,10,5))
h<-data.frame(time="9500",x=rnorm(35,30,2))
i<-data.frame(time="10000",x=rnorm(30,28,4))
a2i<-rbind(a,b,c,d,e,f,g,h,i) 

library(ggplot2)
a2i$time<-as.numeric(levels(a2i$time))[a2i$time] 
ggplot(a2i,aes(time,x))+stat_smooth()+geom_point()

enter image description here

He aquí lo que me gustaría ser más agradecido para algunos consejos sobre...

Q1. Sería válido para hacer el Cambio de Punto de Análisis en un vector como el de los medios o de las medianas de los grupos? Que me permitiría empezar con un punto de datos por unidad de tiempo " formato de entrada que se adapte a la R paquetes, como yo las entiendo. He visto hecho con datos ambientales como mensual de las concentraciones de gas (a partir de las observaciones diarias), pero pensé que me gustaría comprobar.

Q2. Hay un tipo de Cambio de Punto de Análisis que puedo hacer sobre los datos en bruto en a2g que me va a dar algunas medidas de las probabilidades de cambios a través de la secuencia? Por ejemplo, algo que detectará el cambio de tiempo=9000 a tiempo=9500, utilizando todos los puntos de datos en la muestra? Supongo que si era posible que alguien ya se han implementado y sólo necesito un puntero a la función correspondiente.

Q3. En el caso de Q2 puede ser la respuesta es 'sí', el método de cambio si la distribución de cada uno de los valores del grupo no era normal (a diferencia de mis datos de la muestra)?

T4. Si el Cambio de Punto de Análisis es completamente el enfoque equivocado aquí, por favor hágamelo saber. Estoy básicamente curiosos acerca de otros métodos de ANOVA de estos tipos de datos. Cualquier otra sugerencia sería la mayoría de la recepción.

1voto

Heath Puntos 1042

Para contestar una de tus preguntas:

  1. Sí, usted puede hacer changepoint análisis en un vector. Creo que el único paquete que actualmente pueden hacerlo en R es la ecp paquete de:

http://cran.r-project.org/web/packages/ecp/

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X