我正在使用正态分布(概率)对二进制独立变量进行对数似然函数编程。使用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)