-
Notifications
You must be signed in to change notification settings - Fork 1
Home
John Robinson edited this page Nov 1, 2015
·
8 revisions
This is script-driven global non-linear least squares parameter estimation written in Fortran 90. The software was originally written to analyze time resolved fluorescence data (decay). The observed fluorescence decay is the intrinsic decay convolved with the instrument resonse function (IRF).
- One must have global.in file (a .conf or .ini file) in the working directory of the .ans file. This controls the nature of the analysis, like the target change in chi^2 for stopping the analysis. For this reason, the .in file must be local.
- When data needs to be pre- or post-processed in some way during the anslysis, this is specified in the data file itself. Specification for processing the data is in the header of the data file.
- This is compiled softare. Additional models must be compiled.
- This was written using the Digital/Compac/Fujitsu/Intel F90 compiler, which uses digital extensions of the F90 standared. This restricts one to development on DOS.
- Srcripts (usually called .ans files) are written in a custom format. These days, one would write in YAML or similar.
- Software is very unforgiving of errors in data format or script files.
- install globalcurve in
c:\bin\ - install pgplot plotting library (dependency) in
c:\
Set system Path (environement variable) to recognize the globalcurve executable. Value on my computer is
%SystemRoot%\system32;%SystemRoot%;%SystemRoot%\System32\Wbem;%SYSTEMROOT%\System32\WindowsPowerShell\v1.0\;C:\bin\GlobalCurve
install at c:
In windows, use PowerShell. To access documents, which are shared using VMWARE, cd z:. cd to directory of interest.
> globalcurve
3.4100000E-02 50 700 80
0.0 1.0 1.0
0.0 2.0 1.4
1.0 3.0 1.7
0.0 3.0 1.7
0.0 1.0 1.0
0.0 3.0 1.7
0.0 3.0 1.7
- data format is that of a .prn file (can export from excel)
- windows formated EOL: CR+LF
- usually given a .glo extension.
- Format: Columns 1 and 2 are IRF and #counts in the donor channel, respectively. Rows are channels 0, 1, 2, .... Column 3, which specifies the precision of the data (sd) is optional. When omitted, the analysis assumed Poisson noise, where sd = sqrt(n). The header contains, successively,
- time/channel (ns/chan)
- start channel of the convolution
- stop channel of chi^2 calculation
- start channel of the chi^2 calculation
- in default.in set simulation to true.
- run globalcurve as usual
- enter the .ans file name that contains the model you wish to simulate
- give the tcal, lims1, lims2, lims3, which appear in the header of .glo files
Enter answer filename: exp1.ans
Enter simulation pearameters for expt # 1
Enter tcal, lims1, lims2, lims3:
0.016,55,350,70
Enter fwhm of the lamp function (in units of tcal):
0.1
Simulation mode:
Add noise (t,f)? t