albert.dorador en estudiant.upc.edu
2016-Feb-09 09:16 UTC
[R-es] Paquete para elevar un matriz a una potencia real
Saludos a todos! Conoceis algun paquete que permita elevar una matriz a cualquier potencia, sin la limitacion de que la potencia haya de ser un entero? Se que existe el paquete expm que implementa el operando %^%, pero la potencia ha de ser necesariamente un entero (i.e. no permite elevar a decimales, fracciones, etc.). Agradeceria enormemente saberlo ya que estoy ultimando los detalles de una funcion que he programado que permite elevar cualquier matriz (cuadrada) a cualquier potencia real (incluyendo fracciones, etc.), y de no existir ninguna funcion similar trataria de incluirla en algun paquete relacionado o crear uno propio. Muchas gracias de antemano!
Olivier Nuñez
2016-Feb-09 11:15 UTC
[R-es] Paquete para elevar un matriz a una potencia real
require(expm) M <- matrix(c(1,0,1,1),2,2) M [,1] [,2] [1,] 1 1 [2,] 0 1 sqrt.M=expm(.5*logm(M)) sqrt.M%*%sqrt.M [,1] [,2] [1,] 1 1 [2,] 0 1 Un saludo. Olivier ----- Mensaje original ----- De: "albert dorador" <albert.dorador en estudiant.upc.edu> Para: r-help-es en r-project.org Enviados: Martes, 9 de Febrero 2016 10:16:08 Asunto: [R-es] Paquete para elevar un matriz a una potencia real Saludos a todos! Conoceis algun paquete que permita elevar una matriz a cualquier potencia, sin la limitacion de que la potencia haya de ser un entero? Se que existe el paquete expm que implementa el operando %^%, pero la potencia ha de ser necesariamente un entero (i.e. no permite elevar a decimales, fracciones, etc.). Agradeceria enormemente saberlo ya que estoy ultimando los detalles de una funcion que he programado que permite elevar cualquier matriz (cuadrada) a cualquier potencia real (incluyendo fracciones, etc.), y de no existir ninguna funcion similar trataria de incluirla en algun paquete relacionado o crear uno propio. Muchas gracias de antemano! _______________________________________________ R-help-es mailing list R-help-es en r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es