求二元正态分布的极大似然估计

风之伤

[img]file:///C:\Users\li\AppData\Roaming\Tencent\Users\306253398\QQ\WinTemp\RichOle\1IXE_7B[2JTX2A0_97(W7R7.jpg[/img][img=0,1]file:///C:\Users\li\AppData\Roaming\Tencent\Users\306253398\QQ\WinTemp\RichOle\1IXE_7B[2JTX2A0_97(W7R7.jpg[/img]
请问 这个函数对吗？
L<-function(theta){
  mu<-c(theta[1],theta[2]);
  p<-theta[5];
  sigma<-matrix(c(theta[3],theta[5],theta[5],theta[4]),ncol=2);
  f<-function(x){t(x-mu)%*%solve(sigma)%*%(x-mu)};
  L=1/((2*pi)^n*p/2*det(sigma)^n/2)*exp(-0.5*sum(apply(X,1,f)))
}

   如果没问题 为什么我optim不出来呢？