[R] problem with uniroot
Ravi Varadhan
RVaradhan at jhmi.edu
Wed Jun 3 17:42:30 CEST 2009
Can you send a minimal, reproducible code so that we can try to help you?
Ravi.
----------------------------------------------------------------------------
-------
Ravi Varadhan, Ph.D.
Assistant Professor, The Center on Aging and Health
Division of Geriatric Medicine and Gerontology
Johns Hopkins University
Ph: (410) 502-2619
Fax: (410) 614-9625
Email: rvaradhan at jhmi.edu
Webpage:
http://www.jhsph.edu/agingandhealth/People/Faculty_personal_pages/Varadhan.h
tml
----------------------------------------------------------------------------
--------
-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On
Behalf Of daniele riggi
Sent: Wednesday, June 03, 2009 11:34 AM
To: r-help at r-project.org
Subject: [R] problem with uniroot
Hi R-users,I'm trying to solve a non linear equation, to find the degrees of
freedom of a mixture of t student. I'm sure i wrote the minimization
equation in the right way, but when i try to run the EM algorithm to
estimate the parameters of the mixture, the following error will appear:
Error in uniroot(function(z) log(z/2) - digamma(z/2) + 1 - log((z + d)/2) +
:
f() values at end points not of opposite sign
I don't know how solve this problem, i have already tried using the function
optimize, but it doesn't work.
Someone could help me?
thanks a lot
--
Dr. Daniele Riggi, PhD student
University of Milano-Bicocca
Department of Statistics
Building U7, Via Bicocca degli Arcimboldi, 8
20126 Milano, Italy
cell. +39 328 3380690
mailto: daniele.riggi at gmail.com
[[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.
More information about the R-help
mailing list