ソースは、非常に重いのだが… 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))) } return(ret) } Simp…
引用をストックしました
引用するにはまずログインしてください
引用をストックできませんでした。再度お試しください
限定公開記事のため引用できません。