Dear Berend
Thanks for your response. I will write my code and estimate the parameters using functions optim() or nlminb(). For the code that I am trying to build for estimation the parameters using mle via EM ahgorithm, I will rewrite it and reformulated my question about, but before that I need to inquire whether  i can control the number of iteration ( maxit) in optim function to get the estimation at each step (maxit =1), (M step) and then use them to compute the expected values of the log likelihood function ( Estep), estimtng using MLE via EM,.  Thank you very much.
	[[alternative HTML version deleted]]


