使用r中的optim控制参数的位数

Question

我希望减少此示例中预测参数的位数。有什么想法吗？

f <- function(par) {  
  a <- par[1]  
  b <- par[2]  
  (a - 1) + 3.2/b + 3*log(gamma(a)) + 3*a*log(b);  
}     

optim(c(1, .3), fn=f)  
# $par  
# [1] 1.399685 0.762025