For loop for Covariance of two vectors

16 visualizzazioni (ultimi 30 giorni)
I am looking to make a for loop to identify the individual covariance between each point in two lines and plot them as a color map along a line
To make the for loop I have:
x = rand(1, 5521)
y = rand(1, 5521)
for ii = 1:length(x)
covar(ii) = cov(x,y)
end
I understand that I need to make an empty array to enter the new data into and I understand that I am not calling the covariance to do individual points of the vector, but I am not sure how to fix it.
The plot I want to plot the two variables x and y ontop of eachother with a color map on the bottom showing the differences of variance along the line.

Risposta accettata

Walter Roberson
Walter Roberson il 13 Set 2022
Individual covariances between two points is always [0 0; 0 0] because there is a constant relationship between the two scalars. Knowing the one possible input scalar tells you the only possible output for that highly restricted population.
You need at least two x and y in order to get meaningful covariances.
Also remember that covariance between two classes (x, y) is a 2 x 2 array, but you want to color the result, which implies you only want one value per x y rather than 4 values.

Più risposte (1)

Torsten
Torsten il 13 Set 2022
Modificato: Torsten il 13 Set 2022
The covariance between two given vectors is a scalar, namely
cov(x,y) = 1/(n-1) * sum_{i=1}^{i=n} (x(i)-mean(x))*(y(i)-mean(y)):
rng('default')
format long
x = rand(1, 5521);
y = rand(1, 5521);
covariance = 1/5520 * sum((x-mean(x)).*(y-mean(y)))
covariance =
-9.473031297795489e-04
or directly with the matlab function
covariance_matrix = cov(x,y);
covariance = covariance_matrix(2,1)
covariance =
-9.473031297795485e-04
Calculating the covariance between individual data points does not make sense.
Calculating the covariance between random variables is something different.
So I'm not sure what you are trying to do in your code.

Prodotti


Release

R2020a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by