logsumexp

log sum exp

メモ 指数関数を使うと値が小さくなりすぎるものの足し合わせの際に、複数の値の相対値がそれほど小さくならないことを利用して、計算の精度を担保する方法 f <- function(a,b){ #log(exp(a)+exp(b)) b + log(exp(a-b)+1) } g <- function(as){ if(length(as…