Login or Register to Ask a Question and Join Our Community

Sponsored Content
Full Discussion: Grep final set of parameters from fit.log gnuplot file
Top Forums Shell Programming and Scripting Grep final set of parameters from fit.log gnuplot file Post 302783475 by kayak on Wednesday 20th of March 2013 01:40:38 PM
Old 03-20-2013
Grep final set of parameters from fit.log gnuplot file
I would like to grep the final set of fit parameters from a gnuplot log file to form columns that look like this.
Code:
a_1001 b_1001 x_1001
a_1002 b_1002 x_1002
a_1003 b_1003 x_1003
   .          .         .   
   .          .         .
   .          .         .
a_1250  b_1250 c_1250

At present, the 'fit.log' file looks like this (also see attachment).

Code:

*****************************************************************




 Iteration 0
 WSSR        : 13.5219           delta(WSSR)/WSSR   : 0
 delta(WSSR) : 0                 limit for stopping : 1e-05
 lambda	  : 1.1547

initial set of free parameter values

a_1001          = 1
b_1001          = 1
x_1001          = 1

After 8 iterations the fit converged.
final sum of squares of residuals : 0.00106152
rel. change during last iteration : -1.34952e-06

degrees of freedom    (FIT_NDF)                        : 15
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.00841238
variance of residuals (reduced chisquare) = WSSR/ndf   : 7.07681e-05

Final set of parameters            Asymptotic Standard Error
=======================            ==========================

a_1001          = 0.497897         +/- 0.001983     (0.3982%)
b_1001          = 0.00459204       +/- 0.005632     (122.6%)
x_1001          = 1.57077          +/- 0.6132       (39.04%)


correlation matrix of the fit parameters:

               a_1001 b_1001 x_1001 
a_1001          1.000 
b_1001         -0.003  1.000 
x_1001         -0.000  0.002  1.000 


****************************************************************

 Iteration 0
 WSSR        : 13.5246           delta(WSSR)/WSSR   : 0
 delta(WSSR) : 0                 limit for stopping : 1e-05
 lambda	  : 1.1547

initial set of free parameter values

a_1002          = 1
b_1002          = 1
x_1002          = 1

After 8 iterations the fit converged.
final sum of squares of residuals : 0.0019177
rel. change during last iteration : -1.59332e-09

degrees of freedom    (FIT_NDF)                        : 15
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.0113069
variance of residuals (reduced chisquare) = WSSR/ndf   : 0.000127847

Final set of parameters            Asymptotic Standard Error
=======================            ==========================

a_1002          = 0.49692          +/- 0.002665     (0.5363%)
b_1002          = 0.00889775       +/- 0.007585     (85.24%)
x_1002          = 1.57078          +/- 0.4262       (27.13%)


correlation matrix of the fit parameters:

               a_1002 b_1002 x_1002 
a_1002          1.000 
b_1002         -0.006  1.000 
x_1002         -0.000  0.002  1.000

I only need the final set of parameters as a table and I am wondering if anyone could help using awk or grep to get these parameters out.
fit.txt (2.8 KB) 
 

9 More Discussions You Might Find Interesting

1. UNIX for Dummies Questions & Answers

least-square fit in Gnuplot

Does anyone know how to find the best least square fit in Gnuplot? (6 Replies)
Discussion started by: cosmologist
6 Replies

2. Shell Programming and Scripting

Call single function multiple times diff set of parameters

Okay, not sure if it can be done, I would think it could be done and I'm just having a hard time with it. fun_close_wait(){ ipVar="${1} ${2}" portVar=`cat "${5}" | cut -d' ' -f 1` for ip in $ipVar do for port in $portVar do netstatVar=`netstat -n | grep... (4 Replies)
Discussion started by: cbo0485
4 Replies

3. Solaris

Which file is read by kernel to set its default system kernel parameters values?

Hi gurus Could anybody tell me which file is read by kernel to set its default system kernal parameters values in solaris. Here I am not taking about /etc/system file which is used to load kernal modules or to change any default system kernal parameter value Is it /dev/kmem file or something... (1 Reply)
Discussion started by: girish.batra
1 Replies

4. Shell Programming and Scripting

gnuplot set key issue

Hi I have a plotting function like this- plot '../processing_old/time0.txt' using 1:6 title "1.6" with boxes fs solid 0.05 ls 5,\ 'time0.txt' using 1:6 title "1.7" with boxes fs solid 0.05 ls 6,\ '../processing_new/time0.txt' using 1:6 title "2.0" with boxes fs solid... (8 Replies)
Discussion started by: jamie_123
8 Replies

5. Shell Programming and Scripting

Log-log plot in gnuplot and limiting the xrange

I am trying to fit a plot in gnuplot using logscale. I have 50000 data points. At first I fit plot in this way. f(x) = b + m*x fit f(x) "xyMSD-all-mal-cel-iso-bcm-thermo2.dat" using 1:2 via m,b I got slope value. Then I tried to get slope value at different range as below. fit f(x)... (2 Replies)
Discussion started by: vjramana
2 Replies

6. Linux

Gnuplot terminal set to 'unknown'

I installed version 4.6.4 of gnuplot recently on my system running 10.04 ubuntu linux. I am getting error Terminal type set to 'unknown' Also when i try to set term X11 or wxt i get error Terminal type set to 'unknown' ^ unknown or ambiguous terminal... (5 Replies)
Discussion started by: Ujjwal1982
5 Replies

7. Shell Programming and Scripting

How do i check if all parameters are set in bash?

Hi, I have 4 parameters passed to my shell and i validate if all four are entered using the following snippet: if then echo "Not entered" else echo "entered" fi I get the following output as 'Not entered even when i enter the values for all prompts. Please advise. Thanks. (5 Replies)
Discussion started by: Jesshelle David
5 Replies

8. Shell Programming and Scripting

Getopt eval set parameters not happening when the script is called through an application!

Hi, There is an Informatica tool through which unix scripts can be called. Now the requirement in my project is that a parent script calls a child script but this parent script has to be called through the Informatica tool. In Parent script I'm using TEMP=`getopt -o x:y: -l area:,volume:... (1 Reply)
Discussion started by: Panna
1 Replies

9. Shell Programming and Scripting

Different set key multiplot gnuplot

Hello everybody, I am using Windows 10. I need to plot several graphs simultaneously with Gnuplot. The script is: GNUTERM = "wxt" set terminal wxt set termoption enhanced set encoding iso_8859_1 set multiplot layout 2,3 set xtics out set ytics out set xtics nomirror set ytics... (4 Replies)
Discussion started by: supernono06
4 Replies
All times are GMT -4. The time now is 05:59 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy