r - R で polyroot の行列の各行を適用する

Question

3 次多項式の係数のさまざまな値を含む行列mがあり、これらの 3 次関数の最大実根を見つけたいと考えています。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")

行列の各行に polyroot を使用したいのですがm、R でこれを行うにはどうすればよいですか?

Answer 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