gnuplot-info — For questions and discussion about gnuplot

Re: [Gnuplot-info] Regarding plotting stacked histograms with custom colors

[Debostuti G. <deb...@gm...>]

Dear Ethan

I installed 5.4.6 but it is showing terminal type as unknown. pngcairo etc
are no longer available to me. Looks like the reason is the absence of
libgd (checked the config.log file). Could you point me to a resource that
can help me in understanding what packages need to be installed for proper
functioning of the current version? I do not want to do an update of my
entire system for several reasons. Would rather install individual
dependencies. Or is there a way to configure gnuplot so that pngcairo, etc
are made available? I do not recall doing anything special for my earlier
version of gnuplot...or maybe since it's been there since a very long time
in my system I have forgotten.

Debostuti

On Thu, May 18, 2023 at 9:44 PM Ethan Merritt <eam...@gm...> wrote:

> On Thursday, 18 May 2023 05:08:36 PDT Debostuti Ghoshdastidar wrote:
>
> >
> > However the bars are not getting colored according to col i+1, although
> the
> > key is (please see attached image). Could you please help me understand
> > whats going wrong?
> >
> > Debostuti
>
> I am guessing that your gnuplot executable is not current.
> There was a bug that affected this specific kind of plot, fixed in
> the September 2022 gnuplot release (version 5.4 patchlevel 5).
> I should have remembered that and mentioned it before.
>
> The current release level is 5.4.6, with 5.4.7 queued for later
> this month.
>
>         Ethan
>
>
> > Hi Ethan
> >
> > I followed your suggestions and made little tweaks to the script as
> follows:
> > -----------------------------------------------
> > set lt 101 lw 4 lc "red"
> > set lt 102 lw 4 lc "pink"
> > set lt 103 lw 4 lc "dark-blue"
> > set lt 104 lw 4 lc "blue"
> > set lt 105 lw 4 lc "royalblue"
> > .
> > .
> > .
> > set lt 120 lw 4 lc 'white' #*set 20 custom colors*
> > set style data histogram
> > set style histogram rows
> > set style fill solid border lc "black"
> > #set style line 2 lc  'black' lt 1 lw 2
> > set boxwidth 0.80
> > alphabet = "GAVLIPWFYSTCMNQHKRDE"
> >     color(category) = 100 + strstrt(alphabet, category)
> >     plot for [i=3:10:2] 'test.dat' using
> > (column(i)):(color(strcol(i+1))):xtic(1):x2tic(2) lc variable
> >
> -----------------------------------------------------------------------------------------------------------------------
> >
> > for the following test data (the real data will have all 20 amino acids
> in
> > a row):
> >
> > test.dat:
> > 1       Gly     0.03    G       0       A       0       V       0       I
> > 2       Met     0.03    M       0       A       0       V       0       I
> > 3       Arg     0.94    R       0       A       0       V       0       I
> > 5       Arg     0.94    R       0       A       0       V       0       I
> > 6       Tyr     0.97    Y       0       A       0       V       0       I
> > 7       Val     0.03    V       0       G       0       A       0       I
> > 8       Cys     0.88    C       0       G       0       A       0       V
> > 9       Glu     1.00    E       0       G       0       A       0       V
> > 11      Pro     0.03    R       0.09    P       0       A       0       G
> > 12      Ser     0.58    S       0       G       0       A       0       V
> > 13      His     0.03    S       0.15    R       0.21    A       0.61    H
> > 14      Gly     0.03    A       0.06    S       0.61    G       0       V
> > 15      Gly     0.06    S       0.12    A       0.21    G       0       V
> >
> >
> --------------------------------------------------------------------------------------------------
> >
> >
> >
> > On Thu, May 18, 2023 at 2:10 AM Ethan Merritt <eam...@gm...>
> wrote:
> >
> > > On Wednesday, 17 May 2023 09:47:43 PDT Debostuti Ghoshdastidar wrote:
> > > > Hi Ethan
> > > >
> > > > Thank you for the solution. Each Y column will be a score for a
> > > particular
> > > > amino acid. The Y columns will be in increasing order of the score
> value.
> > > > Each Z will always correspond to a particular color we use to
> represent
> > > the
> > > > amino acid that will be in the corresponding y column. We will use a
> set
> > > of
> > > > 20 colors to represent 20 amino Acids.
> > >
> > > Ah. In that case you probably don't want "lc palette", which assumes
> > > a continuous range of values and a smoothly-varying color palette.
> > > Instead define 20 distinct colors in a convenient range of linetypes,
> > > say lt 101 through lt 120, indexed by the value in column Zn.
> > > That matches this demo from the on-line collection very closely:
> > >
> > >     https://gnuplot.sourceforge.net/demo_5.4/histogram_colors.html
> > >
> > > If you have the one letter code in the Z column (P for Proline, etc)
> > > you could use pieces of the demo code directly.
> > >
> > >     alphabet = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
> > >     color(category) = 100 + strstrt(alphabet, category)
> > >     plot for [i=2:*:2] DATA using (column(i)):(color(strcol(i+1)) lc
> > > variable
> > >
> > > Otherwise you could define a slightly more complicated color() function
> > > that
> > > for instance maps from the 3-letter code, e.g.
> > >     array AA = ["Ala", "Cys", "Asp", "Glu", ...]
> > >     color(category) = 100 + sum [1:20] (category eq AA[i]) ? i : 0)
> > >
> > > cheers,
> > >
> > >         Ethan
> > > >
> > > > Debostuti
> > > >
> > > > On Wed, 17 May 2023, 21:36 Ethan Merritt, <eam...@gm...>
> wrote:
> > > >
> > > > > On Tue, May 16, 2023 at 9:52 PM Debostuti Ghoshdastidar <
> > > > > deb...@gm...> wrote:
> > > > >
> > > > >> Dear Users
> > > > >>
> > > > >> I have data in the following format
> > > > >>
> > > > >> X1 Y1 Z1 Y2 Z2...Yn Zn
> > > > >> X2 Y1 Z1 Y2 Z2...Yn Zn
> > > > >> .
> > > > >> .
> > > > >> .
> > > > >> Xn Y1 Z1 Y2 Z2...Yn Zn
> > > > >>
> > > > >> I would like to draw a rowstacked histogram (which I know how to)
> > > where
> > > > >> the
> > > > >> colors for the bars are specified by the corresponding Z columns.
> For
> > > > >> example Y1 should be colored based on Z1 value, Y2 based on Z2
> value,
> > > and
> > > > >> so on...mush like the use of cbrange. Could you suggest how I can
> > > > >> accomplish this?
> > > > >>
> > > > >
> > > > > set style data histogram
> > > > > set style histogram rows
> > > > > set style fill solid border lc "black"
> > > > > set cbrange [0:<something>]
> > > > > plot for [i=2:*:2] DATA using (column(i)):(column(i+1)) lc palette
> > > title
> > > > > "???"
> > > > >
> > > > > It is not clear to me how you would label or construct a legend
> > > describing
> > > > > this in a way that is easy for the viewer to understand.  Perhaps
> > > > > superimpose
> > > > > a label on top of each individual box?
> > > > >
> > > > >     Ethan
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >> Thank you
> > > > >>
> > > > >> Debostuti
> > > > >> --
> > > > >> Debostuti Ghosh Dastidar
> > > > >> Research Associate
> > > > >> PDB-India
> > > > >>
> > > > >> _______________________________________________
> > > > >> gnuplot-info mailing list
> > > > >> gnu...@li...
> > > > >> Membership management via:
> > > > >> https://lists.sourceforge.net/lists/listinfo/gnuplot-info
> > > > >>
> > > > >
> > > >
> > >
> > >
> > >
> > >
> > >
> >
> >
>
>
>
>
>

-- 
Debostuti Ghosh Dastidar
Research Associate
PDB-India

Re: [Gnuplot-info] Regarding plotting stacked histograms with custom colors

[Ethan M. <eam...@gm...>]

On Thursday, 18 May 2023 05:08:36 PDT Debostuti Ghoshdastidar wrote:

> 
> However the bars are not getting colored according to col i+1, although the
> key is (please see attached image). Could you please help me understand
> whats going wrong?
> 
> Debostuti

I am guessing that your gnuplot executable is not current.
There was a bug that affected this specific kind of plot, fixed in
the September 2022 gnuplot release (version 5.4 patchlevel 5).
I should have remembered that and mentioned it before.

The current release level is 5.4.6, with 5.4.7 queued for later
this month.

	Ethan


> Hi Ethan
> 
> I followed your suggestions and made little tweaks to the script as follows:
> -----------------------------------------------
> set lt 101 lw 4 lc "red"
> set lt 102 lw 4 lc "pink"
> set lt 103 lw 4 lc "dark-blue"
> set lt 104 lw 4 lc "blue"
> set lt 105 lw 4 lc "royalblue"
> .
> .
> .
> set lt 120 lw 4 lc 'white' #*set 20 custom colors*
> set style data histogram
> set style histogram rows
> set style fill solid border lc "black"
> #set style line 2 lc  'black' lt 1 lw 2
> set boxwidth 0.80
> alphabet = "GAVLIPWFYSTCMNQHKRDE"
>     color(category) = 100 + strstrt(alphabet, category)
>     plot for [i=3:10:2] 'test.dat' using
> (column(i)):(color(strcol(i+1))):xtic(1):x2tic(2) lc variable
> -----------------------------------------------------------------------------------------------------------------------
> 
> for the following test data (the real data will have all 20 amino acids in
> a row):
> 
> test.dat:
> 1       Gly     0.03    G       0       A       0       V       0       I
> 2       Met     0.03    M       0       A       0       V       0       I
> 3       Arg     0.94    R       0       A       0       V       0       I
> 5       Arg     0.94    R       0       A       0       V       0       I
> 6       Tyr     0.97    Y       0       A       0       V       0       I
> 7       Val     0.03    V       0       G       0       A       0       I
> 8       Cys     0.88    C       0       G       0       A       0       V
> 9       Glu     1.00    E       0       G       0       A       0       V
> 11      Pro     0.03    R       0.09    P       0       A       0       G
> 12      Ser     0.58    S       0       G       0       A       0       V
> 13      His     0.03    S       0.15    R       0.21    A       0.61    H
> 14      Gly     0.03    A       0.06    S       0.61    G       0       V
> 15      Gly     0.06    S       0.12    A       0.21    G       0       V
> 
> --------------------------------------------------------------------------------------------------
> 
> 
> 
> On Thu, May 18, 2023 at 2:10 AM Ethan Merritt <eam...@gm...> wrote:
> 
> > On Wednesday, 17 May 2023 09:47:43 PDT Debostuti Ghoshdastidar wrote:
> > > Hi Ethan
> > >
> > > Thank you for the solution. Each Y column will be a score for a
> > particular
> > > amino acid. The Y columns will be in increasing order of the score value.
> > > Each Z will always correspond to a particular color we use to represent
> > the
> > > amino acid that will be in the corresponding y column. We will use a set
> > of
> > > 20 colors to represent 20 amino Acids.
> >
> > Ah. In that case you probably don't want "lc palette", which assumes
> > a continuous range of values and a smoothly-varying color palette.
> > Instead define 20 distinct colors in a convenient range of linetypes,
> > say lt 101 through lt 120, indexed by the value in column Zn.
> > That matches this demo from the on-line collection very closely:
> >
> >     https://gnuplot.sourceforge.net/demo_5.4/histogram_colors.html
> >
> > If you have the one letter code in the Z column (P for Proline, etc)
> > you could use pieces of the demo code directly.
> >
> >     alphabet = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
> >     color(category) = 100 + strstrt(alphabet, category)
> >     plot for [i=2:*:2] DATA using (column(i)):(color(strcol(i+1)) lc
> > variable
> >
> > Otherwise you could define a slightly more complicated color() function
> > that
> > for instance maps from the 3-letter code, e.g.
> >     array AA = ["Ala", "Cys", "Asp", "Glu", ...]
> >     color(category) = 100 + sum [1:20] (category eq AA[i]) ? i : 0)
> >
> > cheers,
> >
> >         Ethan
> > >
> > > Debostuti
> > >
> > > On Wed, 17 May 2023, 21:36 Ethan Merritt, <eam...@gm...> wrote:
> > >
> > > > On Tue, May 16, 2023 at 9:52 PM Debostuti Ghoshdastidar <
> > > > deb...@gm...> wrote:
> > > >
> > > >> Dear Users
> > > >>
> > > >> I have data in the following format
> > > >>
> > > >> X1 Y1 Z1 Y2 Z2...Yn Zn
> > > >> X2 Y1 Z1 Y2 Z2...Yn Zn
> > > >> .
> > > >> .
> > > >> .
> > > >> Xn Y1 Z1 Y2 Z2...Yn Zn
> > > >>
> > > >> I would like to draw a rowstacked histogram (which I know how to)
> > where
> > > >> the
> > > >> colors for the bars are specified by the corresponding Z columns. For
> > > >> example Y1 should be colored based on Z1 value, Y2 based on Z2 value,
> > and
> > > >> so on...mush like the use of cbrange. Could you suggest how I can
> > > >> accomplish this?
> > > >>
> > > >
> > > > set style data histogram
> > > > set style histogram rows
> > > > set style fill solid border lc "black"
> > > > set cbrange [0:<something>]
> > > > plot for [i=2:*:2] DATA using (column(i)):(column(i+1)) lc palette
> > title
> > > > "???"
> > > >
> > > > It is not clear to me how you would label or construct a legend
> > describing
> > > > this in a way that is easy for the viewer to understand.  Perhaps
> > > > superimpose
> > > > a label on top of each individual box?
> > > >
> > > >     Ethan
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >> Thank you
> > > >>
> > > >> Debostuti
> > > >> --
> > > >> Debostuti Ghosh Dastidar
> > > >> Research Associate
> > > >> PDB-India
> > > >>
> > > >> _______________________________________________
> > > >> gnuplot-info mailing list
> > > >> gnu...@li...
> > > >> Membership management via:
> > > >> https://lists.sourceforge.net/lists/listinfo/gnuplot-info
> > > >>
> > > >
> > >
> >
> >
> >
> >
> >
> 
> 

Re: [Gnuplot-info] Regarding plotting stacked histograms with custom colors

[Debostuti G. <deb...@gm...>]

Hi Ethan

I followed your suggestions and made little tweaks to the script as follows:
-----------------------------------------------
set lt 101 lw 4 lc "red"
set lt 102 lw 4 lc "pink"
set lt 103 lw 4 lc "dark-blue"
set lt 104 lw 4 lc "blue"
set lt 105 lw 4 lc "royalblue"
.
.
.
set lt 120 lw 4 lc 'white' #*set 20 custom colors*
set style data histogram
set style histogram rows
set style fill solid border lc "black"
#set style line 2 lc  'black' lt 1 lw 2
set boxwidth 0.80
alphabet = "GAVLIPWFYSTCMNQHKRDE"
    color(category) = 100 + strstrt(alphabet, category)
    plot for [i=3:10:2] 'test.dat' using
(column(i)):(color(strcol(i+1))):xtic(1):x2tic(2) lc variable
-----------------------------------------------------------------------------------------------------------------------

for the following test data (the real data will have all 20 amino acids in
a row):

test.dat:
1       Gly     0.03    G       0       A       0       V       0       I
2       Met     0.03    M       0       A       0       V       0       I
3       Arg     0.94    R       0       A       0       V       0       I
5       Arg     0.94    R       0       A       0       V       0       I
6       Tyr     0.97    Y       0       A       0       V       0       I
7       Val     0.03    V       0       G       0       A       0       I
8       Cys     0.88    C       0       G       0       A       0       V
9       Glu     1.00    E       0       G       0       A       0       V
11      Pro     0.03    R       0.09    P       0       A       0       G
12      Ser     0.58    S       0       G       0       A       0       V
13      His     0.03    S       0.15    R       0.21    A       0.61    H
14      Gly     0.03    A       0.06    S       0.61    G       0       V
15      Gly     0.06    S       0.12    A       0.21    G       0       V

--------------------------------------------------------------------------------------------------

However the bars are not getting colored according to col i+1, although the
key is (please see attached image). Could you please help me understand
whats going wrong?

Debostuti



On Thu, May 18, 2023 at 2:10 AM Ethan Merritt <eam...@gm...> wrote:

> On Wednesday, 17 May 2023 09:47:43 PDT Debostuti Ghoshdastidar wrote:
> > Hi Ethan
> >
> > Thank you for the solution. Each Y column will be a score for a
> particular
> > amino acid. The Y columns will be in increasing order of the score value.
> > Each Z will always correspond to a particular color we use to represent
> the
> > amino acid that will be in the corresponding y column. We will use a set
> of
> > 20 colors to represent 20 amino Acids.
>
> Ah. In that case you probably don't want "lc palette", which assumes
> a continuous range of values and a smoothly-varying color palette.
> Instead define 20 distinct colors in a convenient range of linetypes,
> say lt 101 through lt 120, indexed by the value in column Zn.
> That matches this demo from the on-line collection very closely:
>
>     https://gnuplot.sourceforge.net/demo_5.4/histogram_colors.html
>
> If you have the one letter code in the Z column (P for Proline, etc)
> you could use pieces of the demo code directly.
>
>     alphabet = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
>     color(category) = 100 + strstrt(alphabet, category)
>     plot for [i=2:*:2] DATA using (column(i)):(color(strcol(i+1)) lc
> variable
>
> Otherwise you could define a slightly more complicated color() function
> that
> for instance maps from the 3-letter code, e.g.
>     array AA = ["Ala", "Cys", "Asp", "Glu", ...]
>     color(category) = 100 + sum [1:20] (category eq AA[i]) ? i : 0)
>
> cheers,
>
>         Ethan
> >
> > Debostuti
> >
> > On Wed, 17 May 2023, 21:36 Ethan Merritt, <eam...@gm...> wrote:
> >
> > > On Tue, May 16, 2023 at 9:52 PM Debostuti Ghoshdastidar <
> > > deb...@gm...> wrote:
> > >
> > >> Dear Users
> > >>
> > >> I have data in the following format
> > >>
> > >> X1 Y1 Z1 Y2 Z2...Yn Zn
> > >> X2 Y1 Z1 Y2 Z2...Yn Zn
> > >> .
> > >> .
> > >> .
> > >> Xn Y1 Z1 Y2 Z2...Yn Zn
> > >>
> > >> I would like to draw a rowstacked histogram (which I know how to)
> where
> > >> the
> > >> colors for the bars are specified by the corresponding Z columns. For
> > >> example Y1 should be colored based on Z1 value, Y2 based on Z2 value,
> and
> > >> so on...mush like the use of cbrange. Could you suggest how I can
> > >> accomplish this?
> > >>
> > >
> > > set style data histogram
> > > set style histogram rows
> > > set style fill solid border lc "black"
> > > set cbrange [0:<something>]
> > > plot for [i=2:*:2] DATA using (column(i)):(column(i+1)) lc palette
> title
> > > "???"
> > >
> > > It is not clear to me how you would label or construct a legend
> describing
> > > this in a way that is easy for the viewer to understand.  Perhaps
> > > superimpose
> > > a label on top of each individual box?
> > >
> > >     Ethan
> > >
> > >
> > >
> > >
> > >
> > >> Thank you
> > >>
> > >> Debostuti
> > >> --
> > >> Debostuti Ghosh Dastidar
> > >> Research Associate
> > >> PDB-India
> > >>
> > >> _______________________________________________
> > >> gnuplot-info mailing list
> > >> gnu...@li...
> > >> Membership management via:
> > >> https://lists.sourceforge.net/lists/listinfo/gnuplot-info
> > >>
> > >
> >
>
>
>
>
>

-- 
Debostuti Ghosh Dastidar
Research Associate
PDB-India

Re: [Gnuplot-info] Regarding plotting stacked histograms with custom colors

[Debostuti G. <deb...@gm...>]

That's great!! Thanks a lot Ethan...

On Thu, 18 May 2023, 02:10 Ethan Merritt, <eam...@gm...> wrote:

> On Wednesday, 17 May 2023 09:47:43 PDT Debostuti Ghoshdastidar wrote:
> > Hi Ethan
> >
> > Thank you for the solution. Each Y column will be a score for a
> particular
> > amino acid. The Y columns will be in increasing order of the score value.
> > Each Z will always correspond to a particular color we use to represent
> the
> > amino acid that will be in the corresponding y column. We will use a set
> of
> > 20 colors to represent 20 amino Acids.
>
> Ah. In that case you probably don't want "lc palette", which assumes
> a continuous range of values and a smoothly-varying color palette.
> Instead define 20 distinct colors in a convenient range of linetypes,
> say lt 101 through lt 120, indexed by the value in column Zn.
> That matches this demo from the on-line collection very closely:
>
>     https://gnuplot.sourceforge.net/demo_5.4/histogram_colors.html
>
> If you have the one letter code in the Z column (P for Proline, etc)
> you could use pieces of the demo code directly.
>
>     alphabet = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
>     color(category) = 100 + strstrt(alphabet, category)
>     plot for [i=2:*:2] DATA using (column(i)):(color(strcol(i+1)) lc
> variable
>
> Otherwise you could define a slightly more complicated color() function
> that
> for instance maps from the 3-letter code, e.g.
>     array AA = ["Ala", "Cys", "Asp", "Glu", ...]
>     color(category) = 100 + sum [1:20] (category eq AA[i]) ? i : 0)
>
> cheers,
>
>         Ethan
> >
> > Debostuti
> >
> > On Wed, 17 May 2023, 21:36 Ethan Merritt, <eam...@gm...> wrote:
> >
> > > On Tue, May 16, 2023 at 9:52 PM Debostuti Ghoshdastidar <
> > > deb...@gm...> wrote:
> > >
> > >> Dear Users
> > >>
> > >> I have data in the following format
> > >>
> > >> X1 Y1 Z1 Y2 Z2...Yn Zn
> > >> X2 Y1 Z1 Y2 Z2...Yn Zn
> > >> .
> > >> .
> > >> .
> > >> Xn Y1 Z1 Y2 Z2...Yn Zn
> > >>
> > >> I would like to draw a rowstacked histogram (which I know how to)
> where
> > >> the
> > >> colors for the bars are specified by the corresponding Z columns. For
> > >> example Y1 should be colored based on Z1 value, Y2 based on Z2 value,
> and
> > >> so on...mush like the use of cbrange. Could you suggest how I can
> > >> accomplish this?
> > >>
> > >
> > > set style data histogram
> > > set style histogram rows
> > > set style fill solid border lc "black"
> > > set cbrange [0:<something>]
> > > plot for [i=2:*:2] DATA using (column(i)):(column(i+1)) lc palette
> title
> > > "???"
> > >
> > > It is not clear to me how you would label or construct a legend
> describing
> > > this in a way that is easy for the viewer to understand.  Perhaps
> > > superimpose
> > > a label on top of each individual box?
> > >
> > >     Ethan
> > >
> > >
> > >
> > >
> > >
> > >> Thank you
> > >>
> > >> Debostuti
> > >> --
> > >> Debostuti Ghosh Dastidar
> > >> Research Associate
> > >> PDB-India
> > >>
> > >> _______________________________________________
> > >> gnuplot-info mailing list
> > >> gnu...@li...
> > >> Membership management via:
> > >> https://lists.sourceforge.net/lists/listinfo/gnuplot-info
> > >>
> > >
> >
>
>
>
>
>

Re: [Gnuplot-info] Regarding plotting stacked histograms with custom colors

[Ethan M. <eam...@gm...>]

On Wednesday, 17 May 2023 09:47:43 PDT Debostuti Ghoshdastidar wrote:
> Hi Ethan
> 
> Thank you for the solution. Each Y column will be a score for a particular
> amino acid. The Y columns will be in increasing order of the score value.
> Each Z will always correspond to a particular color we use to represent the
> amino acid that will be in the corresponding y column. We will use a set of
> 20 colors to represent 20 amino Acids.

Ah. In that case you probably don't want "lc palette", which assumes
a continuous range of values and a smoothly-varying color palette.
Instead define 20 distinct colors in a convenient range of linetypes,
say lt 101 through lt 120, indexed by the value in column Zn.
That matches this demo from the on-line collection very closely:

    https://gnuplot.sourceforge.net/demo_5.4/histogram_colors.html

If you have the one letter code in the Z column (P for Proline, etc)
you could use pieces of the demo code directly.

    alphabet = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
    color(category) = 100 + strstrt(alphabet, category)
    plot for [i=2:*:2] DATA using (column(i)):(color(strcol(i+1)) lc variable

Otherwise you could define a slightly more complicated color() function that
for instance maps from the 3-letter code, e.g.
    array AA = ["Ala", "Cys", "Asp", "Glu", ...]
    color(category) = 100 + sum [1:20] (category eq AA[i]) ? i : 0) 

cheers,

	Ethan
> 
> Debostuti
> 
> On Wed, 17 May 2023, 21:36 Ethan Merritt, <eam...@gm...> wrote:
> 
> > On Tue, May 16, 2023 at 9:52 PM Debostuti Ghoshdastidar <
> > deb...@gm...> wrote:
> >
> >> Dear Users
> >>
> >> I have data in the following format
> >>
> >> X1 Y1 Z1 Y2 Z2...Yn Zn
> >> X2 Y1 Z1 Y2 Z2...Yn Zn
> >> .
> >> .
> >> .
> >> Xn Y1 Z1 Y2 Z2...Yn Zn
> >>
> >> I would like to draw a rowstacked histogram (which I know how to) where
> >> the
> >> colors for the bars are specified by the corresponding Z columns. For
> >> example Y1 should be colored based on Z1 value, Y2 based on Z2 value, and
> >> so on...mush like the use of cbrange. Could you suggest how I can
> >> accomplish this?
> >>
> >
> > set style data histogram
> > set style histogram rows
> > set style fill solid border lc "black"
> > set cbrange [0:<something>]
> > plot for [i=2:*:2] DATA using (column(i)):(column(i+1)) lc palette title
> > "???"
> >
> > It is not clear to me how you would label or construct a legend describing
> > this in a way that is easy for the viewer to understand.  Perhaps
> > superimpose
> > a label on top of each individual box?
> >
> >     Ethan
> >
> >
> >
> >
> >
> >> Thank you
> >>
> >> Debostuti
> >> --
> >> Debostuti Ghosh Dastidar
> >> Research Associate
> >> PDB-India
> >>
> >> _______________________________________________
> >> gnuplot-info mailing list
> >> gnu...@li...
> >> Membership management via:
> >> https://lists.sourceforge.net/lists/listinfo/gnuplot-info
> >>
> >
> 

Re: [Gnuplot-info] Regarding plotting stacked histograms with custom colors

[Debostuti G. <deb...@gm...>]

Hi Ethan

Thank you for the solution. Each Y column will be a score for a particular
amino acid. The Y columns will be in increasing order of the score value.
Each Z will always correspond to a particular color we use to represent the
amino acid that will be in the corresponding y column. We will use a set of
20 colors to represent 20 amino Acids.

Debostuti

On Wed, 17 May 2023, 21:36 Ethan Merritt, <eam...@gm...> wrote:

> On Tue, May 16, 2023 at 9:52 PM Debostuti Ghoshdastidar <
> deb...@gm...> wrote:
>
>> Dear Users
>>
>> I have data in the following format
>>
>> X1 Y1 Z1 Y2 Z2...Yn Zn
>> X2 Y1 Z1 Y2 Z2...Yn Zn
>> .
>> .
>> .
>> Xn Y1 Z1 Y2 Z2...Yn Zn
>>
>> I would like to draw a rowstacked histogram (which I know how to) where
>> the
>> colors for the bars are specified by the corresponding Z columns. For
>> example Y1 should be colored based on Z1 value, Y2 based on Z2 value, and
>> so on...mush like the use of cbrange. Could you suggest how I can
>> accomplish this?
>>
>
> set style data histogram
> set style histogram rows
> set style fill solid border lc "black"
> set cbrange [0:<something>]
> plot for [i=2:*:2] DATA using (column(i)):(column(i+1)) lc palette title
> "???"
>
> It is not clear to me how you would label or construct a legend describing
> this in a way that is easy for the viewer to understand.  Perhaps
> superimpose
> a label on top of each individual box?
>
>     Ethan
>
>
>
>
>
>> Thank you
>>
>> Debostuti
>> --
>> Debostuti Ghosh Dastidar
>> Research Associate
>> PDB-India
>>
>> _______________________________________________
>> gnuplot-info mailing list
>> gnu...@li...
>> Membership management via:
>> https://lists.sourceforge.net/lists/listinfo/gnuplot-info
>>
>

Re: [Gnuplot-info] Regarding plotting stacked histograms with custom colors

[Ethan M. <eam...@gm...>]

On Tue, May 16, 2023 at 9:52 PM Debostuti Ghoshdastidar <
deb...@gm...> wrote:

> Dear Users
>
> I have data in the following format
>
> X1 Y1 Z1 Y2 Z2...Yn Zn
> X2 Y1 Z1 Y2 Z2...Yn Zn
> .
> .
> .
> Xn Y1 Z1 Y2 Z2...Yn Zn
>
> I would like to draw a rowstacked histogram (which I know how to) where the
> colors for the bars are specified by the corresponding Z columns. For
> example Y1 should be colored based on Z1 value, Y2 based on Z2 value, and
> so on...mush like the use of cbrange. Could you suggest how I can
> accomplish this?
>

set style data histogram
set style histogram rows
set style fill solid border lc "black"
set cbrange [0:<something>]
plot for [i=2:*:2] DATA using (column(i)):(column(i+1)) lc palette title
"???"

It is not clear to me how you would label or construct a legend describing
this in a way that is easy for the viewer to understand.  Perhaps
superimpose
a label on top of each individual box?

    Ethan





> Thank you
>
> Debostuti
> --
> Debostuti Ghosh Dastidar
> Research Associate
> PDB-India
>
> _______________________________________________
> gnuplot-info mailing list
> gnu...@li...
> Membership management via:
> https://lists.sourceforge.net/lists/listinfo/gnuplot-info
>

[Gnuplot-info] Regarding plotting stacked histograms with custom colors

[Debostuti G. <deb...@gm...>]

Dear Users

I have data in the following format

X1 Y1 Z1 Y2 Z2...Yn Zn
X2 Y1 Z1 Y2 Z2...Yn Zn
.
.
.
Xn Y1 Z1 Y2 Z2...Yn Zn

I would like to draw a rowstacked histogram (which I know how to) where the
colors for the bars are specified by the corresponding Z columns. For
example Y1 should be colored based on Z1 value, Y2 based on Z2 value, and
so on...mush like the use of cbrange. Could you suggest how I can
accomplish this?

Thank you

Debostuti
-- 
Debostuti Ghosh Dastidar
Research Associate
PDB-India

Re: [Gnuplot-info] Plot with covering

[hchiPer <hc...@gm...>]

I don't see any attachment.
Maybe what you want is something like this?

	curve(x) = (x - 2)**2 + 3
	hline(x, y) = curve(x) > y ? 1/0 : y
	plot curve(x), hline(x, 25), hline(x, 75), hline(x, 125)

There are sometimes gaps visible at the ends of the horizontal lines, 
due to the sampling gnuplot makes within the xrange. You can make them 
almost invisible using

	set sample 1000 # or higher value, at your whish

Le 16/03/23 à 17:23, Patrick Dupre a écrit :
> Hello,
> 
> Is there is an option to plot a similar figure to that I am attaching?
> Typically, I have y values (horizontal lines) that I wish to limit the
> display by a curve
> 
> Thank
> 
> ===========================================================================
>   Patrick DUPRÉ                                 | | email: pd...@gm...
>   Laboratoire interdisciplinaire Carnot de Bourgogne
>   9 Avenue Alain Savary, BP 47870, 21078 DIJON Cedex FRANCE
>   Tel: +33 (0)380395988                    | | Room# D114A
> ===========================================================================
> 
> _______________________________________________
> gnuplot-info mailing list
> gnu...@li...
> Membership management via: https://lists.sourceforge.net/lists/listinfo/gnuplot-info

[Gnuplot-info] Plot with covering

[Patrick D. <pd...@gm...>]

Hello,

Is there is an option to plot a similar figure to that I am attaching?
Typically, I have y values (horizontal lines) that I wish to limit the
display by a curve

Thank

===========================================================================
 Patrick DUPRÉ                                 | | email: pd...@gm...
 Laboratoire interdisciplinaire Carnot de Bourgogne
 9 Avenue Alain Savary, BP 47870, 21078 DIJON Cedex FRANCE
 Tel: +33 (0)380395988                    | | Room# D114A
===========================================================================

[Gnuplot-info] gnuplot 5.4.6 windows binary packages are updated

[Tatsuro M. <tma...@ya...>]

gnuplot 5.4.6 windows binary packages are updated.

Fontconfig configuration files and related files are added to packages.
See
https://sourceforge.net/p/gnuplot/bugs/2512/

Site and files :
https://sourceforge.net/projects/gnuplot/files/gnuplot/5.4.6/
gp546-win64-mingw-2.7z
gp546-win64-mingw-2.exe

Tatsuro

Re: [Gnuplot-info] y2ticlabels

[Rusty T. <rus...@gm...>]

Excellent, thank you so very much. I compiled gnuplot 5.4.6 from 
https://sourceforge.net/projects/gnuplot/files/gnuplot/

On 2/26/23 23:39, hchiPer wrote:
> Hi,
>
> I was not more successful than you in using y2ticlabels with gnuplot 
> 5.4.4, but it works with gnuplot 5.4.5 (binaries compiled by myself).
>
> In this case, the y2labels are the values of the last column, with all 
> 4 decimals, which is not the result you expect.
>
> To have y2tics identical to ytics you can use:
>
> set xdata time
> set timefmt "%Y %b %d"
> set link y2
> set y2tics
> plot "/tmp/fn" using 1:4 with circles
>
> Laurent
>
>
>
>
> Le 25/02/23 à 22:49, Rusty Travis a écrit :
>> On 2/22/23 10:25, Rusty Travis wrote:
>>> Greetings, I'm having difficulty understanding the Help for 
>>> y2ticlabels and would appreciate your assistance.
>>>
>>> The file I'm plotting is quite simple:
>>> ...
>>> 2023 Jan 31  12.3038
>>> 2023 Jan 31  19.0661
>>> 2023 Feb 01  9.8433
>>> 2023 Feb 01  13.5022
>>> 2023 Feb 01  13.5083
>>> 2023 Feb 01  13.5318
>>> 2023 Feb 02  10.9849
>>> 2023 Feb 02  16.0072
>>> ...
>>>
>>> gnuplot> set xdata time
>>> gnuplot> set timefmt "%Y %b %d"
>>> gnuplot> plot "/tmp/fn" using 1:4 with circles
>>>
>>> However, trying to plot y2ticlabels on the right of the plot the 
>>> same as yticlabels seems difficult for me:
>>> gnuplot> plot "/tmp/fn" using 1:4:y2ticlabels(4) with circles
>>>
>>> ... results in the same display as above.
>>>
>>> How do I do a simple y2ticlabels so that the right ordinate axis is 
>>> displayed the same as the left axis, please?
>>>
>>>
>>>
>>
>> _______________________________________________
>> gnuplot-info mailing list
>> gnu...@li...
>> Membership management via: 
>> https://lists.sourceforge.net/lists/listinfo/gnuplot-info
>
>
> _______________________________________________
> gnuplot-info mailing list
> gnu...@li...
> Membership management via: 
> https://lists.sourceforge.net/lists/listinfo/gnuplot-info

Re: [Gnuplot-info] y2ticlabels

[hchiPer <hc...@gm...>]

Hi,

I was not more successful than you in using y2ticlabels with gnuplot 
5.4.4, but it works with gnuplot 5.4.5 (binaries compiled by myself).

In this case, the y2labels are the values of the last column, with all 4 
decimals, which is not the result you expect.

To have y2tics identical to ytics you can use:

set xdata time
set timefmt "%Y %b %d"
set link y2
set y2tics
plot "/tmp/fn" using 1:4 with circles

Laurent




Le 25/02/23 à 22:49, Rusty Travis a écrit :
> On 2/22/23 10:25, Rusty Travis wrote:
>> Greetings, I'm having difficulty understanding the Help for 
>> y2ticlabels and would appreciate your assistance.
>>
>> The file I'm plotting is quite simple:
>> ...
>> 2023 Jan 31  12.3038
>> 2023 Jan 31  19.0661
>> 2023 Feb 01  9.8433
>> 2023 Feb 01  13.5022
>> 2023 Feb 01  13.5083
>> 2023 Feb 01  13.5318
>> 2023 Feb 02  10.9849
>> 2023 Feb 02  16.0072
>> ...
>>
>> gnuplot> set xdata time
>> gnuplot> set timefmt "%Y %b %d"
>> gnuplot> plot "/tmp/fn" using 1:4 with circles
>>
>> However, trying to plot y2ticlabels on the right of the plot the same 
>> as yticlabels seems difficult for me:
>> gnuplot> plot "/tmp/fn" using 1:4:y2ticlabels(4) with circles
>>
>> ... results in the same display as above.
>>
>> How do I do a simple y2ticlabels so that the right ordinate axis is 
>> displayed the same as the left axis, please?
>>
>>
>>
> 
> _______________________________________________
> gnuplot-info mailing list
> gnu...@li...
> Membership management via: 
> https://lists.sourceforge.net/lists/listinfo/gnuplot-info

[Gnuplot-info] y2ticlabels

[Rusty T. <rus...@gm...>]

On 2/22/23 10:25, Rusty Travis wrote:
> Greetings, I'm having difficulty understanding the Help for 
> y2ticlabels and would appreciate your assistance.
>
> The file I'm plotting is quite simple:
> ...
> 2023 Jan 31  12.3038
> 2023 Jan 31  19.0661
> 2023 Feb 01  9.8433
> 2023 Feb 01  13.5022
> 2023 Feb 01  13.5083
> 2023 Feb 01  13.5318
> 2023 Feb 02  10.9849
> 2023 Feb 02  16.0072
> ...
>
> gnuplot> set xdata time
> gnuplot> set timefmt "%Y %b %d"
> gnuplot> plot "/tmp/fn" using 1:4 with circles
>
> However, trying to plot y2ticlabels on the right of the plot the same 
> as yticlabels seems difficult for me:
> gnuplot> plot "/tmp/fn" using 1:4:y2ticlabels(4) with circles
>
> ... results in the same display as above.
>
> How do I do a simple y2ticlabels so that the right ordinate axis is 
> displayed the same as the left axis, please?
>
>
>

Re: [Gnuplot-info] Trying to get boxes to have different colors for border and interior

[hchiPer <hc...@gm...>]

Hi,


The border properties can be set via the 'set fill style' command, where 
a linetype can be specified, so you have either to use the linetype 
which has the desired color as a default, or to redefine a particular 
linetype with the parameters you need.

I managed to get the expected result this way (I'm using gnuplot 5.4.4):

1) redefine one of the default linetypes with the desired parameters 
(line color, line width and even dash type)
      set linetype 1 lc 'blue'  lw 3 # dt 2

2) redefine the fill style in order that it uses the desired linetype 
(fillcolor can't be defined here)
      set style fill border lt 1 # or 'noborder' for a box without border

3) plot your data by specifying the fill color only (the border presence 
and properties depending of the parameters of the previous 'set style 
fill' command:
      plot '-' w boxes fs solid 0.1 fc 'red'
      0.2 0.3 0.1
      0.4 0.8 0.1
      0.6 0.4 0.1
      e


Laurent

Le 19/02/23 à 02:46, Robert Dodier a écrit :
> Hi, I'm working with boxes and I'd like to see if I can arrange for
> the line around the border of the box to be a different color than the
> interior of the box.
>
> I tried the following and I expected it to draw boxes which have a red
> interior and blue border:
>
> plot '-' with boxes fillstyle solid 0.1 linewidth 3 linecolor 'blue'
> fillcolor 'red'
> 0.2 0.3 0.1
> 0.4 0.8 0.1
> 0.6 0.4 0.1
> e
>
> but what I see is that the border is red as well as the interior.
>
> When I try it without fillcolor,
>
> plot '-' with boxes fillstyle solid 0.1 linewidth 3 linecolor 'blue'
> 0.2 0.3 0.1
> 0.4 0.8 0.1
> 0.6 0.4 0.1
> e
>
> the border is blue, as expected (and the interior is also blue, as
> expected, since the documentation for boxes says that the fillcolor
> defaults to the linecolor if not otherwise specified).
>
> Looking at the documentation for boxes, it doesn't seem to say
> anything about how the border color is determined, and also doesn't
> seem to say anything about the linecolor parameter, except to say
> fillcolor defaults to linecolor.
>
> Is the observed behavior for the first example to be expected, or is
> it a bug? Is there a way to get the border and the interior to be
> different colors?
>
> Thanks for any light you can shed on this, I appreciate your help.
>
> Robert Dodier
>
>
> _______________________________________________
> gnuplot-info mailing list
> gnu...@li...
> Membership management via: https://lists.sourceforge.net/lists/listinfo/gnuplot-info

Re: [Gnuplot-info] Trying to get boxes to have different colors for border and interior

[Robert D. <rob...@gm...>]

On Sat, Feb 18, 2023 at 6:46 PM Alan Corey <ala...@gm...> wrote:

> Can't you specify colors as RGB triplets?  I would think that would be
> preferable to color names.

Hi Alan, thanks for your reply. It turns out the behavior is the same
whether the colors are named or specified as rgb triplets.

best,

Robert

Re: [Gnuplot-info] Trying to get boxes to have different colors for border and interior

[Alan C. <ala...@gm...>]

Can't you specify colors as RGB triplets?  I would think that would be
preferable to color names.

On 2/18/23, Robert Dodier <rob...@gm...> wrote:
> Sorry, I forgot to mention. I am working w/ Gnuplot 5.4 patchlevel 2
> on Linux. Terminal type is 'qt'.
>
> On Sat, Feb 18, 2023 at 5:46 PM Robert Dodier <rob...@gm...>
> wrote:
>>
>> Hi, I'm working with boxes and I'd like to see if I can arrange for
>> the line around the border of the box to be a different color than the
>> interior of the box.
>>
>> I tried the following and I expected it to draw boxes which have a red
>> interior and blue border:
>>
>> plot '-' with boxes fillstyle solid 0.1 linewidth 3 linecolor 'blue'
>> fillcolor 'red'
>> 0.2 0.3 0.1
>> 0.4 0.8 0.1
>> 0.6 0.4 0.1
>> e
>>
>> but what I see is that the border is red as well as the interior.
>>
>> When I try it without fillcolor,
>>
>> plot '-' with boxes fillstyle solid 0.1 linewidth 3 linecolor 'blue'
>> 0.2 0.3 0.1
>> 0.4 0.8 0.1
>> 0.6 0.4 0.1
>> e
>>
>> the border is blue, as expected (and the interior is also blue, as
>> expected, since the documentation for boxes says that the fillcolor
>> defaults to the linecolor if not otherwise specified).
>>
>> Looking at the documentation for boxes, it doesn't seem to say
>> anything about how the border color is determined, and also doesn't
>> seem to say anything about the linecolor parameter, except to say
>> fillcolor defaults to linecolor.
>>
>> Is the observed behavior for the first example to be expected, or is
>> it a bug? Is there a way to get the border and the interior to be
>> different colors?
>>
>> Thanks for any light you can shed on this, I appreciate your help.
>>
>> Robert Dodier
>
>
> _______________________________________________
> gnuplot-info mailing list
> gnu...@li...
> Membership management via:
> https://lists.sourceforge.net/lists/listinfo/gnuplot-info
>


-- 
-------------
Education is contagious.

Re: [Gnuplot-info] Trying to get boxes to have different colors for border and interior

[Robert D. <rob...@gm...>]

Sorry, I forgot to mention. I am working w/ Gnuplot 5.4 patchlevel 2
on Linux. Terminal type is 'qt'.

On Sat, Feb 18, 2023 at 5:46 PM Robert Dodier <rob...@gm...> wrote:
>
> Hi, I'm working with boxes and I'd like to see if I can arrange for
> the line around the border of the box to be a different color than the
> interior of the box.
>
> I tried the following and I expected it to draw boxes which have a red
> interior and blue border:
>
> plot '-' with boxes fillstyle solid 0.1 linewidth 3 linecolor 'blue'
> fillcolor 'red'
> 0.2 0.3 0.1
> 0.4 0.8 0.1
> 0.6 0.4 0.1
> e
>
> but what I see is that the border is red as well as the interior.
>
> When I try it without fillcolor,
>
> plot '-' with boxes fillstyle solid 0.1 linewidth 3 linecolor 'blue'
> 0.2 0.3 0.1
> 0.4 0.8 0.1
> 0.6 0.4 0.1
> e
>
> the border is blue, as expected (and the interior is also blue, as
> expected, since the documentation for boxes says that the fillcolor
> defaults to the linecolor if not otherwise specified).
>
> Looking at the documentation for boxes, it doesn't seem to say
> anything about how the border color is determined, and also doesn't
> seem to say anything about the linecolor parameter, except to say
> fillcolor defaults to linecolor.
>
> Is the observed behavior for the first example to be expected, or is
> it a bug? Is there a way to get the border and the interior to be
> different colors?
>
> Thanks for any light you can shed on this, I appreciate your help.
>
> Robert Dodier

[Gnuplot-info] Trying to get boxes to have different colors for border and interior

[Robert D. <rob...@gm...>]

Hi, I'm working with boxes and I'd like to see if I can arrange for
the line around the border of the box to be a different color than the
interior of the box.

I tried the following and I expected it to draw boxes which have a red
interior and blue border:

plot '-' with boxes fillstyle solid 0.1 linewidth 3 linecolor 'blue'
fillcolor 'red'
0.2 0.3 0.1
0.4 0.8 0.1
0.6 0.4 0.1
e

but what I see is that the border is red as well as the interior.

When I try it without fillcolor,

plot '-' with boxes fillstyle solid 0.1 linewidth 3 linecolor 'blue'
0.2 0.3 0.1
0.4 0.8 0.1
0.6 0.4 0.1
e

the border is blue, as expected (and the interior is also blue, as
expected, since the documentation for boxes says that the fillcolor
defaults to the linecolor if not otherwise specified).

Looking at the documentation for boxes, it doesn't seem to say
anything about how the border color is determined, and also doesn't
seem to say anything about the linecolor parameter, except to say
fillcolor defaults to linecolor.

Is the observed behavior for the first example to be expected, or is
it a bug? Is there a way to get the border and the interior to be
different colors?

Thanks for any light you can shed on this, I appreciate your help.

Robert Dodier

Re: [Gnuplot-info] Gnuplot is interpolating missing data

[hchiPer <hc...@gm...>]

Instead of an 'x', use 'Nan' or '1/0' (without the quotes) for 
unavailable data. You will get a gap in the graphs. See gnuplot manual 
under the section "Set datafile missing" for more info. There is a gap 
only where a data is invalid (e.g. NaN or 1/0). When a data is missing, 
the line is drawn anyway.

When gnuplot tries to compute the square root of the products, the 
result is indeed NaN when one of the values is the letter x, and this 
explains why there is a gap in your third line.


Le 18/12/22 à 22:27, Dave Horsfall a écrit :
> Gnuplot 5.4.3 on old MacBook Pro running High Sierra.
>
> I have many graphs which have missing data e.g. I was too lazy to record
> it on that day etc; these are indicated with "x" in the data file.
>
> Now, hie thee to http://www.horsfall.org/bike.pdf which is a graph of the
> progress (if any) that I am making on my exercise bike; the "brake" is the
> friction adjustment and does not seem to be in any particular unit.
>
> For laughs I compute an overall "effort" value which is the product of the
> distance and the brake (I was forced to rescale it by taking the square
> root to both make it fit and to make it look better, but that's irrelevant
> here).
>
> Notice how the recorded data is interpolated over the missing data, yet
> the computed plot shows a gap (the desired behaviour).
>
> (I do not use `set datafile missing "x"' anywhere)
>
> Oddly enough my other graphs such as www.horsfall.org/health.pdf etc are
> fine; the only difference is the linetype where I am using "lt rgb" and no
> "pt", whereas with my exercise bike I am using a specific point type to be
> compatible with both a colour inkjet and a mono laser.
>
> Is this expected behaviour?  It looks like a bug to me...
>
> Attached are the relevant *.dat and *.gp files.
>
> Thanks.
>
>   -- Dave
>
>
> _______________________________________________
> gnuplot-info mailing list
> gnu...@li...
> Membership management via: https://lists.sourceforge.net/lists/listinfo/gnuplot-info

[Gnuplot-info] Gnuplot is interpolating missing data

[Dave H. <da...@ho...>]

Gnuplot 5.4.3 on old MacBook Pro running High Sierra.

I have many graphs which have missing data e.g. I was too lazy to record 
it on that day etc; these are indicated with "x" in the data file.

Now, hie thee to http://www.horsfall.org/bike.pdf which is a graph of the 
progress (if any) that I am making on my exercise bike; the "brake" is the 
friction adjustment and does not seem to be in any particular unit.

For laughs I compute an overall "effort" value which is the product of the 
distance and the brake (I was forced to rescale it by taking the square 
root to both make it fit and to make it look better, but that's irrelevant 
here).

Notice how the recorded data is interpolated over the missing data, yet 
the computed plot shows a gap (the desired behaviour).

(I do not use `set datafile missing "x"' anywhere)

Oddly enough my other graphs such as www.horsfall.org/health.pdf etc are 
fine; the only difference is the linetype where I am using "lt rgb" and no 
"pt", whereas with my exercise bike I am using a specific point type to be 
compatible with both a colour inkjet and a mono laser.

Is this expected behaviour?  It looks like a bug to me...

Attached are the relevant *.dat and *.gp files.

Thanks.

 -- Dave

Re: [Gnuplot-info] Possible optimal way to handle extremum points problem

[Dan H. <dan...@ph...>]

On Tue, 29 Nov 2022, Ethan Merritt wrote:

> The program already does something of the sort to insure that a function
> is evaluated at x=0 even if regular subdivision of the x-axis range would
> not generate a sample point there.  If the location of the function's
> extrema are known in advance, those special values of x could be added
> to the sample set.  Unfortunately the program does not in general know
> how many or where such extrema are expected.

For that, you probably want a computer algebra system that can
interface with Gnuplot for plotting purposes (e.g. Maxima) rather than
Gnuplot on its own.

Re: [Gnuplot-info] invalid context

[theozh <th...@gm...>]

Apparently, for some reason (which I don't know) you cannot define an inline datablock within an if-else loop.
However, you can *print* to datablocks.

flag_satur=1

if (flag_satur==1) {
     set print $EXT_CEN
         print "0.0 1e10 1e-13"
         print "0.0 1e10 1"
     set print $EXT_p1
         print "5 1e10 1e-13"
         print "5 1e10 1"
     set print $EXT_m1
         print "-5 1e10 1e-13"
         print "-5 1e10 1"
     set print
     }
else {
     # For Linear absorption
     set print $EXT_CEN
         print "0.0 1e4 1e-11"
         print "0.0 1e4 1"
     set print $EXT_p1
         print "5 1e4 1e-11"
         print "5 1e4 1"
     set print $EXT_m1
         print "-5 1e4 1e-11"
         print "-5 1e4 1"
     set print
}


> Maybe I missed the point, there is a
>
> replot $EXT_CEN  w l lt 0 lw 1 dt "." lc rgb "gray"
>
>

Re: [Gnuplot-info] invalid context

[Patrick D. <pd...@gm...>]

Maybe I missed the point, there is a

replot $EXT_CEN  w l lt 0 lw 1 dt "." lc rgb "gray"


> Sent: Tuesday, December 06, 2022 at 4:49 PM
> From: "Patrick Dupre" <pd...@gm...>
> To: "theozh" <th...@gm...>
> Cc: gnu...@li...
> Subject: Re: [Gnuplot-info] invalid context
>
> Thank for the feedback,
> 
> Of course not,
> There is an else
> but the issue is with the if !
> 
> if (flag_satur==1) {
> $EXT_CEN << EOD
>  0.0 1e10 1e-13
>  0.0 1e10 1
> EOD
> $EXT_p1 << EOD
>  5 1e10 1e-13
>  5 1e10 1
> EOD
> $EXT_m1 << EOD
>  -5 1e10 1e-13
>  -5 1e10 1
> EOD
>  }
> else {
> # For Linear absorption
> $EXT_CEN << EOD
>  0.0 1e4 1e-11
>  0.0 1e4 1
> EOD
> $EXT_p1 << EOD
>  5 1e4 1e-11
>  5 1e4 1
> EOD
> $EXT_m1 << EOD
>  -5 1e4 1e-11
>  -5 1e4 1
> EOD
>  }

replot $EXT_CEN  w l lt 0 lw 1 dt "." lc rgb "gray"


> 
> ===========================================================================
>  Patrick DUPRÉ                                 | | email: pd...@gm...
>  Laboratoire interdisciplinaire Carnot de Bourgogne
>  9 Avenue Alain Savary, BP 47870, 21078 DIJON Cedex FRANCE
>  Tel: +33 (0)380395988                    | | Room# D114A
> ===========================================================================
> 
> 
> > Sent: Tuesday, December 06, 2022 at 4:27 PM
> > From: "theozh" <th...@gm...>
> > To: gnu...@li...
> > Subject: Re: [Gnuplot-info] invalid context
> >
> > What is the larger context of your script snippet?
> > 
> > What should happen if flag_satur!=0 ?
> > Then $EXT_CEN should become (or stay) undefined or become something else?
> > 
> > > Hello,
> > >
> > > It seems that I cannot do this
> > > if (flag_satur==1) {
> > > $EXT_CEN << EOD
> > >   0.0 1e10 1e-13
> > >   0.0 1e10 1
> > >   EOD
> > > }
> > >
> > 
> > 
> > _______________________________________________
> > gnuplot-info mailing list
> > gnu...@li...
> > Membership management via: https://lists.sourceforge.net/lists/listinfo/gnuplot-info
> >
> 
> 
> _______________________________________________
> gnuplot-info mailing list
> gnu...@li...
> Membership management via: https://lists.sourceforge.net/lists/listinfo/gnuplot-info
>

Re: [Gnuplot-info] invalid context

[Patrick D. <pd...@gm...>]

Thank for the feedback,

Of course not,
There is an else
but the issue is with the if !

if (flag_satur==1) {
$EXT_CEN << EOD
 0.0 1e10 1e-13
 0.0 1e10 1
EOD
$EXT_p1 << EOD
 5 1e10 1e-13
 5 1e10 1
EOD
$EXT_m1 << EOD
 -5 1e10 1e-13
 -5 1e10 1
EOD
 }
else {
# For Linear absorption
$EXT_CEN << EOD
 0.0 1e4 1e-11
 0.0 1e4 1
EOD
$EXT_p1 << EOD
 5 1e4 1e-11
 5 1e4 1
EOD
$EXT_m1 << EOD
 -5 1e4 1e-11
 -5 1e4 1
EOD
 }


===========================================================================
 Patrick DUPRÉ                                 | | email: pd...@gm...
 Laboratoire interdisciplinaire Carnot de Bourgogne
 9 Avenue Alain Savary, BP 47870, 21078 DIJON Cedex FRANCE
 Tel: +33 (0)380395988                    | | Room# D114A
===========================================================================


> Sent: Tuesday, December 06, 2022 at 4:27 PM
> From: "theozh" <th...@gm...>
> To: gnu...@li...
> Subject: Re: [Gnuplot-info] invalid context
>
> What is the larger context of your script snippet?
> 
> What should happen if flag_satur!=0 ?
> Then $EXT_CEN should become (or stay) undefined or become something else?
> 
> > Hello,
> >
> > It seems that I cannot do this
> > if (flag_satur==1) {
> > $EXT_CEN << EOD
> >   0.0 1e10 1e-13
> >   0.0 1e10 1
> >   EOD
> > }
> >
> 
> 
> _______________________________________________
> gnuplot-info mailing list
> gnu...@li...
> Membership management via: https://lists.sourceforge.net/lists/listinfo/gnuplot-info
>