使用log likelihoood函数时,如何修复“ vmmin中的初始值不固定”?

时间:2019-04-27 21:37:45

标签: r

我正在使用正态分布(概率)对二进制独立变量进行对数似然函数编程。使用optim时,出现错误:

  

$ optim中的$ Error(起始值,probitll,gr = grad,Y = Ydoc,X =   Xbiais,:“ vmmin”中的初始值不是$$

probitll<-function( par, X, Y){ 
 Y<-as.matrix(Y)
 X<-as.matrix(X)
 K<-ncol(X)
 b<-matrix(1:K, ncol = 1)
 R<-as.vector(dnorm(X%*%b))
 -sum(Y*log(R)+(1-Y)*log(1-R))
}

grad<-function(SV,  X, Y){
  X<-as.matrix(X)
  Y<-as.matrix(Y)
  K<-ncol(X)
  b<-matrix(1:K, ncol = 1)
  R<-as.vector(dnorm(X%*%b))
  apply(R*X, 2, sum)
}

startvalues<-as.vector(modeltestMCO[,1])

resultprobit<-optim( startvalues, probitll, gr= grad,   Y=Ydoc, X=Xbiais, method="BFGS", hessian=TRUE)

0 个答案:

没有答案