pintos-os.org Git - pspp/blob - src/output/charts/plot-hist-cairo.c

   1 /* PSPP - a program for statistical analysis.
   2    Copyright (C) 2009, 2011, 2014 Free Software Foundation, Inc.
   3
   4    This program is free software: you can redistribute it and/or modify
   5    it under the terms of the GNU General Public License as published by
   6    the Free Software Foundation, either version 3 of the License, or
   7    (at your option) any later version.
   8
   9    This program is distributed in the hope that it will be useful,
  10    but WITHOUT ANY WARRANTY; without even the implied warranty of
  11    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12    GNU General Public License for more details.
  13
  14    You should have received a copy of the GNU General Public License
  15    along with this program.  If not, see <http://www.gnu.org/licenses/>. */
  16
  17 #include <config.h>
  18
  19 #include "output/charts/plot-hist.h"
  20
  21 #include <float.h>
  22 #include <gsl/gsl_randist.h>
  23
  24 #include "data/val-type.h"
  25 #include "output/cairo-chart.h"
  26
  27 #include "gl/xvasprintf.h"
  28
  29 #include "gettext.h"
  30 #define _(msgid) gettext (msgid)
  31
  32 /* Write the legend of the chart */
  33 static void
  34 histogram_write_legend (cairo_t *cr, const struct xrchart_geometry *geom,
  35                         double n, double mean, double stddev)
  36 {
  37   double y = geom->axis[SCALE_ORDINATE].data_min;
  38   cairo_save (cr);
  39
  40   if (n != SYSMIS)
  41     {
  42       char *buf = xasprintf (_("N = %.2f"), n);
  43       cairo_move_to (cr, geom->legend_left, y);
  44       xrchart_label (cr, 'l', 'b', geom->font_size, buf);
  45       y += geom->font_size * 1.5;
  46       free (buf);
  47     }
  48
  49   if (mean != SYSMIS)
  50     {
  51       char *buf = xasprintf (_("Mean = %.1f"), mean);
  52       cairo_move_to (cr,geom->legend_left, y);
  53       xrchart_label (cr, 'l', 'b', geom->font_size, buf);
  54       y += geom->font_size * 1.5;
  55       free (buf);
  56     }
  57
  58   if (stddev != SYSMIS)
  59     {
  60       char *buf = xasprintf (_("Std. Dev = %.2f"), stddev);
  61       cairo_move_to (cr, geom->legend_left, y);
  62       xrchart_label (cr, 'l', 'b', geom->font_size, buf);
  63       free (buf);
  64     }
  65
  66   cairo_restore (cr);
  67 }
  68
  69 static void
  70 hist_draw_bar (cairo_t *cr, const struct xrchart_geometry *geom,
  71                const gsl_histogram *h, int bar, bool label)
  72 {
  73   double upper;
  74   double lower;
  75   double height;
  76
  77   const size_t bins = gsl_histogram_bins (h);
  78
  79   const double x_pos =
  80     (geom->axis[SCALE_ABSCISSA].data_max - geom->axis[SCALE_ABSCISSA].data_min) *
  81     bar / (double) bins ;
  82
  83   const double width =
  84     (geom->axis[SCALE_ABSCISSA].data_max - geom->axis[SCALE_ABSCISSA].data_min) / (double) bins ;
  85
  86   assert ( 0 == gsl_histogram_get_range (h, bar, &lower, &upper));
  87
  88   assert ( upper >= lower);
  89
  90   height = geom->axis[SCALE_ORDINATE].scale * gsl_histogram_get (h, bar);
  91
  92   cairo_rectangle (cr,
  93                    geom->axis[SCALE_ABSCISSA].data_min + x_pos,
  94                    geom->axis[SCALE_ORDINATE].data_min,
  95                    width, height);
  96   cairo_save (cr);
  97   cairo_set_source_rgb (cr,
  98                         geom->fill_colour.red / 255.0,
  99                         geom->fill_colour.green / 255.0,
 100                         geom->fill_colour.blue / 255.0);
 101   cairo_fill_preserve (cr);
 102   cairo_restore (cr);
 103   cairo_stroke (cr);
 104
 105   if (label)
 106     draw_tick (cr, geom, SCALE_ABSCISSA, bins > 10,
 107                x_pos + width / 2.0, "%.*g",
 108                DBL_DIG + 1, (upper + lower) / 2.0);
 109 }
 110
 111 void
 112 xrchart_draw_histogram (const struct chart_item *chart_item, cairo_t *cr,
 113                         struct xrchart_geometry *geom)
 114 {
 115   struct histogram_chart *h = to_histogram_chart (chart_item);
 116   int i;
 117   int bins;
 118
 119   xrchart_write_title (cr, geom, _("HISTOGRAM"));
 120
 121   xrchart_write_ylabel (cr, geom, _("Frequency"));
 122   xrchart_write_xlabel (cr, geom, chart_item_get_title (chart_item));
 123
 124   if (h->gsl_hist == NULL)
 125     {
 126       /* Probably all values are SYSMIS. */
 127       return;
 128     }
 129
 130   bins = gsl_histogram_bins (h->gsl_hist);
 131
 132   xrchart_write_yscale (cr, geom, 0, gsl_histogram_max_val (h->gsl_hist), 5);
 133
 134   for (i = 0; i < bins; i++)
 135     {
 136       hist_draw_bar (cr, geom, h->gsl_hist, i, true);
 137     }
 138
 139   histogram_write_legend (cr, geom, h->n, h->mean, h->stddev);
 140
 141   if (h->show_normal
 142       && h->n != SYSMIS && h->mean != SYSMIS && h->stddev != SYSMIS)
 143     {
 144       /* Draw the normal curve */
 145       double d;
 146       double x_min, x_max, not_used;
 147       double abscissa_scale;
 148       double ordinate_scale;
 149       double range;
 150
 151       gsl_histogram_get_range (h->gsl_hist, 0, &x_min, &not_used);
 152       range = not_used - x_min;
 153       gsl_histogram_get_range (h->gsl_hist, bins - 1, &not_used, &x_max);
 154
 155       abscissa_scale = (geom->axis[SCALE_ABSCISSA].data_max - geom->axis[SCALE_ABSCISSA].data_min) / (x_max - x_min);
 156       ordinate_scale = (geom->axis[SCALE_ORDINATE].data_max - geom->axis[SCALE_ORDINATE].data_min) /
 157         gsl_histogram_max_val (h->gsl_hist);
 158
 159       cairo_move_to (cr, geom->axis[SCALE_ABSCISSA].data_min, geom->axis[SCALE_ORDINATE].data_min);
 160       for (d = geom->axis[SCALE_ABSCISSA].data_min;
 161            d <= geom->axis[SCALE_ABSCISSA].data_max;
 162            d += (geom->axis[SCALE_ABSCISSA].data_max - geom->axis[SCALE_ABSCISSA].data_min) / 100.0)
 163         {
 164           const double x = (d - geom->axis[SCALE_ABSCISSA].data_min) / abscissa_scale + x_min;
 165           const double y = h->n * range *
 166             gsl_ran_gaussian_pdf (x - h->mean, h->stddev);
 167
 168           cairo_line_to (cr, d, geom->axis[SCALE_ORDINATE].data_min  + y * ordinate_scale);
 169
 170         }
 171       cairo_stroke (cr);
 172     }
 173 }