search for: q_guess

Displaying 2 results from an estimated 2 matches for "q_guess".

Did you mean: a_guess
2009 Jul 30
0
Constrained MLE of fixed mean Singh-Maddala distribution
...q_iter<-exp(p[3]) #NEED TO PUT IN CONTITION ENSURING THAT 'q_iter-1/a_iter>0'. z<-(-sum(log(dsinmad(y, a=a_iter, scale=b_iter, q.arg=q_iter, log=FALSE)))) } #Set values for initial guess of the parameters. a_guess<-2 b_guess<-3 q_guess<-7 q_guess-1/a_guess p0<-c(log(a_guess), log(b_guess), log(q_guess)) #Optimize to find minimum of the negative likelihood function est_p_A<-optim(par=p0, fn=neg.LL_sinmad, y=y) est_p_A$par a_hat1<-exp(est_p_A$par[1]) b_hat1<-exp(est_p_A$par[2]) q_hat1<-exp(est_...
2009 Jul 31
0
MLE estimation of constrained mean Singh-Maddala distribution
...a_iter<-exp(p[1]) b_iter<-exp(p[2]) q_iter<-exp(p[3]) #NEED TO PUT IN CONTITION ENSURING THAT 'q_iter-1/a_iter>0'. z<-(-sum(log(dsinmad(y, a=a_iter, scale=b_iter, q.arg=q_iter, log=FALSE)))) } #Set values for initial guess of the parameters. a_guess<-2 b_guess<-3 q_guess<-7 q_guess-1/a_guess p0<-c(log(a_guess), log(b_guess), log(q_guess)) #Optimize to find minimum of the negative likelihood function est_p_A<-optim(par=p0, fn=neg.LL_sinmad, y=y) est_p_A$par a_hat1<-exp(est_p_A$par[1]) b_hat1<-exp(est_p_A$par[2]) q_hat1<-exp(est_p_A$par[3])...