Rank: 4537 based on 20 downloads (last 30 days) and 1 file submitted
photo

cheart

E-mail

Personal Profile:
Professional Interests:

 

Watch this Author's files

 

Files Posted by cheart
Updated   File Tags Downloads
(last 30 days)
Comments Rating
24 Jun 2008 Screenshot Cosinor Analysis Fits cosine curve to a time series using least squares Author: cheart statistics, probability, cosinor anlaysis cosi..., sinusoid, sine, harmonic 20 8
Comments and Ratings on cheart's Files View all
Updated File Comment by Comments Rating
10 Aug 2014 Cosinor Analysis Fits cosine curve to a time series using least squares Author: cheart Stanislav

It seems that unequal sampling problem can be resolved by addition ~isnan(y) to every value in NE matrix. Thus 2) is not an issue. :)

06 Aug 2014 Cosinor Analysis Fits cosine curve to a time series using least squares Author: cheart Stanislav

UPD: 3) added by mistake

06 Aug 2014 Cosinor Analysis Fits cosine curve to a time series using least squares Author: cheart Stanislav

I am by no means a specialist in Matlab, but there are several comments:
1) In section "Confidence interval for the mesor" is written 'sigma^2', whereas according to Nelson et al, 1979 and Bingham et al., 1982 it should be 'sigma'.
This leads to inflated estimation of mesor CI.
2) This code is can not be used for data sets with missing values - only equally spaced and full time series thus simplified estimation of M, beta and gamma can be applied as shown in Nelson et al.,1979.
3) According to Nelson et al., 1979 t should be computed for two-tailed distribution, thus 1-alpha/2 should be used.

20 Dec 2011 Cosinor Analysis Fits cosine curve to a time series using least squares Author: cheart Bulger, David

Wouldn't the p-value for the zero-amplitude test have to come from the F-distribution's CDF, rather than its PDF?

13 Mar 2009 Cosinor Analysis Fits cosine curve to a time series using least squares Author: cheart David

Undefined function or variable "phi".

Error in ==> cosinor at 94
fprintf(1,'Mesor = %g \nAmplitude = %g \nAcrophase = %g \n\n',M,Amp,phi);

Contact us