gamma(300，1.5，0.5)<br>
求救，关于伽马函数的用极大似然，然后关于两个参数求偏导等于0，然后解出来，如何用solve函数<br>
关于这个步骤如何用solve代码写出来，求求各位能不能写一下，下面是我这个初学者写的代码，也不太对<br>
library(BB)<br>
function(x)  { <br>
  f<deriv(f(theta), c("alpha", "beta"), func = TRUE)<br>
  nu0 <-0<br>
  tau0 <-0<br>
   f[1] <-n*alpha*1/beta-sum(dato)-nu0<br>
  f[2] <-n*log(beta)-n*1/gamma(alpha)*gamma(alpha)-tau0<br>
  f<br>
}<br>
p0 <- c(1.5, 0.5)          <br>
BBsolve(par=p0, fn=func)