Analyse de cosinor statistique,

Question

Hey je suis en train de calculer une analyse de cosinor STATISTICA, mais je suis à une perte quant à la façon de le faire. Je dois calculer la mésor, AMPLITUDE et acrophase des données de rythme ciracadian.

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

il y a un lien qui montre comment le faire, les formules et autres, mais il ne m'a pas donné beaucoup d'aide. Est-ce que quelqu'un sait le code pour elle, que ce soit dans ?? STATISTICA ou SPSS

J'ai vraiment besoin pour obtenir ce fait parce qu'il est un document important

Answer 1

Je n'ai pas SPSS ou Statistica, donc je ne peux pas vous dire exactement « push-ce-bouton » genre d'étapes, mais peut-être cette volonté aide.

analyse Cosinor est ajustement d'une courbe cosinus (ou sinusoïdale) avec une période connue. L'idée principale est que le problème non-linéaire de mise en place d'une fonction cosinus peut être réduite à un problème qui est linéaire dans ses paramètres si la période est connue. Je suppose que votre période T = 24 heures.

1. Vous devez déjà avoir deux variables: Heure à laquelle la mesure est prise, et Valeur de la mesure (ceux-ci, bien sûr, peut être appelé autre chose).
2. Maintenant, créez deux nouvelles variables: SINTIME = sin (2 x pi x Heure / 24) et CosTime = cos (2 x pi x Heure / 24) - c'est desribed sur la page 11 de la présentation que vous avez lié (x est la multiplication). Utilisez pi = 3,1415 si la valeur exacte n'est pas intégrée.
3. Exécuter la régression multiple linéaire avec Valeur comme résultat et SINTIME et CosTime comme deux facteurs prédictifs. Vous devriez obtenir des estimations de leurs coefficients, que nous appellerons A et B .
4. Le terme d'interception du modèle de régression est le mésor.
5. l'amplitude est sqrt ( A ^ 2 + B ^ 2) [racine carrée de A carré plus B au carré]
6. Le acrophase est arctan (- B / A ), où arctan est la fonction inverse de tan. Les deux dernières formules sont de p.14 de la présentation.
7. Le modèle de régression devrait également vous donner une valeur R-carré pour voir comment le modèle circadien de 24 heures correspond aux données, et une p-valeur globale que les tests pour détecter la présence d'un composant circadien avec période de 24 heures.
8. On peut obtenir des erreurs standard sur l'amplitude et phase en utilisant des formules d'erreur standard Propogation, mais qui ne figure pas dans la présentation.