2013-10-10 3 views
0

에 적용합니다. 입방 다항식의 계수 값이 다른 amatrix m이 있습니다.이 큐빅 함수의 최대 실수 근원을 찾고 싶습니다. R에서 이것을 어떻게 할 수 있습니까?다원 행렬에 대한 행렬의 각 행을 R

 m = matrix(c(1, 3, 4, 6, 7, 8, 9, 2, -3, -5, 6, 2), nrow = 3) 
     colnames(m) <- c("d", "c", "b", "a") 
나는 m 행렬의 각 행에 대해 polyroot을 사용하고 싶습니다

, 어떻게이 R에서 할 수 있습니까?

답변

1

이 불과하지 않을까요 :

apply(m , 1, polyroot) 

... 결과는 열 순서에있을 것입니다 기억하는, 물론있다.

apply(apply(m, 1,polyroot) , 2, function(x) ifelse(abs(Im(x)) < 0.00001, x,NA)) 
      [,1]   [,2]   [,3] 
[1,]   NA -0.416256-0i -0.4167669+0i 
[2,]   NA   NA   NA 
[3,] 2.34748+0i   NA   NA 
+0

감사합니다.하지만 어떻게 최대 실제 루트를 얻을 수 있습니까? – rose

+0

글쎄, 그건 원래의 질문이 아니 었지, 그렇지? –

+0

아니, soory, 나는 그것을 언급하는 것을 잊었다. – rose

관련 문제