statistics - 统计余弦分析，

Question

嘿，我正在尝试在 statistica 中计算余弦分析，但不知道如何计算。我需要计算昼夜节律数据的 MESOR、AMPLITUDE 和 ACROPHASE。

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

有一个链接显示了如何做到这一点，公式等，但它并没有给我太多帮助。有谁知道它的代码，无论是在 statistica 还是 SPSS 中？

我真的需要完成这项工作，因为这是一篇重要的论文

Answer 1

我没有 SPSS 或 Statistica，所以我无法告诉您确切的“按此按钮”类型的步骤，但也许这会有所帮助。

余弦分析是拟合具有已知周期的余弦（或正弦）曲线。主要思想是，如果周期已知，拟合余弦函数的非线性问题可以简化为参数线性问题。我会假设你的时间 T=24 小时。

1. 您应该已经有两个变量：进行测量的时间和测量值（当然，这些变量可能被称为其他变量）。
2. 现在创建两个新变量：SinTime = sin(2 x pi x Time / 24) 和CosTime = cos(2 x pi x Time / 24) - 这在您链接的演示文稿的第 11 页有所描述（x 是乘法）。如果没有内置精确值，请使用 pi=3.1415。
3. 以Value作为结果，SinTime和CosTime作为两个预测变量运行多元线性回归。你应该得到它们系数的估计值，我们称之为A和B。
4. 回归模型的截距项是 MESOR。
5. AMPLITUDE 是 sqrt( A ^2 + B ^2) [ A 的平方根加上B的平方]
6. ACROPHASE 是 arctan(- B / A )，其中 arctan 是 tan 的反函数。最后两个公式来自演示文稿的第 14 页。
7. 回归模型还应该为您提供一个 R 平方值，以了解 24 小时昼夜节律模式与数据的拟合程度，以及一个整体 p 值，用于测试是否存在 24 小时周期的昼夜节律成分。
8. 可以使用标准误差传播公式获得幅度和相位的标准误差，但这不包括在演示文稿中。