我能';I don’我不能让我喜欢的函数在R中工作

我能';I don’我不能让我喜欢的函数在R中工作,r,R,我在R中编程该似然函数时遇到问题: 我已经设法对函数h和lambda进行了编码。这是我需要计算的总数,到目前为止,我的尝试如下: log.lik<-function(x){ sigma<-x[1]; tau<-x[2] for(j in (1:6)){ t<-rep(0,7) t[j]<- -0.5((I-j)*log((1/p[j])*sigma)-log(1+lambda1(j,sigma,tau)*sum1[j])-h_j.1(j,s

我在R中编程该似然函数时遇到问题:

我已经设法对函数
h
lambda
进行了编码。这是我需要计算的总数,到目前为止,我的尝试如下:

log.lik<-function(x){
  sigma<-x[1]; tau<-x[2]
  for(j in (1:6)){
    t<-rep(0,7)
    t[j]<- -0.5((I-j)*log((1/p[j])*sigma)-log(1+lambda1(j,sigma,tau)*sum1[j])-h_j.1(j,sigma,tau)
}
return(-0.5*sum(t))
}

log.lik在不知道自己实际问题的情况下,我猜你是想改变立场

t<-rep(0,7)

t什么类型的对象是
x
?x是一个包含两个参数的向量:xClose vote,因为没有定义
h_j.1
lambda
或任何其他向量。hmm实际上给出了更好的结果,至少我的hessian矩阵现在有值了。只是不确定我的参数是否正确,但谢谢!