covariance matrix

Question

Мне было интересно, может ли кто-нибудь объяснить мне, как пакет geoR вычисляет функцию ковариации? Я имею в виду, как вы это сделаете вручную?

library(geoR) 
#suppose I have the following coordinates 
X = c(60,30,20,40) 
Y = c(50,20,50,50) 
my_coordinates = cbind(X,Y) 
print(my_coordinates) 

#computing covariance 
my_cov= varcov.spatial(my_coordinates,cov.model="exp", cov.pars=c(0.2,25)) 
print(my_cov) 

И вы получите:

  [,1]  [,2]  [,3]  [,4] 
[1,] 0.20000000 0.03664442 0.04037930 0.08986579 
[2,] 0.03664442 0.20000000 0.05645288 0.05645288 
[3,] 0.04037930 0.05645288 0.20000000 0.08986579 
[4,] 0.08986579 0.05645288 0.08986579 0.20000000 

Однако, один, возможно, захотите сделать это в Matlab, а также.

Answer 1

Лучший способ узнать, как пакет или функция что-то делает, чтобы посмотреть исходный код. Это одна из потрясающих особенностей проектов с открытым исходным кодом, вы можете это сделать.

попробуй ввести varcov.spatial или поиск через распакованный пакет деготь шар для определения функции

Для вычисления ковариации (которое зависит от расстояния между точками), вам нужно вычислить

• на расстояние между вашими точками (вам действительно нужен только нижний треугольник, так как он будет симметричным
• Значение функции ковариации на каждом расстоянии
• сформировать полную симметричную v матрицы ковариации ариесности из этих вычисленных ковариаций.

Ковариационные функции определены в ?cov.spatial. Вы можете позвонить по телефону cov.spatial, чтобы рассчитать это значение R (что именно такое geoR::varcov.spatial)