análisis cosinor estadístico,

Question

Hola estoy tratando de calcular un análisis cosinor en statistica pero estoy en una pérdida en cuanto a cómo hacerlo. Necesito calcular el MESOR, amplitud y Acrophase de los datos ritmo ciracadian.

http://www.wepapers.com/Papers/73565/Cosinor_analysis_of_accident_risk_using__SPSS % 27s_regression_procedures.ppt

hay un enlace que muestra cómo lo hacen, las fórmulas y tal, pero no me ha dado mucha ayuda. ¿Alguien sabe el código para ello, ya sea en statistica o SPSS ??

Es realmente necesario para hacer esto porque es para un papel importante

Answer 1

No tengo SPSS o Statistica, así que no puedo decir el "push-este-botón" exacta tipo de medidas, pero tal vez esto ayudará.

análisis Cosinor se ajuste de una curva coseno (o seno) con un periodo conocido. La idea principal es que el problema no lineal de ajuste de una función coseno puede reducirse a un problema que es lineal en sus parámetros si se conoce el periodo. Voy a suponer que su periodo T = 24 horas.

1. Usted ya debe tener dos variables: Tiempo en la que se toma la medida, y Valor de la medición (éstos, por supuesto, podría ser llamado otra cosa).
2. Ahora crear dos nuevas variables: SINTIME = sen (2 x pi x Tiempo / 24) y CosTime = cos (2 x pi x Tiempo / 24) - esto se desribed en la página 11 de la presentación se enlazó (x es la multiplicación). Uso pi = 3,1415 si el valor exacto no está incorporada.
3. Ejecutar regresión múltiple lineal con Valor como resultado y SINTIME y CosTime como dos predictores. Usted debe obtener estimaciones de sus coeficientes, que llamaremos y B .
4. El término de intersección del modelo de regresión es la MESOR.
5. La amplitud es sqrt ( A ^ 2 + B ^ 2) [raíz cuadrada de A cuadrado más B al cuadrado]
6. El acrophase es arctan (- B / A ), donde arctan es la función inversa de bronceado. Las dos últimas fórmulas son de p.14 de la presentación.
7. El modelo de regresión debe también le dará un valor de R cuadrado para ver lo bien que el patrón circadiano 24 horas ajusta a los datos, y un valor de p global que las pruebas para la presencia de un componente circadiano con un período de 24 hrs.
8. Uno puede obtener errores estándar en amplitud y fase utilizando fórmulas propagación de error estándar, pero que no está incluido en la presentación.