статистический анализ косинора,

Question

Эй, я пытаюсь вычислить анализ косиноров в статистике, но я не понимаю, как это сделать. Мне нужно вычислить MESOR, AMPLITUDE и ACROPHASE данных о циракадном ритме.

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

есть ссылка, которая показывает, как это сделать, формулы и такие, но это не дало мне большую помощь. Кто-нибудь знает код для него, либо в статистике, либо в SPSS?

мне действительно нужно, чтобы это было сделано потому, что это для важной работы

Answer 1

у меня нет SPSS или Statistica, поэтому я не могу сказать вам точно «PUSH-это кнопки» подобные шаги , но, возможно, это поможет.

Анализ косиноров подходит для кривой косинуса (или синуса) с известным периодом. Основная идея заключается в том, что нелинейная задача подгонки косинусной функции сводится к задаче, линейной по своим параметрам, если известен период. Я буду считать, что ваш период T = 24 часа.

1. Вы уже должны иметь две переменные: Time, при которой берется измерение, и Value измерения (это, конечно, можно было бы назвать что-то еще).
2. Теперь создайте две новые переменные: SinTime = sin (2 х пи х Время/24) и CosTime = Cos (2 х пи х Время/24) - это на стр.11 мы проделали из представленной вами презентации (x - умножение). Используйте pi = 3.1415, если точное значение не встроено.
3. Run множественная линейная регрессия с Value в качестве результата и SinTime и CosTime как два предикторов. Вы должны получить оценки своих коэффициентов, которые мы будем называть A и B.
4. Термин перехвата регрессионной модели - MESOR.
5. амплитудно является SQRT (^2 + В^2) [квадратный корень из квадрат плюс В квадрат]
6. акрофаза является арктангенс (- Б/), где arctan - обратная функция tan. Последние две формулы взяты из п. 14 презентации.
7. Модель регрессии также должна дать вам значение R-квадрата, чтобы увидеть, насколько хорошо подходит 24-часовой циркадный образец, и общее значение p, которое проверяет наличие циркадного компонента с периодом 24 часа.
8. Можно получить стандартные ошибки по амплитуде и фазе с использованием стандартных формул пропозиции ошибок, но это не включено в презентацию.

Answer 2

Пожалуйста, обратите внимание, что существует потенциальная ошибка в противном случае, отличную реакцию - если А коэффициент регрессии для SinTime и B является коэффициент CosTime, чем акрофаза является арктангенс (-A/B), а не как подразумевалось.