--- a/src/modules/glm/samplers/AMMethod.h
+++ b/src/modules/glm/samplers/AMMethod.h
@@ -7,6 +7,19 @@
 
     class AuxMix;
 
+    /**
+     * @short Auxiliary Mixture sampler for Poisson and Binomial GLMs
+     *
+     * Sampler for Poisson and logistic regression proposed by Sylvia
+     * Fruhwirth-Schnatter, Rudolf Fruhwirth, Leonhard Held, Havard
+     * Rue (2009) Improved auxiliary mixture sampling for hierarchical
+     * models of non-Gaussian data. Statistics and Computing 19: 479-492
+     *
+     * The auxiliary mixture method uses data augmentation to replace
+     * the discrete outcome with a latent continous auxiliary
+     * variable, and then uses a finite mixture approximation to treat
+     * the auxiliary variables as normal.
+     */
     class AMMethod : public GLMMethod {
 	std::vector<AuxMix*> _aux;
       public:
@@ -14,9 +27,25 @@
 		 std::vector<GraphView const *> const &sub_views,
 		 unsigned int chain);
 	~AMMethod();
+	/**
+	 * Returns "AuxiliaryMixture"
+	 */
 	std::string name() const;
+	/**
+	 * Returns the precision parameter of the auxiliary variable
+	 * for the outcome variable with index i, according to the
+	 * current normal mixture approximation.
+	 */
 	double getPrecision(unsigned int i) const;
+	/**
+	 * Returns the value of the auxiliary variable for the outcome
+	 * with index i, according to the current normal mixture
+	 * approximation.
+	 */
 	double getValue(unsigned int i) const;
+	/**
+	 * Updates the auxiliary variables and then calls updateLM
+	 */
 	void update(RNG *rng);
     };