fade
的作者似乎已将lambda
暴露为可以指定的参数。寻找与showMethods(fade, includeDef=TRUE)
代码,一切都在rEMM:::.fade
,在这里我们看到
> rEMM:::.fade
function (x, t = 1, lambda = NULL)
{
if (is.null(lambda))
lambda_factor <- [email protected]_d$lambda_factor
else lambda_factor <- 2^(-lambda)
[email protected]_d$counts <- [email protected]_d$counts * lambda_factor^t
[email protected]_d$mm <- smc_fade([email protected]_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))
您是否希望修改程序包本身,或者您是否愿意拥有一个“fade”函数来执行您指定的操作? – 2012-07-27 04:05:25
我真的不想修改这个包,但可能只是在这个函数上有一个包装。 – notrockstar 2012-07-27 16:28:51