Performs a correlated fit of a constant to data generated with
bootstrap.effectivemass.
fit.effectivemass(cf, t1, t2, useCov = FALSE, replace.na = TRUE,
boot.fit = TRUE, autoproceed = FALSE, every)An object of class effectivemass generated by a call to
bootstrap.effectivemass.
The fit range. If several correlators are fitted, this is
automatically replicated accordingly. The fit range is adjusted such that
NAs are removed from the fit. They must fulfill \(t_1<t_2\).
For symmetric correlators, they must both run from 0 to Time/2-1,
otherwise from 0 to Time-1.
Use the correlated chisquare. This works only for not too noisy data.
The functions inverted to determine the effective mass
values might, due to fluctuations, return NA. If
replace.na=TRUE, these are reaplaced in the bootstrap samples by
randomly chosen values from the distribution that are not NA.
Otherwise the fits in which the NA values occur will fail.
If set to FALSE, the effective mass fit is not
bootstrapped, even though bootstrap samples are still used to estimate the
variance-covariance matrix for the correlated fit. This is a useful
time-saver if error information is not strictly necessary. Of course, this
affects the return values related to the bootstrap, which are set to
NA.
When the inversion of the variance-covariance matrix
fails, the default behaviour is to abort the fit. Setting this to
TRUE means that the fit is instead continued with a diagonal inverse
of the variance-covariance matrix.
Fit only a part of the data points. Indices that are not
multiples of every are skipped. If no value is provided, all points
are taken into account.
An object with class effectivemassfit is returned. It
contains all the data of the input object effMass with the following
additional member objects:
opt.res: the object returned by the optim on the original
data.
massfit.tsboot: the bootstrap values of the mass and the chisquare
function.
ii: the index array of data used in the fit.
invCovMatrix: the inverse covariance matrix.
dof: the degrees of freedom of the fit.
t1,t2: the fit range.
A correlated chisquare minimisation is performed on the original data as
well as on all bootstrap samples generated by
bootstrap.effectivemass. The inverse covariance matrix is generated
as described in hep-lat/9412087 in case of too little data to relibably
estimate it.
C.Michael, A.McKerrell, Phys.Rev. D51 (1995) 3745-3750, hep-lat/9412087
bootstrap.effectivemass,
bootstrap.gevp, gevp2cf,
invertCovMatrix
# NOT RUN {
data(samplecf)
samplecf <- bootstrap.cf(cf=samplecf, boot.R=99, boot.l=2, seed=1442556)
effmass <- fit.effectivemass(bootstrap.effectivemass(cf=samplecf), t1=15, t2=23)
summary(effmass)
plot(effmass, ylim=c(0.14,0.15))
# }
Run the code above in your browser using DataCamp Workspace