library(fAsianOptions) cgamma(1i) my.cbeta <- function(a,b){ ret <- cgamma(a)*cgamma(b)/cgamma(a+b) return(ret) } my.dcbeta <- function(x,a,b){ x^(a-1) * (1-x)^(b-1) / my.cbeta(a,b) } x <- seq(from=-1,to=2,length=100) a <- -0.5 b <- -0.5 y …

自然数について、は、偶数のときは自然数、奇数のときは、を含む リンク public static double GammaLN(int df){ double chiProb=0; int evenodd=(int)(df)%2; if(evenodd==0){ for(double i=(double)df/2-1;i>0;i--){ chiProb+=Math.log(i); } }else{ chiPr…