analisi statistica Cosinor,

Question

Ciao sto provando a calcolare un'analisi Cosinor in statistica, ma sono in perdita su come farlo. Ho bisogno di calcolare il MESOR, ampiezza e acrofase dei dati ritmo ciracadian.

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

v'è un link che mostra come lo fanno, le formule e simili, ma non mi ha dato molto aiuto. Qualcuno sa il codice per esso, sia in statistica o SPSS ??

ho davvero bisogno di arrivare a questo fatto perché è per un documento importante

Answer 1

Non ho SPSS o Statistica, quindi non posso dirvi il "push-questo-tasto" tipo esatto di passi, ma forse questo ti aiuterà.

analisi Cosinor conviene una curva coseno (o sinusoidale) con un periodo noto. L'idea principale è che il problema non lineare di montare una funzione coseno può essere ridotto ad un problema che è lineare nei suoi parametri se si conosce il periodo. I assumerà che il periodo T = 24 ore.

1. Si dovrebbe già avere due variabili: Ora in cui si esegue la misura, e Valore della misurazione (questi, naturalmente, potrebbe essere chiamato qualcos'altro).
2. Ora create due nuove variabili: SinTime = sin (2 x pi x Ora / 24) e CosTime = cos (2 x pi x Ora / 24) - questo è desribed sulla p.11 della presentazione si è collegato (x è la moltiplicazione). Usa pi = 3.1415 se il valore esatto non è built-in.
3. Esegui regressione lineare multipla con Valore come risultato e SinTime e CosTime come due predittori. Si dovrebbe ottenere stime dei loro coefficienti, che chiameremo A e B .
4. Il termine di intercetta del modello di regressione è la MESOR.
5. l'ampiezza è sqrt ( A ^ 2 + B ^ 2) [radice quadrata di A quadrato più B quadrato]
6. L'acrofase è arctan (- B / A ), dove arctan è la funzione inversa di abbronzatura. Le ultime due formule sono da p.14 della presentazione.
7. Il modello di regressione dovrebbe anche dare un valore di R-quadro per vedere come il ritmo circadiano 24 ore si adatta ai dati, e un p-value complessivo che i test per la presenza di un componente circadiano con un periodo di 24 ore.
8. Si può ottenere errori standard in ampiezza e fase utilizzando formule errore di propagazione standard, ma che non è incluso nella presentazione.