FORCES
FORtran lib for Comp. Env. Sys.
|
Fits a straight line to input data by minimizing chi^2. More...
Public Member Functions | |
real(sp) function, dimension(:), allocatable | linfit_sp (x, y, a, b, siga, sigb, chi2, model2) |
real(dp) function, dimension(:), allocatable | linfit_dp (x, y, a, b, siga, sigb, chi2, model2) |
Fits a straight line to input data by minimizing chi^2.
Given a set of data points x(1:ndata), y(1:ndata), fit them to a straight line \( y = a+bx \) by minimizing chi2.
Model I minimizes y vs. x while Model II takes the geometric mean of y vs. x and x vs. y. Returned is the fitted line at x. Optional returns are a, b and their respective probable uncertainties siga and sigb, and the chi-square chi2.
[in] | real(sp/dp) :: x(:) | 1D-array with input x |
[in] | real(sp/dp) :: y(:) | 1D-array with input y |
[in] | logical, optional :: model2 | If present, use geometric mean regression instead of ordinary least square |
[out] | real(sp/dp), dimension(M) :: a | intercept |
[out] | real(sp/dp), dimension(M) :: b | slope |
[out] | real(sp/dp), dimension(M) :: siga | error on intercept |
[out] | real(sp/dp), dimension(M) :: sigb | error on slope |
[out] | real(sp/dp) :: chisq | Minimum chi^2 |
Definition at line 84 of file mo_linfit.f90.
real(dp) function, dimension(:), allocatable mo_linfit::linfit::linfit_dp | ( | real(dp), dimension(:), intent(in) | x, |
real(dp), dimension(:), intent(in) | y, | ||
real(dp), intent(out), optional | a, | ||
real(dp), intent(out), optional | b, | ||
real(dp), intent(out), optional | siga, | ||
real(dp), intent(out), optional | sigb, | ||
real(dp), intent(out), optional | chi2, | ||
logical, intent(in), optional | model2 | ||
) |
Definition at line 98 of file mo_linfit.f90.
real(sp) function, dimension(:), allocatable mo_linfit::linfit::linfit_sp | ( | real(sp), dimension(:), intent(in) | x, |
real(sp), dimension(:), intent(in) | y, | ||
real(sp), intent(out), optional | a, | ||
real(sp), intent(out), optional | b, | ||
real(sp), intent(out), optional | siga, | ||
real(sp), intent(out), optional | sigb, | ||
real(sp), intent(out), optional | chi2, | ||
logical, intent(in), optional | model2 | ||
) |
Definition at line 189 of file mo_linfit.f90.