YujiSODE / crvFit

Curve fitting tool based on the least absolute value method and the Monte Carlo method

crvFit

Curve fitting tool based on the least absolute value method and the Monte Carlo method.
GitHub: https://github.com/YujiSODE/crvFit

Copyright (c) 2020 Yuji SODE <[e-mail]>
This software is released under the MIT License.
See LICENSE or http://opensource.org/licenses/mit-license.php

---

1. Description

Namespace ::crvFit has procedures to estimate curve fitting parameters.

• ::crvFit::setFunction ?formula?;
  Procedure that adds a mathematical function that returns a result of f(x) defined at namespace ::crvFit.
  Generated function is called crvFit_F(x), which is expressed with variable $x, parameters (in upper case only), numerical values and mathematical functions in Tcl expressions.

  • $formula: a formula for f(x), which is expressed with variable $x, parameters (in upper case only), numerical values and mathematical functions in Tcl expressions, default value is {$A*$x+$B}

• ::crvFit::setRange parameterName range;
  Procedure that sets range of random parameters.

  • $parameterName: a name of parameter to set
  • $range: a range for a random number that is expressed as "min,max"

• ::crvFit::setParameter parameterName value;
  Procedure that sets parameter for generated function based on given values.

  • $parameterName: a name of parameter to set
  • $value: a numerical value

• ::crvFit::loadXY ?xyList?;
  Procedure that loads list of x-y data.
  Stored data is returned when $xyList is not specified.
  Format of x-y data is "x,y".

  • $xyList: a list of x-y data, and every element is expressed as "x,y"

• ::crvFit::estimate ?n?;
  Procedure that estimates parameters.
  Returned value is named list.

  • $n: sample size used for a single estimation based on the least absolute value method, default size is 100

• ::crvFit::estimateMC N ?n?;
  Procedure returns result of the least absolute value method with the Monte Carlo approximation.
  Returned value is named list.

  • $N: a number of estimated parameter sets in order to estimate average
  • $n: sample size used for a single estimation based on the least absolute value method, default size is 100

• ::crvFit::outputLog namedList fileName;
  Procedure that outputs estimation log and estimated tcl math function.
  Output function is called crvFitLog_F(x).

  • $namedList: a named list that is returned by ::crvFit::estimate ?n?; or ::crvFit::estimateMC N ?n?;
  • $fileName: a name of file to output

2. Script

Tcl

It requires Tcl 8.6+.

Main script

• crvFit.tcl

Additional script

• varCSV.tcl: CSV file dealing interface

Sample script and sample data

• sample01.tcl,
• sample02.tcl
• sample03.tcl
• sample04.tcl
• sample05.tcl
• csvSample.csv

3. Library list

• Sode, Y. 2018. lSum_min.tcl: https://gist.github.com/YujiSODE/1f9a4e2729212691972b196a76ba9bd0
• Sode, Y. 2019. nonparametricP.tcl: the MIT License; https://gist.github.com/YujiSODE/4d5d004a6e98f2d17b42940fcc97c7d6