analyseMCMC_1dpdfs.f90 File Reference

(AnalyseMCMC version 0.7.1, revision 448, hash 2f54722 (2013-11-19))

Routines to compute and plot one-dimensional PDFs. More...

Functions/Subroutines
subroutine	pdfs1d (exitcode)
	Plot 1D marginalised PDFs. More...
subroutine	bin_data_1d_a (ndat, xdat, ydat, norm, nbin, xmin, xmax, xbin, ybin)
	"Bin data" 1D by measuring the amount of likelihood in each bin More...
subroutine	verthist (n, x1, y1, style)
	Plot a 1D vertical histogram. x is the left of the bin! More...
subroutine	horzhist (n, x, y)
	Plot a 1D horizontal histogram. More...
subroutine	smoothpdf1d (ybin, nbin, smooth)
	Smooth a 1D PDF. More...

Detailed Description

Routines to compute and plot one-dimensional PDFs.

Function/Subroutine Documentation

subroutine pdfs1d

(

integer, intent(out)

exitcode

)

Plot 1D marginalised PDFs.

Return values

exitcode

Exit code: 0=ok

References pginitl(), report_undefined_parameter(), rev180(), rev24(), rev360(), smoothpdf1d(), and verthist().

Referenced by analysemcmc().

subroutine bin_data_1d_a	(	integer, intent(in)	ndat,
		real, dimension(ndat), intent(in)	xdat,
		real, dimension(ndat), intent(in)	ydat,
		integer, intent(in)	norm,
		integer, intent(in)	nbin,
		real, intent(inout)	xmin,
		real, intent(inout)	xmax,
		real, dimension(nbin+1), intent(out)	xbin,
		real, dimension(nbin+1), intent(out)	ybin
	)

"Bin data" 1D by measuring the amount of likelihood in each bin

Parameters

ndat	Number of data points
xdat	Data to be binned (n points)
ydat	Likelihoods of the n points
norm	Normalise histogram (1) or not (0)
nbin	Desired number of bins
xmin	Minimum value of the binning range. Set xmin=xmax to auto-determine (I/O)
xmax	Maximum value of the binning range. Set xmin=xmax to auto-determine (I/O)

Return values

xbin	Binned data, location of the bins. The x values are the left side of the bin!
ybin	Binned data, height of the bins.

subroutine verthist	(	integer, intent(in)	n,
		real, dimension(n+1), intent(in)	x1,
		real, dimension(n+1), intent(in)	y1,
		integer, intent(in)	style
	)

Plot a 1D vertical histogram. x is the left of the bin!

Parameters

n	Number of bins
x1	Location of the bins
y1	Height of the bins
style	Histogram style: 0: don't drop vertical lines to 0 at each bin, 1: do this, 2: fill the histogram

Referenced by pdfs1d().

subroutine horzhist	(	integer, intent(in)	n,
		real, dimension(n), intent(in)	x,
		real, dimension(n), intent(in)	y
	)

Plot a 1D horizontal histogram.

Parameters

n	Number of bins
x	Locations of the bins
y	Heights of the bins

subroutine smoothpdf1d	(	real, dimension(nbin), intent(inout)	ybin,
		integer, intent(in)	nbin,
		integer, intent(in)	smooth
	)

Smooth a 1D PDF.

Parameters

ybin	Height of the bins
nbin	Number of bins
smooth	Smoothing factor

References savgol().

Referenced by animation(), and pdfs1d().