行列の指数関数というのがある 簡単に言うと、固有値を指数の肩に乗せて、それを固有ベクトルの行列でサンドイッチする exp.m <- function(A,n=1,eigen=FALSE){ # 固有値分解 eigen.out<-eigen(A) # P=V,P^{-1}=U V<-eigen.out[[2]] U<-solve(V) B<-diag(exp…

単体の体積の計算は、行列式の計算で一発 こちら SimplexVolume<-function(x,Factorial=FALSE){ n<-length(x[,1]) #d<-t(x[2:n,])-x[1,] d<-apply(x,2,FUN="diff") if(Factorial){ ret<-log(abs(det(d))) - lfactorial(n-1) }else{ ret<-log(abs(det(d))) } …