wektor korelacji współczynników

Procesy stochastyczne. Sposoby racjonalizowania wielkich ilości informacji. Matematyka w naukach społecznych.
magda2291
Użytkownik
Użytkownik
Posty: 82
Rejestracja: 31 maja 2009, o 14:37
Płeć: Kobieta
Podziękował: 30 razy

wektor korelacji współczynników

Post autor: magda2291 »

Ogromnie prosze o pomoc, nie mam pojecia jak utworzyc wektor z tej macierzy:
W wyniku przeprowadzonych sześciu pomiarów doświadczenia otrzymano następujace dane:
\(\displaystyle{ \begin{tabular}{ccccccc}
y_{t} & x_{1t} & x_{2t} & x_{3t} & x_{4t}\\
20 & 1 & 5 & 18 & 5\\
21 & 0 & 10 & 16 & 2\\
24 & 0 & 12 & 15 & 6\\
26 & 1 & 18 & 12 & 1\\
32 & 0 & 20 & 8 & 0\\
39 & 1 & 25 & 3 & 4\\
\end{tabular}}\)

Wyznaczyć wektor oraz macierz współczynników korelacji między zmiennymi \(\displaystyle{ y}\) , \(\displaystyle{ x_{1}}\),\(\displaystyle{ x_{2}}\),\(\displaystyle{ x_{3}}\),\(\displaystyle{ x_{4}}\)
szw1710

wektor korelacji współczynników

Post autor: szw1710 »

Wektor ma zawierać współczynniki korelacji pomiędzy \(\displaystyle{ y}\) a kolejno zmiennymi \(\displaystyle{ x_1,x_2,x_3,x_4}\). Macierz ma zawierać liczby \(\displaystyle{ r_{ij}}\) będące współczynnikami korelacji pomiędzy zmiennymi \(\displaystyle{ x_i}\) a \(\displaystyle{ x_j}\).

Uruchom w R następujący skrypt:

Kod: Zaznacz cały

y=c(20,21,24,26,32,39)
x1=c(1,0,0,1,0,1)
x2=c(5,10,12,18,20,25)
x3=c(18,16,15,12,8,3)
x4=c(5,2,6,1,0,4)
Wektor=c(cor(y,x1),cor(y,x2),cor(y,x3),cor(y,x4))
Macierz=matrix(
 c(
   cor(x1,x1),cor(x1,x2),cor(x1,x3),cor(x1,x4),
   cor(x2,x1),cor(x2,x2),cor(x2,x3),cor(x2,x4),
   cor(x3,x1),cor(x3,x2),cor(x3,x3),cor(x3,x4),
   cor(x4,x1),cor(x4,x2),cor(x4,x3),cor(x4,x4)
  ),
nrow=4,ncol=4)
Wektor
Macierz
Oto wyniki jego działania:

Kod: Zaznacz cały

> Wektor
[1]  0.2010076  0.9436357 -0.9939523 -0.2093589
> Macierz
           [,1]       [,2]       [,3]       [,4]
[1,]  1.0000000  0.1496264 -0.1948709  0.1543033
[2,]  0.1496264  1.0000000 -0.9670683 -0.4040375
[3,] -0.1948709 -0.9670683  1.0000000  0.3006924
[4,]  0.1543033 -0.4040375  0.3006924  1.0000000
> 
Dopiero uczę się R, więc wszelkie sugestie co do lepszego napisania tego skryptu są mile widziane.
ODPOWIEDZ