Làm thế nào để làm phép tính ma trận trong R

Ngoài các toán tử số học cổ điển, R chứa một bộ các toán tử và các hàm lớn để thực hiện một tập hợp rộng các hoạt động ma trận. Nhiều trong số những hoạt động này được sử dụng trong toán học tiên tiến, do đó bạn có thể không bao giờ cần đến chúng. Một số trong số chúng có thể có khá tiện dụng, mặc dù, nếu bạn cần phải lật xung quanh dữ liệu hoặc bạn muốn tự mình tính một số thống kê.

Chuyển một ma trận trong R

Xoay quanh ma trận để các hàng trở thành các cột và ngược lại rất dễ dàng trong R. Chức năng t ()

>> t (first matrix) [1] [2] [3] [1,] 1 2 3 [2,] 4 5 6 [3,] 7 8 9 [4,] 10 11 12

Bạn có thể thử điều này với một vector, quá. Khi các ma trận được đọc và điền vào cột, không nên ngạc nhiên khi thấy hàm t () xem vector như một ma trận một cột. Sự chuyển vị của một vector là một ma trận một hàng: >> t (1: 10) [1] [2] [3] [4] [5] [6], 7] [8] [9] [10] [1,] 1 2 3 4 5 6 7 8 9 10

Bạn có thể nói đây là một ma trận theo kích thước. Thông tin này có vẻ tầm thường nhưng bạn hãy tưởng tượng bạn đang chỉ chọn một hàng từ ma trận và chuyển nó. Không giống như những gì bạn mong đợi, bạn sẽ nhận được một hàng thay vì một cột:

 >> t (first matrix [2,]) [1] [2] [3] [4] [1,] 2 5 8 11

Đảo ngược ma trận trong R

Trái ngược với trực giác của bạn, đảo ngược một ma trận không được thực hiện bằng cách nâng nó lên sức mạnh của -1, R thường áp dụng các toán tử số học yếu tố trên ma trận. Vì vậy, lệnh đầu tiên. ma trận ^ (- 1) không cung cấp cho bạn các nghịch đảo của ma trận; thay vào đó, nó mang lại cho bạn các nghịch đảo của các yếu tố. Để đảo ngược một ma trận, bạn sử dụng hàm giải quyết () như sau:

 >> vuông. ma trận giải (ma trận vuông) [1] [2] [3] [1,] 0. 5 -0. 8333333 0,1666667 [2,] -0. 5 0.6666667 0. 1666667 [3,] 0. 5 -0. 1666667 -0. 1666667

Hãy cẩn thận đảo ngược một ma trận như thế này do nguy cơ lỗi tròn. R tính toán hầu hết các thống kê dựa trên các phân rã như phân hủy QR, phân hủy đơn giá trị và sự phân hủy Cholesky. Bạn có thể làm điều đó bằng cách sử dụng các chức năng qr (), svd () và chol (), tương ứng. Kiểm tra các trang trợ giúp tương ứng để biết thêm thông tin.

Nhân hai ma trận trong R

Toán tử nhân (*) làm việc theo nguyên tố trên các ma trận. Để tính sản phẩm bên trong của hai ma trận, bạn sử dụng toán tử đặc biệt% *%, như sau:

đầu tiên. Ma trận% *% t (ma trận thứ hai) [1] [2] [3] [1] 22 44 66 [2,] 26 52 78 [3,] 30 60 90

Bạn phải chuyển đổi thư hai.ma trận đầu; nếu không, cả hai ma trận có kích thước không phù hợp. Nhân một ma trận với một vector là một chút của một trường hợp đặc biệt; miễn là các kích thước phù hợp, R sẽ tự động chuyển đổi các vector vào một hàng hoặc một ma trận cột, bất cứ điều gì có thể áp dụng trong trường hợp đó. Bạn có thể kiểm tra cho mình trong ví dụ sau: >> đầu tiên. ma trận% *% 1: 4 [1] [1,] 70 [2,] 80 [3,] 90> 1: 3% *% đầu tiên. ma trận [1] [2] [3] [4] [1,] 14 32 50 68