Re: [R] Estimates at each iteration

From: Jonathan P Daily <jdaily_at_usgs.gov>
Date: Fri, 08 Apr 2011 11:57:37 -0400

I'm assuming that English is not your first language. In light of that, I would suggest that you post your question in a language that we both speak: R. Please include:

  1. E-step and M-step functions.
  2. The proposed optim/maxLik call.

Mostly because (at least optim) does not have a parameter called "limiter".

HtoH,
Jon



Jonathan P. Daily
Technician - USGS Leetown Science Center 11649 Leetown Road
Kearneysville WV, 25430
(304) 724-4480
"Is the room still a room when its empty? Does the room,  the thing itself have purpose? Or do we, what's the word... imbue it."

r-help-bounces_at_r-project.org wrote on 04/08/2011 08:16:36 AM:

> [image removed]
>
> [R] Estimates at each iteration
>
> hello bye
>
> to:
>
> r-help
>
> 04/08/2011 11:45 AM
>
> Sent by:
>
> r-help-bounces_at_r-project.org
>
> Dear Sir\ Madam
> I am trying to maximise a complicated loglikelihood function via
> EM algorithm
> which consists of two step "E-step and M-step"and in this case I need to

> maximize the expected of like lihood function " which I get from E-
step" and
> take those estimates of parameter to update the E-step and repate these
till
> convargence has been met. So I need to know if I control the number
> of iterlim
> in maxLik or optim Package by puting this number equals to one to
maximize
> M-step and getting just the estimates of parameter for one
> iteration and hence
> use them to update the E-step and again maximaize M-step by using maxLik
or
> optim by contoling the number of limiter and putting the number equals
to 1
> is it true and acceptance? I mean if I do this my results can be
> reasonable?. please I am waiting your response.
> Regards
> Fauzia Taweab
> phd.student.
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Fri 08 Apr 2011 - 15:59:37 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 08 Apr 2011 - 16:00:28 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive