आर

Question

में मैट्रिक्स के विपरीत मैं सोच रहा था कि मैट्रिक्स के विपरीत की गणना करने के लिए आपका अनुशंसित तरीका क्या है?

मुझे मिले तरीकों को संतोषजनक नहीं लगता है। उदाहरण के लिए,

> 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 

धन्यवाद!

Answer 1

solve(c) सही उलटा देता है। आपके कोड के साथ समस्या यह है कि आप मैट्रिक्स गुणा के लिए गलत ऑपरेटर का उपयोग कर रहे हैं। आपको solve(c) %*% c का उपयोग आर

आर पर गुणांक गुणा करने के लिए solve(c) %*% c का उपयोग करना चाहिए जब आप solve(c) * c का आह्वान करते हैं तो तत्व गुणा द्वारा तत्व निष्पादित करता है।

Answer 2

आप मास पैकेज में समारोह ginv() (मूर-Penrose सामान्यीकृत व्युत्क्रम) का उपयोग कर सकते

Answer 3

ध्यान दें कि यदि आप गति के बारे में परवाह है और विशिष्टता के बारे में चिंता करने की जरूरत नहीं है, 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 

Answer 4

मैट्रिक्स संकेतन में यह एक बड़ा अंतर ऑपरेटर "०१२३२२७५९२४१ बनाता है"और ऑपरेटर" %*% "। पहला तत्व द्वारा गुणा तत्व करता है, दूसरा मैट्रिक्स गुणा के लिए सही सूत्र है। क्या Hou किया जाना चाहिए है:

c = rbind(c(1, -1/4), c(-1/4, 1)) 

solve(c) %*% c