Plot range for line of fit using date on the x-axis in gnuplot -

-

i trying plot fit function (col 3 versus col 7) dataset given below 

 28-08-1991 20-12-1992  24-04-1992  -263347200  -221875200  -242611200  0.060859   20-12-1992 02-09-1993  27-04-1993  -221875200  -199756800  -210816000  0.064681   10-09-1996 13-09-1997  13-03-1997  -104371200  -72576000   -88473600   0.095728 ......

i script plot fit range of dates, between 1995 2003. unable range function right since gnuplot handle dates past year 2000.

please find code below. 

 set xdata time  set timefmt "%d-%m-%y"  set xrange ["01-01-1991":"01-01-2015"]   set yrange [0.0:0.4]  offset = 24*60*60   f(x) = 1.e-7*m*(x-offset)+c   fit f(x) "file"  using 3:7  via m,c  set fit logfile 'file.log'   plot '(x < "24-02-1993" ) ? f(x) : 1/0' title "best fit" lines ls 5

the variable x in plot contain time in form of seconds since 1/1/1970. in order compare date, need convert date same format strptime.

if additionally remove ' around function plot in mwe, get:

plot (x<strptime("%d-%m-%y","24-02-1993"))?f(x):1/0 title "best fit" lines ls 5

note you've fitted f(x) on whole dataset, not on dates before 24-02-1993, it's plot restricted these dates.


Comments

Post a Comment

Popular posts from this blog

sequelize.js - Sequelize group by with association includes id -

-
so when requesting group sequelize follows: return models.workingcalendar .findall({ attributes: [ 'workingcalendar.periodid', 'workingcalendar.date', 'period.name' ], include: [ { model: models.period, attributes: [] } ], where: { getsudoid: currentgetsudo.id, unitplantid: unitplantid }, group: ['workingcalendar.periodid', 'workingcalendar.date', 'period.name'], }); sequelize run query: select [workingcalendar].[id], [workingcalendar].[periodid], [workingcalendar].[date], [period].[name] [workingcalendars] [workingcalendar] left outer join [periods] [period] on [wor...
Read more

c++ - Migration from QScriptEngine to QJSEngine -

-
i'm migrating qscriptengine code on qjsengine. now, have: class pars { public: static qscriptvalue printmainlog(qscriptcontext* c, qscriptengine* e); }; qscriptvalue pars::printmainlog(qscriptcontext* c, qscriptengine* e) { //some actions return e->globalobject().property(""); } ... qscriptengine engine; ... engine.globalobject().setproperty("printlog",engine.newfunction(pars::printmainlog)); so, user can put printlog("what ever"); in application in , example, qlineedit , function pars::printmainlog evalute. is there way qjsengine? so, user put same printlog("what ever");? way find here , user should put logger.printlog("what ever"); logger class inherited qobject printlog slot.
Read more

delphi - Take screenshot in webcam using VFrames in Console Application -

-
i'm doing program take picture of webcam using delphi xe2 , vframes achieve this, problem have figured out, in graphic application works fine, when use unit in console application, returns me error saying first chance exception @ $76b6b727. exception class eaccessviolation message 'access violation @ address 004a271b in module 'console.exe'. read of address 00000260'. process console.exe (3676) my unit : unit webcam; interface uses sysutils, windows, vcl.imaging.jpeg, vcl.graphics, vsample, vframes, classes; type twebcam = class private procedure newvideoframeevent(sender: tobject; width, height: integer; dataptr: pointer); public constructor create; destructor destroy; override; procedure capture_webcam(take_name: string); end; var web_image: tvideoimage; name_screen: string; implementation constructor twebcam.create; begin inherited create; end; destructor twebcam.destroy; begin inherited destroy; end; ...
Read more