CosmoMC: 3D Likelihood Grids?
-
- Posts: 15
- Joined: September 14 2005
- Affiliation: Univ. at Buffalo, SUNY
- Contact:
CosmoMC: 3D Likelihood Grids?
Hi Antony --
I am interested in producing gridded likelihoods for three parameters, i.e. the equivalent of the <file_root>_2D_x_y_likes files produced by getdist, except for three parameters instead of two.
Any easy way to do this?
-- Will
I am interested in producing gridded likelihoods for three parameters, i.e. the equivalent of the <file_root>_2D_x_y_likes files produced by getdist, except for three parameters instead of two.
Any easy way to do this?
-- Will
-
- Posts: 15
- Joined: September 14 2005
- Affiliation: Univ. at Buffalo, SUNY
- Contact:
Re: 3D Likelihood Grids?
This should have had "COSMOMC" in the subject line!Will Kinney wrote:Hi Antony --
I am interested in producing gridded likelihoods for three parameters, i.e. the equivalent of the <file_root>_2D_x_y_likes files produced by getdist, except for three parameters instead of two.
Any easy way to do this?
Looks relatively straightforward to hack a new subroutine for getdist to accomplish this, but there must be a more elegant solution...
-- Will
-
- Posts: 1944
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: 3D Likelihood Grids?
It is not in GetDist, but I think you can do it using CosmoloGui or the matlab scripts that that uses (e.g. see the 3D contour plots). Sarah may be able to help further here.
-
- Posts: 144
- Joined: September 24 2004
- Affiliation: University College London (UCL)
- Contact:
CosmoMC: 3D Likelihood Grids?
Hi, yes, the relevant routines are in the cosmomc/mscripts directory (cosmologui)
e.g. could call it like this:
[xim,xn]=samp3imw([x' y' z']',samples,resolution); % grid
handle=plot3dconf(xn(1,:),xn(2,:),xn(3,:),xim,0.68,5,1); % isosurface
e.g. could call it like this:
[xim,xn]=samp3imw([x' y' z']',samples,resolution); % grid
handle=plot3dconf(xn(1,:),xn(2,:),xn(3,:),xim,0.68,5,1); % isosurface
-
- Posts: 15
- Joined: September 14 2005
- Affiliation: Univ. at Buffalo, SUNY
- Contact:
Re: CosmoMC: 3D Likelihood Grids?
Thanks much!Sarah Bridle wrote:Hi, yes, the relevant routines are in the cosmomc/mscripts directory (cosmologui)
e.g. could call it like this:
[xim,xn]=samp3imw([x' y' z']',samples,resolution); % grid
handle=plot3dconf(xn(1,:),xn(2,:),xn(3,:),xim,0.68,5,1); % isosurface
-- Will