[R-es] Paquete para elevar un matriz a una potencia real

[Olivier Nuñez]

 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