#include <stdio.h>
#include <math.h>
-#include <gsl/gsl_histogram.h>
#include <gsl/gsl_randist.h>
#include <assert.h>
#include <output/charts/plot-hist.h>
-#include <output/charts/plot-chart.h>
-#include <output/chart-provider.h>
+#include <output/chart-item-provider.h>
-#include <data/variable.h>
#include <libpspp/cast.h>
-#include <libpspp/hash.h>
-#include <output/chart.h>
#include <math/histogram.h>
#include <math/moments.h>
#include "gettext.h"
#define _(msgid) gettext (msgid)
-static const struct chart_class histogram_chart_class;
-
-/* Write the legend of the chart */
-static void
-histogram_write_legend (cairo_t *cr, const struct chart_geometry *geom,
- double n, double mean, double stddev)
-{
- double y = geom->data_bottom;
- cairo_save (cr);
-
- if (n != SYSMIS)
- {
- char *buf = xasprintf ("N = %.2f", n);
- cairo_move_to (cr, geom->legend_left, y);
- chart_label (cr, 'l', 'b', geom->font_size, buf);
- y += geom->font_size * 1.5;
- free (buf);
- }
-
- if (mean != SYSMIS)
- {
- char *buf = xasprintf ("Mean = %.1f", mean);
- cairo_move_to (cr,geom->legend_left, y);
- chart_label (cr, 'l', 'b', geom->font_size, buf);
- y += geom->font_size * 1.5;
- free (buf);
- }
-
- if (stddev != SYSMIS)
- {
- char *buf = xasprintf ("Std. Dev = %.2f", stddev);
- cairo_move_to (cr, geom->legend_left, y);
- chart_label (cr, 'l', 'b', geom->font_size, buf);
- free (buf);
- }
-
- cairo_restore (cr);
-}
-
-static void
-hist_draw_bar (cairo_t *cr, const struct chart_geometry *geom,
- const gsl_histogram *h, int bar)
-{
- double upper;
- double lower;
- double height;
-
- const size_t bins = gsl_histogram_bins (h);
- const double x_pos = (geom->data_right - geom->data_left) * bar / (double) bins ;
- const double width = (geom->data_right - geom->data_left) / (double) bins ;
-
- assert ( 0 == gsl_histogram_get_range (h, bar, &lower, &upper));
-
- assert ( upper >= lower);
-
- height = gsl_histogram_get (h, bar) *
- (geom->data_top - geom->data_bottom) / gsl_histogram_max_val (h);
-
- cairo_rectangle (cr, geom->data_left + x_pos, geom->data_bottom,
- width, height);
- cairo_save (cr);
- cairo_set_source_rgb (cr,
- geom->fill_colour.red / 255.0,
- geom->fill_colour.green / 255.0,
- geom->fill_colour.blue / 255.0);
- cairo_fill_preserve (cr);
- cairo_restore (cr);
- cairo_stroke (cr);
-
- draw_tick (cr, geom, TICK_ABSCISSA,
- x_pos + width / 2.0, "%g", (upper + lower) / 2.0);
-}
-
-struct histogram_chart
- {
- struct chart chart;
- gsl_histogram *gsl_hist;
- char *label;
- double n;
- double mean;
- double stddev;
- bool show_normal;
- };
-
/* Plots a histogram of the data in HIST with the given LABEL.
Labels the histogram with each of N, MEAN, and STDDEV that is
not SYSMIS. If all three are not SYSMIS and SHOW_NORMAL is
true, also draws a normal curve on the histogram. */
-struct chart *
-histogram_chart_create (const struct histogram *hist, const char *label,
+struct chart_item *
+histogram_chart_create (const gsl_histogram *hist, const char *label,
double n, double mean, double stddev,
bool show_normal)
{
struct histogram_chart *h;
h = xmalloc (sizeof *h);
- chart_init (&h->chart, &histogram_chart_class);
- h->gsl_hist = hist->gsl_hist ? gsl_histogram_clone (hist->gsl_hist) : NULL;
- h->label = xstrdup (label);
+ chart_item_init (&h->chart_item, &histogram_chart_class, label);
+ h->gsl_hist = hist != NULL ? gsl_histogram_clone (hist) : NULL;
h->n = n;
h->mean = mean;
h->stddev = stddev;
h->show_normal = show_normal;
- return &h->chart;
-}
-
-static void
-histogram_chart_draw (const struct chart *chart, cairo_t *cr,
- struct chart_geometry *geom)
-{
- struct histogram_chart *h = UP_CAST (chart, struct histogram_chart, chart);
- int i;
- int bins;
-
- chart_write_title (cr, geom, _("HISTOGRAM"));
-
- chart_write_ylabel (cr, geom, _("Frequency"));
- chart_write_xlabel (cr, geom, h->label);
-
- if (h->gsl_hist == NULL)
- {
- /* Probably all values are SYSMIS. */
- return;
- }
-
- bins = gsl_histogram_bins (h->gsl_hist);
-
- chart_write_yscale (cr, geom, 0, gsl_histogram_max_val (h->gsl_hist), 5);
-
- for (i = 0; i < bins; i++)
- hist_draw_bar (cr, geom, h->gsl_hist, i);
-
- histogram_write_legend (cr, geom, h->n, h->mean, h->stddev);
-
- if (h->show_normal
- && h->n != SYSMIS && h->mean != SYSMIS && h->stddev != SYSMIS)
- {
- /* Draw the normal curve */
- double d;
- double x_min, x_max, not_used;
- double abscissa_scale;
- double ordinate_scale;
- double range;
-
- gsl_histogram_get_range (h->gsl_hist, 0, &x_min, ¬_used);
- range = not_used - x_min;
- gsl_histogram_get_range (h->gsl_hist, bins - 1, ¬_used, &x_max);
-
- abscissa_scale = (geom->data_right - geom->data_left) / (x_max - x_min);
- ordinate_scale = (geom->data_top - geom->data_bottom) /
- gsl_histogram_max_val (h->gsl_hist);
-
- cairo_move_to (cr, geom->data_left, geom->data_bottom);
- for (d = geom->data_left;
- d <= geom->data_right;
- d += (geom->data_right - geom->data_left) / 100.0)
- {
- const double x = (d - geom->data_left) / abscissa_scale + x_min;
- const double y = h->n * range *
- gsl_ran_gaussian_pdf (x - h->mean, h->stddev);
-
- cairo_line_to (cr, d, geom->data_bottom + y * ordinate_scale);
-
- }
- cairo_stroke (cr);
- }
+ return &h->chart_item;
}
-
static void
-histogram_chart_destroy (struct chart *chart)
+histogram_chart_destroy (struct chart_item *chart_item)
{
- struct histogram_chart *h = UP_CAST (chart, struct histogram_chart, chart);
+ struct histogram_chart *h = UP_CAST (chart_item, struct histogram_chart,
+ chart_item);
if (h->gsl_hist != NULL)
gsl_histogram_free (h->gsl_hist);
- free (h->label);
free (h);
}
-static const struct chart_class histogram_chart_class =
+const struct chart_item_class histogram_chart_class =
{
- histogram_chart_draw,
histogram_chart_destroy
};
projects / pspp-builds.git / blobdiff