1

我正在尝试优化一个涉及伽马函数的函数。我的数据是经过审查的数据。我遇到的 错误是:“fn(par,...)中的错误:尝试应用非功能” R代码是:

logB<-function(theta, data){
#theta=(lambda, rho, tau)
#hpa = the data with first column X, second column Age, third column t,fourth column group and fifth column delta
data<-hpa
n<-length(data[,1])
t<-data[,2]
ep<-(t<=theta[3])
delta<-data[,4]
D<-sum(delta*ep)
d2<-sum(ep)
d3<-sum(delta)
w1<-sum(delta*ep*t+(1-delta)*ep*t)
w2<-sum(delta*(1-ep)*t)+ sum((1-delta)*(1-ep)*t)
d<-(d3-d1)
b<-((lgamma(d)+ lgamma(D))- d*log(w2-theta[3](n-d2))- D*log(w1+theta[3](n-d2)))
return(b)
  }
optim(c(0.4,0.9,96),logB, method="Nelder-Mead")

fn(par, ...) 中的错误:尝试应用非函数

预先感谢您的帮助。

4

1 回答 1

3

d*log(w2-theta[3](n-d2))- D*log(w1+theta[3](n-d2))

你可能想要

d*log(w2-theta[3]*(n-d2))- D*log(w1+theta[3]*(n-d2))

即,您缺少乘法运算符。

于 2013-07-07T18:13:38.570 回答