下面这个式子反解alpha的值？或者怎么写出alpha的表达式啊

Lnlx_61

第一张图：Q函数对alpha求一阶偏导数，第二张图：将第一张图中后面的上不完全伽马函数对alpha的一阶偏导展开。 为求alpha的值，令Q函数对alpha的一阶偏导数为0，但是我不知道该怎么写出alpha的确切表达式，如alpha=***这样的。有没有大神救救孩子啊？ 能写出alpha的确切表达式或者用R语言能反解出alpha 以下个别部分的R语言程序：(nu_i是0-1变量)x=rgamma(20,4,3) tau=rgamma(20,6,4) #upgamma function   上gamma函数 upgamma <- function(shape, x) {   pgamma(x, shape=shape, lower=F)*gamma(shape) } #对表达式中出现的无穷级数求和的函数 s.1 = function(y,N,a,b){   s1=0   for(m in 0:N){     s1 = s1 + (y*b)^(a-m)/gamma(a-m+1)   }   return(s1) } s.2 = function(y,N,a,b){   s1=0   for(m in 0:N){     s1 = s1 + (y*b)^(a*-m)*digamma(a-m+1)/gamma(a-m+1)   }   return(s1) } ##作为独立功能的所需期望 #E1是log(x)的条件期望值--E(log(x)|x>y) E1 = function(a,b, y){   -log(b)+1/upgamma(a,y*b)*( digamma(a)*gamma(a)*(1-exp(-y*b)*s.1(y,4,a,b))                              -gamma(a)*exp(-y*b)*log(y*b)*s.1(y,4,a,b)                              +gamma(a)*exp(-y*b)*s.2(y,4,a,b))/(pgamma(y*b,a,lower.tail=FALSE)*gamma(a))

负红颜

不懂数学，哈哈哈哈哈