Hi Steve, Actually I figured that out after I emailed you but I forgot to email you before I left office. Anyway thanks a lot.And thanks for the note. I'll do that from now on.
Subodh On Fri, Jun 4, 2010 at 4:54 PM, Steve Lianoglou < mailinglist.honey...@gmail.com> wrote: > Hi Subodh, > > Minor note: please keep replies on list so everyone benefits from > answers/questions. > > Now: > > > On Fri, Jun 4, 2010 at 2:49 PM, Subodh Acharya <shoeb...@gmail.com> wrote: > > Thanks a lot Steve, > > It worked. I appreciate. But I have another question, may be thats > trivial. > > Say, it doesn't converge at itermax. I need to display error message > saying > > "values don't converge at at itermax". Is there another statement that > > accompanies while for this? > > Perhaps you can use "warning(...)". You can check whether (or not) > your Fpt - Fpi is > tolerance after you loop. If it is, you know your > loop terminated because you hit itermax, and not because your algo > converged: > ... > ... > iter <- 0 > while (((Fpt - Fpi) > tolerance) && (iter < itermax)) { > Fpi = Fpt > Fpt = K*Time + M*S*log(1+ Fpi/(M*S)) > Fp0 = Fpt > iter <- iter + 1 > } > > if (Fpt - Fpi > tolerance) { > warning("Algorithm reached itermax and did not converge") > } > > ... > > -- > Steve Lianoglou > Graduate Student: Computational Systems Biology > | Memorial Sloan-Kettering Cancer Center > | Weill Medical College of Cornell University > Contact Info: http://cbio.mskcc.org/~lianos/contact > -- Acharya, Subodh [[alternative HTML version deleted]] ______________________________________________ R-help@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.