の作者は指定できるパラメータとしてfade
公開しているようです。lambda
でコードを見るとshowMethods(fade, includeDef=TRUE)
、すべてが で終わりrEMM:::.fade
ます。
> rEMM:::.fade
function (x, t = 1, lambda = NULL)
{
if (is.null(lambda))
lambda_factor <- x@tnn_d$lambda_factor
else lambda_factor <- 2^(-lambda)
x@tnn_d$counts <- x@tnn_d$counts * lambda_factor^t
x@tracds_d$mm <- smc_fade(x@tracds_d$mm, lambda_factor^t)
invisible(x)
}
<environment: namespace:rEMM>
したがって、新しい関数を作成する代わりに、この関数で使用される変換と一致するlambda
ように変換して、好きな方法で計算してみませんか?lambda_factor
mylambda <- function(t, lambda, T, c)
-log(1/(1+(T-t/c)^lambda)), 2) / t
そして、次fade
のように呼び出します
fade(x, t, mylambda(t, lambda, T, c))
? これをさらに次のようにラップします。
myfade <- function(x, t, lambda)
fade(x, t, mylambda(t, lambda, T, c))