Displaying 2 results from an estimated 2 matches for "a_guess".
Did you mean:
q_guess
2009 Jul 30
0
Constrained MLE of fixed mean Singh-Maddala distribution
...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_...
2009 Jul 31
0
MLE estimation of constrained mean Singh-Maddala distribution
...mad<-function(p, y, x)
{
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_ha...