Usage
schechter.fit(data, vmax = NA, knee, slope, norm, knee.alt = NA, slope.alt = NA, norm.alt = NA, kneelo = -Inf, slopelo = -Inf, normlo = 0, kneehi = Inf, slopehi = Inf, normhi = Inf, fixk1 = FALSE, fixs1 = FALSE, fixn1 = FALSE, fixk2 = FALSE, fixs2 = FALSE, fixn2 = FALSE, range = range(data), lim1 = NA, lim2 = NA, numlim = 1, method = "nlminb", volume = max(vmax), bw = 0.1, mag = FALSE, log = FALSE, null = 1E-9, error = "jack", subvol = 10, sampnum = subvol, msun = solar("r"))
Arguments
vmax
vector of maximum comoving volumes within which object could lie
knee
the knees(s) of the Schechter function (L_star/M_star)
slope
the slope(s) of the Schechter function (alpha)
norm
the normalisation(s) of the Schechter function (phi_star)
knee.alt
alternative knees(s) of the Schechter function (L_star/M_star)
slope.alt
alternative slope(s) of the Schechter function (alpha)
norm.alt
alternative normalisation(s) of the Schechter function (phi_star)
kneelo
a lower bound on the knee parameter
slopelo
a lower bound on the slope parameter
normlo
a lower bound on the norm parameter
kneehi
an upper bound on the knee parameter
slopehi
an upper bound on the slope parameter
normhi
an upper bound on the norm parameter
fixs1
fix the first slope?
fixn1
fix the first normalisation?
fixk2
fix the second knee?
fixs2
fix the second slope?
fixn2
fix the second normalisation?
range
data range of interest
lim1
lower data limit for fitting
lim2
upper data limit for fitting
numlim
lower number (per bin) limit for fitting
method
choice of 'nlminb' (recommended) or one of 'optim's minimisation methods (e.g., 'Nelder-Mead')
volume
total volume across which the data has been collected (default 1 if vmax = NA)
mag
are the input data magnitudes?
log
are the input data logged?
error
parameter error estimation method [jack/boot]
subvol
number of sub-volumes to split the input data into
sampnum
number of samplings to be made for bootstrapping method
msun
absolute magnitude of the sun (to be used in calculation of the luminosity density should mag=TRUE)