Biện pháp phù hợp để tìm ma trận hiệp phương sai nhỏ nhất


10

Trong sách giáo khoa tôi đang đọc họ sử dụng tính xác định dương (xác định bán tích cực) để so sánh hai ma trận hiệp phương sai. Ý tưởng được rằng nếu là PD sau đó là nhỏ hơn . Nhưng tôi đang đấu tranh để có được trực giác của mối quan hệ này?ABBA

Có một chủ đề tương tự ở đây:

/math/239166/what-is-the-intuition-for-USE-definitiness-to-compare-matrices

Trực giác cho việc sử dụng dứt khoát để so sánh ma trận là gì?

Mặc dù câu trả lời rất hay nhưng họ không thực sự đề cập đến trực giác.

Đây là một ví dụ tôi thấy khó hiểu:

[1612129][1224]

Bây giờ ở đây xác định của sự khác biệt là -25 vì vậy mối quan hệ không phải là pd hoặc thậm chí psd và vì vậy ma trận đầu tiên không lớn hơn đầu tiên?

Tôi chỉ đơn giản muốn so sánh hai ma trận hiệp phương sai 3 * 3 để xem cái nào nhỏ nhất? Nó có vẻ trực quan hơn đối với tôi để sử dụng một cái gì đó giống như tiêu chuẩn euclide để so sánh chúng? Tuy nhiên, điều này có nghĩa là ma trận đầu tiên ở trên lớn hơn ma trận thứ hai. Hơn nữa, tôi chỉ thấy tiêu chí pd / psd được sử dụng để so sánh ma trận hiệp phương sai.

Ai đó có thể giải thích tại sao pd / psd tốt hơn so với việc sử dụng một biện pháp khác như định mức euclide?

Tôi cũng đã đăng câu hỏi này lên diễn đàn toán học (không chắc điều gì là tốt nhất) hy vọng điều này không trái với bất kỳ quy tắc nào.

/math/628135/compared-two-covariance-matrices


2
Bạn có thể muốn đọc điều này trong đó trực giác đằng sau sự dứt khoát (bán) tích cực được xem xét. Khi bạn so sánh 2 phương sai ab, nếu a-blà dương thì chúng ta sẽ nói rằng khi loại bỏ biến thiên bra khỏi ađó vẫn còn một số biến thiên "thực" còn lại a. Tương tự như vậy là một trường hợp phương sai đa biến (= ma trận hiệp phương sai) AB. Nếu A-Blà xác định dương thì điều đó có nghĩa là A-Bcấu hình của vectơ là "thực" trong không gian euclide: nói cách khác, khi loại bỏ Bkhỏi A, cái sau vẫn là một biến thiên khả thi.
ttnphns

2
Làm gì bạn có nghĩa là do "nhỏ nhất" của hai ma trận hiệp phương sai?
whuber

Xin chào, các ma trận hiệp phương sai liên quan đến các công cụ ước tính cạnh tranh, tôi muốn chọn công cụ ước tính có phương sai nhỏ nhất. (Điều này có làm rõ mọi thứ không?)
Baz

2
Baz: Tại sao không so sánh trực tiếp phương sai của các công cụ ước tính?
Glen_b -Reinstate Monica

Xin chào, phương thức được đặt, biểu thức cho cái mà họ gọi là phương sai (bao gồm hiệp phương sai) được đưa ra. Tuy nhiên, ngay cả khi tôi chỉ so sánh các phương sai, điều này vẫn liên quan đến việc so sánh các giá trị vectơ sẽ có vấn đề tương tự như so sánh các giá trị ma trận?
Baz

Câu trả lời:


8

Thứ tự các ma trận mà bạn đề cập được gọi là thứ tự Loewner và là một thứ tự một phần được sử dụng nhiều trong nghiên cứu các ma trận xác định dương. Một điều trị độ dài cuốn sách của hình học trên đa tạp của ma trận dương-xác định (posdef) là ở đây .

Trước tiên tôi sẽ cố gắng giải quyết câu hỏi của bạn về trực giác . Ma trận (đối xứng) là posdef nếu cho tất cả . Nếu là biến ngẫu nhiên (rv) với ma trận hiệp phương sai , thì là (tỷ lệ với) phép chiếu của nó trên một không gian con một chiều và . Áp dụng điều này cho trong Q của bạn, thứ nhất: đó là ma trận hiệp phương sai, thứ hai: Một biến ngẫu nhiên với ma trận covar dự án theo mọi hướng có phương sai nhỏ hơn rv với ma trận hiệp phương saiAcTAc0cRnXAcTXVar(cTX)=cTAcABBA. Điều này làm cho trực giác rõ ràng rằng thứ tự này chỉ có thể là một phần, có nhiều rv sẽ chiếu theo các hướng khác nhau với các phương sai khác nhau cực kỳ. Đề xuất của bạn về một số chỉ tiêu Euclide không có cách giải thích thống kê tự nhiên như vậy.

"Ví dụ khó hiểu" của bạn là khó hiểu vì cả hai ma trận đều có số không xác định. Vì vậy, đối với mỗi người, có một hướng (eigenvector với eigenvalue zero) trong đó họ luôn chiếu về 0 . Nhưng hướng này là khác nhau cho hai ma trận, do đó chúng không thể được so sánh.

Thứ tự Loewner được định nghĩa sao cho , có giá trị xác định dương hơn , nếu là posdef. Đây là một thứ tự từng phần, đối với một số ma trận posdef, cả và đều không phải là posdef. Một ví dụ là: Một cách hiển thị đồ họa này là vẽ một đồ thị có hai hình elip, nhưng tập trung ở điểm gốc, được liên kết theo cách chuẩn với ma trận (khi đó khoảng cách xuyên tâm theo mỗi hướng tỷ lệ thuận với phương sai của chiếu theo hướng đó):ABBABABAAB

A=(10.50.51),B=(0.5001.5)

Hai ma trận posdef hiển thị dưới dạng elip

Trong những trường hợp này, hai hình elip đồng dạng, nhưng xoay khác nhau (thực tế là góc 45 độ). Điều này tương ứng với thực tế là các ma trận và có cùng giá trị riêng, nhưng các hàm riêng được quay.AB

Vì câu trả lời này phụ thuộc nhiều vào các thuộc tính của dấu chấm lửng, nên sau đây Trực giác đằng sau các phân phối Gaussian có điều kiện là gì? giải thích các hình elip về mặt hình học, có thể hữu ích.

Bây giờ tôi sẽ giải thích làm thế nào các hình elip liên quan đến ma trận được xác định. Ma trận posdef xác định dạng bậc hai . Điều này có thể được vẽ như là một hàm, đồ thị sẽ là một bậc hai. Nếu thì đồ thị của sẽ luôn ở trên đồ thị của . Nếu chúng ta cắt các biểu đồ với mặt phẳng nằm ngang ở độ cao 1, thì các vết cắt sẽ mô tả các hình elip (thực tế đó là một cách xác định các hình elip). Hình elip cắt này được cho bởi các phương trình và chúng ta thấy rằngAQA(c)=cTAcABQBQA

QA(c)=1,QB(c)=1
ABtương ứng với hình elip của B (bây giờ có phần bên trong) được chứa trong hình elip của A. Nếu không có thứ tự, sẽ không có ngăn chặn. Chúng tôi quan sát rằng thứ tự bao gồm ngược với thứ tự một phần Loewner, nếu chúng tôi không thích rằng chúng tôi có thể vẽ các hình elip của các nghịch đảo. Điều này là do tương đương với . Nhưng tôi sẽ ở lại với các dấu chấm lửng như được định nghĩa ở đây.ABB1A1

Một hình elip có thể được mô tả với dấu chấm phẩy và chiều dài của chúng. Chúng ta sẽ chỉ thảo luận về ma trận ở đây, vì chúng là những thứ chúng ta có thể vẽ ... Vì vậy, chúng ta cần hai trục chính và chiều dài của chúng. Điều này có thể được tìm thấy, như được giải thích ở đây với một sự xuất hiện của ma trận posdef. Sau đó, các trục chính được đưa ra bởi các hàm riêng và độ dài có thể được tính từ các giá trị riêng bởi Chúng ta cũng có thể thấy rằng diện tích của hình elip đại diện cho là .2×2một , b λ 1 , λ 2 một = a,bλ1,λ2

a=1/λ1,b=1/λ2.
Aπab=π1/λ11/λ2=πdetA

Tôi sẽ đưa ra một ví dụ cuối cùng trong đó các ma trận có thể được đặt hàng:

Hai ma trận có thể được đặt hàng được vẽ như hình elip

Hai ma trận trong trường hợp này là:

A=(2/31/51/53/4),B=(11/71/71)


3

@kjetil b halvorsen đưa ra một cuộc thảo luận tốt đẹp về trực giác hình học đằng sau tính bán chính xác tích cực như một thứ tự từng phần. Tôi sẽ đưa ra một cách dễ thương hơn với cùng một trực giác. Một trong số đó xuất phát từ loại tính toán nào bạn có thể muốn thực hiện với ma trận phương sai của mình.

Giả sử bạn có hai biến ngẫu nhiên và . Nếu chúng là vô hướng, thì chúng ta có thể tính toán phương sai của chúng là vô hướng và so sánh chúng theo cách rõ ràng bằng cách sử dụng các số thực vô hướng và . Vì vậy, nếu và , chúng ta nói rằng biến ngẫu nhiên có phương sai nhỏ hơn .xyV(x)V(y)V(x)=5V(y)=15xy

Mặt khác, nếu và là các biến ngẫu nhiên có giá trị véc tơ (giả sử chúng là hai vectơ), thì cách chúng ta so sánh phương sai của chúng không quá rõ ràng. Giả sử phương sai của chúng là: Làm thế nào để chúng ta so sánh phương sai của hai vectơ ngẫu nhiên này? Một điều chúng ta có thể làm chỉ là so sánh phương sai của các yếu tố tương ứng. Vì vậy, chúng ta có thể nói rằng phương sai của nhỏ hơn phương sai của bằng cách chỉ so sánh các số thực, như: vàxy

V(x)=[10.50.51]V(y)=[8336]
x1y1V(x1)=1<8=V(y1)V(x2)=1<6=V(y2). Vì vậy, có lẽ chúng ta có thể nói rằng phương sai của là phương sai của nếu phương sai của mỗi phần tử của là phương sai của phần tử tương ứng của . Điều này sẽ giống như nói nếu mỗi phần tử đường chéo của là phần tử đường chéo tương ứng của .xyxyV(x)V(y)V(x)V(y)

Định nghĩa này có vẻ hợp lý lúc đầu đỏ mặt. Hơn nữa, miễn là các ma trận phương sai mà chúng ta đang xem xét là đường chéo (tức là tất cả các hiệp phương sai là 0), nó cũng giống như sử dụng bán chính xác. Đó là, nếu các phương sai trông giống như sau đó nói là dương-bán xác định (nghĩa là ) giống như nói và . Tất cả có vẻ tốt cho đến khi chúng tôi giới thiệu hiệp phương sai. Xem xét ví dụ này:

V(x)=[V(x1)00V(x2)]V(y)=[V(y1)00V(y2)]
V(y)V(x)V(x)V(y)V(x1)V(y1)V(x2)V(y2)
V(x)=[10.10.11]V(y)=[1001]
Bây giờ, sử dụng phép so sánh chỉ xem xét các đường chéo, chúng ta sẽ nói và, thực sự, vẫn đúng là từng phần tử . Điều có thể bắt đầu làm phiền chúng tôi về điều này là nếu chúng tôi tính tổng các phần tử có trọng số của các vectơ, như và , thì chúng tôi gặp phải sự thật là mặc dù chúng tôi đang nói .V(x)V(y)V(xk)V(yk)3x1+2x23y1+2y2V(3x1+2x2)>V(3y1+2y2)V(x)V(y)

Điều này thật kỳ lạ đúng không? Khi và là vô hướng, thì đảm bảo rằng với mọi , cố định, không ngẫu nhiên .xyV(x)V(y)aV(ax)V(ay)

Nếu vì lý do nào đó, chúng tôi quan tâm đến sự kết hợp tuyến tính của các yếu tố của các biến ngẫu nhiên như thế này, thì chúng tôi có thể muốn củng cố định nghĩa của chúng tôi về cho các ma trận phương sai. Có lẽ chúng ta muốn nói khi và chỉ khi nó là sự thật rằng , không có vấn đề gì cố định số và chúng tôi chọn. Lưu ý, đây là định nghĩa mạnh hơn định nghĩa chỉ có đường chéo vì nếu chúng ta chọn thì nó nói và nếu chúng ta chọn thì nó nói .V(x)V(y)V(a1x1+a2x2)V(a1y1+a2y2)a1a2a1=1,a2=0V(x1)V(y1)a1=0,a2=1V(x2)V(y2)

Định nghĩa thứ hai này, định nghĩa khi và chỉ khi cho mọi vectơ cố định có thể , là phương pháp thông thường để so sánh phương sai ma trận dựa trên bán chính xác dương: Nhìn vào biểu thức cuối cùng và định nghĩa bán xác định dương để thấy rằng định nghĩa của cho ma trận phương sai được chọn chính xác để đảm bảo rằng khi và chỉ khi cho bất kỳ lựa chọn nào của , tức là khi là bán tích cực -xác định.V(x)V(y)V(ax)V(ay)a

V(ay)V(ax)=aV(x)aaV(y)a=a(V(x)V(y))a
V(x)V(y)V(ax)V(ay)a(V(y)V(x))

Vì vậy, câu trả lời cho câu hỏi của bạn là mọi người nói rằng ma trận phương sai nhỏ hơn ma trận phương sai nếu là bán xác định dương họ quan tâm đến việc so sánh phương sai của các tổ hợp tuyến tính của các vectơ ngẫu nhiên cơ bản. Định nghĩa nào bạn chọn tuân theo những gì bạn quan tâm khi tính toán và cách định nghĩa đó giúp bạn với những tính toán đó.VWWV

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.