Could anyone recommend a package for visualizing a likelihood function of two scalar parameters? I would prefer a three-dimensional plot similar to the kind Mathematica is known for, perhaps generated by a package not specific to likelihood. David ______________________________ David R. Bickel Ottawa Institute of Systems Biology BMI Dept., University of Ottawa http://www.oisb.ca/members.htm
See ?persp Uwe Ligges David Bickel wrote:> Could anyone recommend a package for visualizing a likelihood function > of two scalar parameters? I would prefer a three-dimensional plot > similar to the kind Mathematica is known for, perhaps generated by a > package not specific to likelihood. > > David > > ______________________________ > David R. Bickel > Ottawa Institute of Systems Biology > BMI Dept., University of Ottawa > > http://www.oisb.ca/members.htm > > ______________________________________________ > 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.
On 12/10/2007 11:26 AM, David Bickel wrote:> Could anyone recommend a package for visualizing a likelihood function > of two scalar parameters? I would prefer a three-dimensional plot > similar to the kind Mathematica is known for, perhaps generated by a > package not specific to likelihood.The basic strategy is to calculate a matrix of values of the likelihood, corresponding to all combinations of values of two vectors. Then persp() (in graphics) or persp3d() (in rgl) can display the surface. I don't know of anything that corresponds to curve() to do the matrix calculations automatically, but it may well exist. Duncan Murdoch> > David > > ______________________________ > David R. Bickel > Ottawa Institute of Systems Biology > BMI Dept., University of Ottawa > > http://www.oisb.ca/members.htm > > ______________________________________________ > 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.
On Mon, 10 Dec 2007, Uwe Ligges wrote:> See ?perspIf you have a suitable OS, I would recommend using package rgl and persp3d: the ability to change viewpoint interactively is very useful.> > Uwe Ligges > > David Bickel wrote: >> Could anyone recommend a package for visualizing a likelihood function >> of two scalar parameters? I would prefer a three-dimensional plot >> similar to the kind Mathematica is known for, perhaps generated by a >> package not specific to likelihood. >> >> David >> >> ______________________________ >> David R. Bickel >> Ottawa Institute of Systems Biology >> BMI Dept., University of Ottawa >> >> http://www.oisb.ca/members.htm >> >> ______________________________________________ >> 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. >-- Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595