ผกผันของเมทริกซ์ใน R

ฉันสงสัยว่าคุณแนะนำวิธีใดในการคำนวณผกผันของเมทริกซ์

วิธีที่ฉันพบดูเหมือนจะไม่น่าพอใจ ตัวอย่างเช่น,

> c=rbind(c(1, -1/4), c(-1/4, 1))  
> c  
      [,1]  [,2]  
[1,]  1.00 -0.25  
[2,] -0.25  1.00  
> inv(c)  
Error: could not find function "inv"  
> solve(c)    
          [,1]      [,2]  
[1,] 1.0666667 0.2666667  
[2,] 0.2666667 1.0666667  
> solve(c)*c  
            [,1]        [,2]  
[1,]  1.06666667 -0.06666667  
[2,] -0.06666667  1.06666667  
> qr.solve(c)*c  
            [,1]        [,2]  
[1,]  1.06666667 -0.06666667  
[2,] -0.06666667  1.06666667  

ขอบคุณ!

solve(c)ให้ผกผันที่ถูกต้อง ปัญหาเกี่ยวกับรหัสของคุณคือคุณใช้ตัวดำเนินการผิดสำหรับการคูณเมทริกซ์ คุณควรใช้solve(c) %*% cเพื่อเรียกใช้การคูณเมทริกซ์ใน R

องค์ประกอบ R solve(c) * cดำเนินการโดยคูณองค์ประกอบเมื่อคุณเรียก

คุณสามารถใช้ฟังก์ชันginv () (Moore-Penrose generalized inverse) ในแพ็คเกจMASS

โปรดทราบว่าหากคุณสนใจความเร็วและไม่จำเป็นต้องกังวลเกี่ยวกับความเป็นเอกฐานsolve()ควรเลือกใช้ginv()เพราะเร็วกว่ามากเนื่องจากคุณสามารถตรวจสอบได้ดังนี้

require(MASS)
mat <- matrix(rnorm(1e6),nrow=1e3,ncol=1e3)

t0 <- proc.time()
inv0 <- ginv(mat)
proc.time() - t0 

t1 <- proc.time()
inv1 <- solve(mat)
proc.time() - t1