Search All of the Math Forum:
Views expressed in these public forums are not endorsed by
NCTM or The Math Forum.



Re: Non linear association
Posted:
Jan 20, 2013 1:24 PM


Hi Greg,
Thanks for the feedback  very helpful!
If I understand your suggestions for my third query correctly, you are using an approach similar to standard or absolute deviations to quantify the amount of variance in a signal? This could be useful for quantifying the consistency of signals; however I am looking to see if the patterns of change in two signals follow a similar style and timescale.
I have been playing with some code to calculate the first derivative of the signals which can then be analysed statistically. Would you say this is a valid approach?
Thanks again Dan
"Greg Heath" <heath@alumni.brown.edu> wrote in message <kdgp44$o93$1@newscl01ah.mathworks.com>... > "Daniel Robbins" <d.w.e.robbins@gmail.com> wrote in message <kdcenb$s4n$1@newscl01ah.mathworks.com>... > > Hi, > > > > I am trying to establish the best way to check the similarity of the pattern of change in two signals/vectors. For a simple example I generated four datasets. One positive parabolic curve, one equal negative curve, one smaller positive curve with a time lag and one vector with a fixed value e.g. [6 6 6 6 6 6 ...etc]. All vectors were the same lengths. > > > > My thinking is that I should use some sort of cross correlation e.g. > > > > [Rxy, Lag] = xcorr(data(:,4),data(:,3)); > > plot(Lag,Rxy) > > > > By changing the two signals compared I note that the peak of the graph changes location. Therefore I could use the x axis values to calculate the lag in the signal. However I don't understand the values on the Yaxis. I also note that if I input the fixed value vector the graph becomes platykurtic, I'm not sure how to interpret this. > > > > Please can someone help me with the following queries: > > > > 1.) How do I interpret the Yvalues of the graph? > > If you standardize those columns of data, then the two autocorrelation functions > will have unity peaks at zero lag. > > Then you can compare the peak of the crosscorrelation function at zero lag with unity > and how many lags it takes the function to decay to some fraction of the max. > > > 2.) Is there some way of quantifying the level of agreement in the form of liner correlations e.g. from 1 to 1? Or via pvalues? > > Most likely. However, that probably isn't the best approach. > > > 3.) Is this the best approach to analyse the pattern of change in two non linear signals or is there a better approach? > > I would choose sqrt(mean((y1y2).^2)), mean(abs(y1y2)) or max(abs(y1y2)) > > > Many thanks > > You're welcome. > > Greg



