From: Jason Stover <jhs@math.gcsu.edu>
Date: Thu, 31 May 2007 23:09:21 +0000 (+0000)
Subject: Initial version of interaction
X-Git-Tag: v0.6.0~466
X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=5715e8846c6e3a8dd539b2dc9c9c65640f84af76;p=pspp-builds.git

Initial version of interaction
---

diff --git a/src/math/ChangeLog b/src/math/ChangeLog
index e351fbbb..acc40388 100644
--- a/src/math/ChangeLog
+++ b/src/math/ChangeLog
@@ -1,3 +1,8 @@
+2007-05-31  Jason Stover  <jhs@math.gcsu.edu>
+
+	* interaction.c: New file.
+	* interaction.h : New file.
+
 Mon Feb  5 15:42:14 2007  Ben Pfaff  <blp@gnu.org>
 
 	* moments.c (moments_pass_two): Reduce number of multiplications.
diff --git a/src/math/automake.mk b/src/math/automake.mk
index 45543fc3..2a8f4408 100644
--- a/src/math/automake.mk
+++ b/src/math/automake.mk
@@ -13,8 +13,10 @@ src_math_libpspp_math_a_SOURCES = \
 	src/math/coefficient.c \
 	src/math/coefficient.h \
 	src/math/group.c  src/math/group.h \
-	src/math/histogram.c src/math/histogram.h \
 	src/math/group-proc.h \
+	src/math/histogram.c src/math/histogram.h \
+	src/math/interaction.c \
+	src/math/interaction.h \
 	src/math/levene.c \
 	src/math/levene.h \
 	src/math/moments.c  src/math/moments.h \
diff --git a/src/math/interaction.c b/src/math/interaction.c
new file mode 100644
index 00000000..be3d3edb
--- /dev/null
+++ b/src/math/interaction.c
@@ -0,0 +1,103 @@
+/* PSPP - Creates data structures to store interactions for
+   statistical routines.  
+
+   Copyright (C) 2007 Free Software Foundation, Inc.
+
+   This program is free software; you can redistribute it and/or
+   modify it under the terms of the GNU General Public License as
+   published by the Free Software Foundation; either version 2 of the
+   License, or (at your option) any later version.
+
+   This program is distributed in the hope that it will be useful, but
+   WITHOUT ANY WARRANTY; without even the implied warranty of
+   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+   General Public License for more details.
+
+   You should have received a copy of the GNU General Public License
+   along with this program; if not, write to the Free Software
+   Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
+   02110-1301, USA. */
+
+/*
+  An interaction is a gsl_vector containing a "product" of other
+  variables. The variables can be either categorical or numeric.
+  If the variables are all numeric, the interaction is just the
+  scalar product. If any of the variables are categorical, their
+  product is a vector containing 0's in all but one entry. This entry
+  is found by combining the vectors corresponding to the variables'
+  OBS_VALS member. If there are K categorical variables, each with
+  N_1, N_2, ..., N_K categories, then the interaction will have
+  N_1 * N_2 * N_3 *...* N_K - 1 entries.
+ */
+#include <assert.h>
+#include <libpspp/alloc.h>
+#include <gsl/gsl_math.h>
+#include <gsl/gsl_vector.h>
+#include <data/category.h>
+#include <data/variable.h>
+
+/*
+  Convert a list of values to a binary vector. The order of VALS must
+  correspond to the order of V.
+ */
+gsl_vector *
+get_interaction (union value **vals, const struct variable **v, size_t n_vars)
+{
+  gsl_vector *result = NULL;
+  size_t *subs = NULL;
+  size_t length = 1;
+  size_t i;
+  size_t j;
+  double tmp = 1.0;
+
+  assert (n_vars > 0);
+  for (i = 0; i < n_vars; i++)
+    {
+      if (var_is_alpha (v[i]))
+	{
+	  length *= cat_get_n_categories (v[i]);
+	}
+      else
+	{
+	  length = (length > 0) ? length : 1;
+	}
+    }
+  if (length > 0)
+    {
+      length--;
+    }
+
+  result = gsl_vector_calloc (length);
+  subs = xnmalloc (n_vars, sizeof (*subs));
+  for (j = 0; j < n_vars; j++)
+    {
+      if (var_is_alpha (v[j]))
+	{
+	  subs[j] = cat_value_find (v[j], vals[j]);
+	}
+    }
+  j = subs[0];
+  for (i = 1; i < n_vars; i++)
+    {
+      j = j * cat_get_n_categories (v[i]) + subs[i];
+    }
+  gsl_vector_set (result, j, 1.0);
+  /*
+     If any of the variables are numeric, the interaction of that
+     variable with another is just a scalar product.
+   */
+  for (i = 1; i < n_vars; i++)
+    {
+      if (var_is_numeric (v[i]))
+	{
+	  tmp *= vals[i]->f;
+	}
+    }
+  if (fabs (tmp - 1.0) > GSL_DBL_EPSILON)
+    {
+      gsl_vector_set (result, j, tmp);
+    }
+  free (subs);
+
+  return result;
+}
diff --git a/src/math/interaction.h b/src/math/interaction.h
new file mode 100644
index 00000000..4e709818
--- /dev/null
+++ b/src/math/interaction.h
@@ -0,0 +1,21 @@
+/* PSPP - Creates data structures to store interactions for
+   statistical routines.  
+
+   Copyright (C) 2007 Free Software Foundation, Inc.
+
+   This program is free software; you can redistribute it and/or
+   modify it under the terms of the GNU General Public License as
+   published by the Free Software Foundation; either version 2 of the
+   License, or (at your option) any later version.
+
+   This program is distributed in the hope that it will be useful, but
+   WITHOUT ANY WARRANTY; without even the implied warranty of
+   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+   General Public License for more details.
+
+   You should have received a copy of the GNU General Public License
+   along with this program; if not, write to the Free Software
+   Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
+   02110-1301, USA. */
+gsl_vector *
+get_interaction (union value **, const struct variable **, size_t);