Revealing components of the galaxy population through nonparametric techniques

Abstract: The distributions of galaxy properties vary with environment, and are often multimodal, suggesting that the galaxy population may be a combination of multiple components. The behaviour of these components versus environment holds details about the processes of galaxy development. To release this information we apply a novel, nonparametric statistical technique, identifying four components present in the distribution of galaxy H$\alpha$ emission-line equivalent-widths. We interpret these components as passive, star-forming, and two varieties of active galactic nuclei. Independent of this interpretation, the properties of each component are remarkably constant as a function of environment. Only their relative proportions display substantial variation. The galaxy population thus appears to comprise distinct components which are individually independent of environment, with galaxies rapidly transitioning between components as they move into denser environments.

\section{Components of the galaxy population}
It has long been recognised that galaxies may be divided into at
least two distinct sub-populations. Originally this division was based
on visual appearance.  Most galaxies can be morphologically
classified as either elliptical or spiral.  Finer classification is
possible, discretizing an apparently continuous variation in galaxy
appearance.  However the dichotomy between elliptical and spiral
morphology is more pronounced than the variations within each
class.  Subsequently, it has been discovered that several other, more
quantitative, galaxy properties are distributed unevenly or in a
multi-modal manner.

The colour distribution of SDSS galaxies is strongly bimodal
\citep{2001AJ....122.1861S}. Galaxies in the ``red'' and ``blue'' modes
can be roughly identified as those with elliptical and spiral
morphology, respectively
\citep{2002AJ....124..646H,2006MNRAS.368..414D}.  Whereas morphology
reflects the dynamical state of galaxies, colour is related to their
star-formation history, particularly over the last $\la 10^9$
years.  The colour bimodality thus implies a division of the galaxy
population into blue galaxies, which have recently formed stars, and
red galaxies, which have not.  Such a bimodality in the star-formation
properties of galaxies has also been observed using more direct
measures of current star-formation, such as emission-line strength
\citep{2004MNRAS.348.1355B}.

The position of the red and blue galaxy sequences in the
colour--luminosity or colour--stellar mass planes display only a weak
dependence on environment.  However, the relative proportions of
galaxies in the two sequences vary strongly. In regions with a higher
local galaxy density the fraction of galaxies on the red sequence is
higher
\citep{2004ApJ...615L.101B,2004AIPC..743..106B,2006MNRAS.373..469B}.

It remains a matter of debate whether colour is more closely related to
environment than morphology.  Some claim that trends in morphology
versus environment can be mostly explained via a morphology--colour
relation which is almost independent of environment
\citep{2006MNRAS.366....2W,2006astro.ph..8353B,2006astro.ph.10171B,2007MNRAS.376L...1W}.
However, other studies oppose this view \citep{2007ApJ...658..898P},
and it has been clearly shown that the colour and morphology
bimodalities behave differently with respect to environment and
stellar mass \citep{2008arXiv0805.2612B}.

There are growing indications that, in a fraction of the galaxy
population, star-formation must be terminated rapidly
\citep{2004MNRAS.348.1355B,2006MNRAS.373..469B}.  The emission lines
in galaxy spectra provide a way of measuring the level of current star
formation on a timescale of $\la 10^7$ years.  They therefore trace
rapid star formation variations more sensitively than colour.  Another
important property of emission lines is that they are produced by
active galactic nuclei (AGN), in addition to star formation.  AGN are
present in many galaxies, and are thought to be produced by accretion
of material onto the super-massive black holes which appear to reside
at the centre of most, if not all, galaxies
\citep{1998Natur.395A..14R}.  Recently a variety of studies have
suggested that AGN strongly influence star-formation in their host
galaxies, and thus play an important role in defining the galaxy
population
\citep{2003MNRAS.346.1055K,2005MNRAS.364.1337S,2006MNRAS.365...11C,2006MNRAS.370..645B}.
The potential presence of an AGN contribution complicates the
traditional usage of emission-lines as an indicator of star
formation rate (SFR).  However, it also presents an opportunity to
study these two interdependent processes, star-formation and AGN,
through the distribution of a single quantity.

Galaxies with contrasting properties are found to be distributed
differently in space.  Elliptical galaxies cluster together more
strongly than spirals \citep{2000ApJ...545....6B,2001ApJ...554..857G}.
Similarly, red galaxies are preferentially found in denser
environments than blue galaxies \citep{2005ApJ...630....1Z}.  We have a
well developed theory for how structure forms in the cosmos, at least
in terms of the underlying cold dark matter which dominates the mass
density \citep{2005Natur.435..629S}.  Baryonic matter is expected to be
similarly distributed, in broad terms.  This theory thus explains the
range of galaxy environments observed.  However, the properties of
galaxies as a function of environment is a much more complicated
issue, depending on the detailed physics of galaxy formation and
evolution.  By studying trends in the galaxy population with
environment we can learn about these physical processes.

There has been a logical progression in studies of galaxy properties
as a function of environment.  Early work was based on dividing
galaxies into simple classes and looking at variations in the
fractions of galaxies of each class in bins of environment
\citep{1985ApJ...288..481D}.  As galaxy samples grew, this moved on to
examining trends in the mean properties of galaxies as a smooth
function of local galaxy density
\citep{2002MNRAS.334..673L,2003ApJ...584..210G}. A significant
development was fitting to the data functions that describe the
distribution of galaxies in two classes
\citep{2004ApJ...615L.101B,2004AIPC..743..106B,2006MNRAS.373..469B}.
Most of the approaches employed so far have relied upon enforcing a
predefined view of how to divide or classify the galaxy population in
increasingly complex ways.  However, our understanding of the physical
processes at work is highly uncertain and does not provide a
sufficient basis to make this decision.  Our only guide is the data
itself.  A natural next step is thus to turn to nonparametric methods,
where the components of the population are deduced consistently from
the data itself.

Recently, several studies have performed multivariate statistical
analyses on datasets containing a wide variety of galaxy properties,
in order to identify components of the galaxy population, and
determine which properties are most important for identifying to which
component a galaxy belongs
\citep{2005MNRAS.363.1257E,2006MNRAS.373.1389C}.  Such studies are
highly informative, but become complicated when one wishes to
determine the behaviour of the identified components versus another
variable.  In this paper we are primarily concerned with variation in
the components of the galaxy population as a function of environment.
The statistical method we present below may be straightforwardly
applied to multivariate datasets.  However, for simplicity, in the
present work we consider the environmental dependence of just one
galaxy property.  Nevertheless, even with this elementary approach, we
are able to learn much about the galaxy population.

\begin{figure*}
\includegraphics[clip=True,trim = 3cm 22.3cm 9.5cm 3.3cm,width=1.0\textwidth]{fig1.ps}
\caption{\label{fig:w13}
  The distribution of transformed H$\alpha$ equivalent width
  ($\W$) for (left) low and (right) high density environments.  The
  histogram displays the data, with Poisson uncertainties indicated by
  the grey shading.  The red, purple, green and blue lines show the
  components derived by applying the NMR technique.  The brown line
  gives the sum of these components, which is clearly a good
  representation of the data.}
\end{figure*}

\section{Conditional density estimation}
A common problem in astronomy, and statistical sciences in general, is
that one wishes to understand how the behaviour of one variable depends
upon another.  This is relatively straightforward in the case where
there is a single relationship between the variables, albeit with
some, possibly variable, scatter or width to the distribution.  Much
statistical and astronomical literature has been devoted to the
development of such regression methods \citep{weisberg}.  However, in
the case where multiple components may be present in the overall
distribution, each with a different functional dependence on the
variables, the situation becomes substantially more difficult.  One
can still attempt to apply single-component statistical tools, for
example nonparametric quantile regression \citep{QR,lQR}, on the whole
distribution, but the understanding one gains from such an exercise is
limited and sometimes misleading.  Alternatively one may individually
analyse subsamples selected by defining regions in the parameter
space, or preferably using additional information\citep{MMR}.  This
approach, however, is unsuitable when the multiple components
significantly overlap, or when it is unclear how many components are
present.

Most regression techniques focus on estimating the conditional mean,
the average value of one variable as a function of another variable;
for example, a line through a set of scattered points.  However, one
may get a better understanding of the relationship between a response
variable and a set of covariates by considering the estimation of the
conditional density as a whole; the \emph{distribution} of one
variable as a function of another.  (Note that \emph{density} here
refers to probability density as a function of the parameter set, not
a measure of environmental local galaxy density as elsewhere in this
paper.)  We use a new conditional density estimator based on finite
mixture models and local likelihood estimation, which describes the
underlying relationship between two variables by a set of
parameterised functions. This feature gives the proposed procedure the
advantage of being easily interpretable. This method is called
nonparametric mixture regression (NMR), and is described in detail in
Appendix \ref{sec:nmr}.

The NMR technique has the potential to aid the understanding of many
datasets, across all fields of science.  In the present work, it
allows us to determine the environmental dependence for individual
components of the galaxy population, with minimal prior assumptions on
the number and properties of these components.

\section{\boldmath Galaxy H$\alpha$ equivalent widths}
\label{sec:Halpha}
The strongest emission line in a galaxy optical spectrum is H$\alpha$.
The luminosity of H$\alpha$ is approximately proportional to the rate
of ongoing star-formation \citep{2006ApJ...642..775M}, when
uncontaminated by additional emission, such as from an AGN.  A
commonly employed quantity is the equivalent width (EW) of a spectral
line, the line flux normalised by the continuum flux at the same
wavelength.  The EW measurement has the advantages of being
approximately independent of uncertainties in the spectral flux
calibration and any extinction present in both the observed galaxy and
our own.  The H$\alpha$ line is in the red region of the spectrum,
where the continuum is dominated by the light from old stars.  The
H$\alpha$ continuum flux is therefore roughly proportional to stellar
mass, and hence $\EW$ is approximately proportional to the SFR
per unit stellar mass.

\begin{figure*}
\includegraphics[clip=True,trim = 3cm 22.3cm 9.5cm 3.3cm,width=1.0\textwidth]{fig2.ps}
\caption{\label{fig:ew13}As \reffig{fig:w13}, but shown here in terms of the
  untransformed equivalent width, $\EW$.  The inset shows the same
  plot with axis-ranges chosen to better show the behaviour at small $\EW$.}
\end{figure*}

The overall distribution of galaxy H$\alpha$ luminosity, equivalent
width, and hence absolute and normalised SFR, are
found to move to lower levels with increasing environmental density
\citep{2002MNRAS.334..673L,2003ApJ...584..210G}.  This generally agrees
with the colour and morphology trends described above, and the
variation of H$\alpha$ emission with morphological type
\citep{2004AJ....127.2511N}.  However, if the galaxy population is
separated into galaxies which are star-forming and those which are
not, the distribution of $\EW$ for each component
does not depend significantly on environment.  Only the relative
proportion of star-forming galaxies changes strongly
\citep{2004MNRAS.348.1355B} with environment.  This finding, of
distinguishable components in the galaxy population with properties
independent of environment but proportions which vary strongly,
mirrors the behaviour found in the colour distribution.  It also
motivates us to perform a more rigorous evaluation of the components
present in the galaxy population in this work.

As mentioned earlier, an important feature of emission lines is that,
in addition to star formation, they are also produced by AGN.
Galaxies whose emission lines are dominated by star-formation or AGN
activity can be separated using various diagnostic diagrams.  The most
common of these plots the emission line ratios
${\rm{[OIII]}\lambda5007}/{\rm{H}\beta}$ versus
${\rm{[NII]}\lambda6583}/{\rm{H}\alpha}$, and is known as the BPT
diagram \citep{1981PASP...93....5B}. The usual approach is to use these
diagrams to reject objects inappropriate to the particular study.
Thus a study of galaxy star formation properties would exclude all
galaxies with signs of AGN contamination.  However, classifying a
galaxy using the BPT diagram requires multiple emission lines to be
detected, resulting in a fraction of objects which cannot be
classified.  In addition, the separation between galaxies dominated by
star-formation and AGN is not clear, and there appears to be a large
population of galaxies which host both star-formation and an AGN.
Roughly 20\% of all galaxies are unambiguously AGN-dominated, while it
is estimated that a further 20\% are star-forming galaxies with a
significant AGN contribution \citep{2003ApJ...597..142M}. This
ambiguity means a variety of SFR--AGN demarcations exist
\citep{2001ApJ...556..121K,2003MNRAS.346.1055K,2006MNRAS.371..972S}.
Star formation studies based on emission lines have therefore rejected
widely varying fractions of galaxies from their samples.  This
fraction is usually low, so significant numbers of AGN-contaminated
galaxies remain.  More importantly, if our aim is to gain knowledge of
star-formation properties across the whole galaxy population, then we
may be rejecting an important fraction of the population.  If there
are any intrinsic correlations between AGN and star-formation, as has
been suggested by other studies \citep{2003MNRAS.346.1055K}, then
information about these will be lost.

A number of classes of AGN have been identified.  A primary
distinction is between Type 1 and Type 2 AGN.  In Type 1 objects our
viewing angle is such that we see the region immediately around the
central black hole directly, and thus the galaxy's light is dominated
by the AGN emission.  In this case the properties of the host galaxy
are generally very difficult to determine.  In Type 2 AGN, the central
region is obscured by a dusty torus surrounding it.  The observed AGN
emission is therefore due to material further removed from the central
ionising source, and mostly confined to emission lines.  Most
photometric and structural galaxy properties may therefore be reliably
measured, despite the presence of a Type 2 AGN.  In this work we
exclude all Type 1 AGN, identified by the large widths of their
emission lines, and consider only the more common Type 2 objects.  A
further subdivision within Type 2 AGN is between LINER and Seyfert 2
objects.  These are similar, and may simply be two parts of a
continuum of objects, with Seyfert 2 AGN being more powerful and
highly ionised.  However, there are signs that LINERs and Seyfert 2
AGN are truly physically distinct classes \citep{2006MNRAS.372..961K}.

In this work we examine the components in the distribution of galaxy
$\EW$, interpretable as a proxy for star formation rate and
nuclear activity per unit stellar mass.  It is possible to estimate
the true star-formation rate and stellar mass, for galaxies which do
not host an AGN, using a combination of several spectral features.
However, such estimates are sensitive to the details of the assumed
model.  There is therefore a concern that any finding concerning the
components of the resulting distribution may be attributable to the
model.  The $\EW$, on the other hand, is a single, robust,
model-independent measurement.

The data we use in our study is from Data Release 4 of the SDSS
\citep{2006ApJS..162...38A}. The emission line fluxes, continua and
resulting EW used in this study are those provided for DR4 by the
MPA-Garching group \citep{2004ApJ...613..898T}\footnote{available from
  http://www.mpa-garching.mpg.de/SDSS/DR4}.  All quantities used in
this paper were obtained from the CMU-PITT SDSS DR4 Value Added
Catalog\footnote{available from
  http://nvogre.phyast.pitt.edu/dr4\_value\_added} (VAC).  The SQL code for
the selection of each of our samples is given in \reftab{tab:sql}.
We construct a volume-limited sample by selecting galaxies with $0.05
< z < 0.095$ and $M_r < -20.4$.  In this work we thus focus on the
behaviour of fairly bright galaxies.  The lower redshift limit ensures
the spectra are based on a reasonable fraction of the galaxies' light;
at $z=0.05$ the $3$~arcsec diameter of each spectroscopic fibre
corresponds to $3$~kpc.  Throughout we convert to physical scales
assuming a flat Friedman-Robertson-Walker cosmology with $\Omega_m =
0.3$, $\Omega_{\lambda} = 0.7$ and $H_0 = 70$ km~s$^{-1}$~Mpc$^{-1}$.

\begin{table}
  \caption{\label{tab:sql}Definitions of the galaxy samples used in this
    study, given as `where' clauses of the SQL queries of the CMU-PITT SDSS
    DR4 VAC}
\begin{tabular}{p{0.06\textwidth}p{0.31\textwidth}c}
\hline
\centering sample & \centering SQL selection & $n$ \\
\hline \hline
density defining sample &
\texttt{!z between 0.02 and 0.10 and absolute\_Petro\_r <= -20.4 and Sort~=~0}
&
117873\\
\hline
$\rho_{1.3}$ \mbox{sample} &
\texttt{!z between 0.05 and 0.095 and absolute\_Petro\_r <= -20.4 and
 2.4 < Dist\_right\_edge and 2.4 < Dist\_left\_edge and
 2.4 < Dist\_upper\_edge and 2.4 < Dist\_lower\_edge and
 H\_ALPHA\_FLUX > -99 and H\_ALPHA\_CONT > 0.0001 and
 H\_ALPHA\_FLUX/H\_ALPHA\_CONT > -0.4 and
 absolute\_Petro\_u > -990 and absolute\_Petro\_r > -990
 and Sort~=~0}
&
76420\\
\hline
$\rho_{5.5}$ \mbox{sample} &
\texttt{!z between 0.05 and 0.095 and absolute\_Petro\_r <= -20.4 and
 11 < Dist\_right\_edge and 11 < Dist\_left\_edge and
 11 < Dist\_upper\_edge and 11 < Dist\_lower\_edge and
 H\_ALPHA\_FLUX > -99 and H\_ALPHA\_CONT > 0.0001 and
 H\_ALPHA\_FLUX/H\_ALPHA\_CONT > -0.4 and
 absolute\_Petro\_u > -990 and absolute\_Petro\_r > -990
 and Sort~=~0}
&
46998\\
\hline
\end{tabular}
\end{table}

\section{\boldmath Measuring galaxy environment}
Galaxy environment can be characterised in many ways, but a commonly
adopted value is the local number density of galaxies brighter than a
given luminosity, averaged over some volume or kernel.  We estimate
the local galaxy number density, $\rho_b$, within a fixed-scale,
spherical kernel with a Gaussian radial profile and bandwidth $b$.
Our local galaxy densities are thus simple to interpret physically.

To select the bandwidth, or scale, $b$ of the kernel, we apply
leave-one-out cross-validation; that is, we select the value of $b$
which minimizes the estimated integrated mean squared error, $CV(b)$. This
error is obtained by estimating the density function $n$ times, each
time leaving out one galaxy from the estimation:
\begin{equation}
CV(b) = \int \widehat f_{n,b}^{\;2}(\bmath{x}) d\bmath{x} -
\frac{2}{n}\sum_{i=1}^n \widehat f_{(-i),b}(\bmath{X_i})
\end{equation}
where $\{\bmath{X_i}\}$ is the set of galaxy positions, and $\widehat
f_{n,b}$ and $\widehat f_{(-i),b}$ are the kernel density estimators
with bandwidth $b$, using all $n$ galaxies and after removing the
$i^{\rmn{th}}$ galaxy, respectively.  We compute $CV(b)$ for a range
of different bandwidth values to find that which minimizes the error.
Applying this cross-validation method we determine an optimum
bandwidth value of $1.3$~Mpc.  A similar optimum bandwidth for local
galaxy density estimation was found using cross-validation by
\citet{2004MNRAS.348.1355B}.

Interestingly, this scale corresponds to the size of
galaxy clusters, and is thus highly appropriate for characterising
density from a physical, as well as a statistical, point of view.
However, while cross-validation provides the statistically optimum
bandwidth for the whole sample, any choice of bandwidth has its
limitations.  This density estimator loses resolution at low
densities, where there are no neighbouring galaxies within the kernel
bandwidth, and is thus unable to discriminate between densities lower
than $\rho_{1.3} \sim 0.03$~Mpc$^{-3}$, comprising 17\% of the sample.
In order to probe environments less dense than this, but necessarily
on larger physical scales, we additionally perform the analysis with
local densities measured using a larger bandwidth of $5.5$~Mpc.
Almost all galaxies have a neighbour within this radius.  One could
also consider estimating densities with a kernel bandwidth
significantly smaller than $1.3$~Mpc. However, such an estimator would
lose resolution below even moderate densities, where galaxies are
typically separated by more than the bandwidth.  It would also be less
able to discriminate between high density environments, because the
densities are estimated using galaxy positions uncorrected for
redshift-space distortions, and hence an increase in true-space
density no longer results in a higher redshift-space density within
the kernel.  We mostly show results based on the
statistically-motivated $1.3$~Mpc bandwidth in the main body of this
article, but provide figures using the $5.5$~Mpc bandwidth in Appendix
\ref{sec:55Mpc}, to demonstrate that we find similar results on larger
scales and to lower densities.

We avoid biased density estimates for galaxies at the edges of our
sample volume by determining the densities using a larger volume
sample of galaxies with $0.02 < z < 0.10$ and $M_r < -20.4$.
We then limit the analysis sample to galaxies with $0.05 < z < 0.095$
and further than approximately twice the bandwidth from a survey
boundary.  We reject a further 3\% of galaxies with unreliable $\EW$
or $(u-r)$ rest-frame colour measurements.  The exact selections, and
corresponding sample sizes, are given in \reftab{tab:sql}.

\begin{figure*}
\includegraphics[clip=True,trim = 3cm 22cm 5cm 3.3cm,width=1.0\textwidth]{fig3.ps}
\caption{\label{fig:wrho13}The behaviour of the NMR components versus
  environment.  The left panel plots the data as dots, along with the
  location of each component, indicated by thick, solid lines, and additionally
  their widths via the coloured shading and dashed lines.  These widths
  are shown explicitly in the middle panel.  The right panel displays
  the variation in the proportion of each component.  While the
  location and width of the components do not change significantly
  with environment, the proportions vary strongly.}
\end{figure*}

\section{Applying the NMR technique}

\begin{figure*}
\includegraphics[clip=True,trim = 5cm 7.5cm 10cm 2cm,width=1.0\textwidth]{fig4.ps}
\caption{\label{fig:3Dwrho55}A three-dimensional view of the NMR estimate of the
  $\W$--$\rho_{5.5}$ distribution, shown by the grey, transparent
  surface, and its constituent components, colour-coded is in the
  previous figures.  It can be clearly seen that the positions and widths of the
  components do not change significantly, while their relative
  proportions vary substantially.}
\end{figure*}

A brief inspection of the sample $\EW$ distribution reveals a peak
around zero EW, with a long, asymmetric tail to high EW.  The NMR
technique is more computationally efficient when using symmetrical,
Gaussian functions to model the distribution.  Gaussians are also an
obvious choice due to their exceptional richness and flexibility.  For
convenience we therefore wish to transform the equivalent width
quantity to a space where its natural components appear to take a more
symmetrical, Gaussian, form.  Better matching the shape of the true
distribution components to that assumed in the NMR technique will also
naturally result in fewer NMR components being required to model the
distribution (but see Appendix \ref{sec:nmr}).  The EW extend slightly to
negative values, proscribing a simple logarithmic transformation.  We
therefore choose the transformation $\W = \log_{10}(\EW + \lambda)$.
The zero offset parameter, $\lambda$, must be large enough to make the
logarithm argument positive for the most negative EW value in our
sample. In constructing our sample we remove outliers by requiring
$\EW > -0.4$, thereby clipping the lowest 0.1\% of the sample.
Therefore, we must have $\lambda > 0.4$.  We have examined the
behaviour of our NMR fits and their likelihood with variations in
$\lambda$.  The chosen value has only a relatively small effect,
slightly altering the shape of the Gaussian basis functions once they
are transformed back into EW space, but not changing our results
significantly.  Here we adopt $\lambda = 1.4$ as a compromise between
maximising the fit likelihood and ensuring stable behaviour.  We must
also choose a reasonable bandwidth for the regression kernel in
$\rho$.  Following extensive tests we adopt an adaptive bandwidth
enclosing the nearest 5000 points (also see discussion in Appendix
\ref{sec:nmr}).

We apply the NMR technique to the distribution of $\W$, and determine
the optimum number of components using the Bayesian Information
Criterion \citep[BIC;][]{BIC}.  Four components are strongly preferred
by the data, by $\Delta$BIC $>$ 7 (see Appendix \ref{sec:nmr} for more
details).  In \reffig{fig:w13} (\reffig{fig:w55}) we show the NMR
components we obtain for the $\rho_{1.3}$ ($\rho_{5.5}$) sample, at
two values of local galaxy density.  The components are plotted in
$\W$-space, in which the technique is applied.  We also show the
components and data transformed back into $\EW$-space in
\reffig{fig:ew13} (\reffig{fig:ew55}).  The properties of these
components as a function of environmental density are shown in
\reffig{fig:wrho13} (\reffig{fig:wrho55}).  In \reffig{fig:3Dwrho55}
we show a three-dimensional view of the components and their sum for
the $\rho_{5.5}$ sample, which includes all the relevant information
(location, width and relative proportion of each component) in a
single plot.  We show the results for the $\rho_{5.5}$ here simply
because they are smoother than those for $\rho_{1.3}$, and the
individual components are more clearly visible in this
three-dimensional view.  It is critical to note that the only data
which has been used to determine these components is the $\EW$
distribution.

At this stage we make no attempt at interpreting the components as
physically distinct populations. Nevertheless, Figs. \ref{fig:w13},
\ref{fig:w55}, \ref{fig:ew13}, \ref{fig:ew55} indicate that the $\EW$
distribution can be well described by multiple components.  The
hypothesis that the galaxy population comprises distinct components,
or types, is strongly supported by the various property bimodalities
described earlier.  We find that the locations and widths of the
components of the $\EW$ distribution are independent of environment.
Only the relative proportions of the components are found to vary
strongly.  This implies that the variations with environment are
primarily the result of differences in the relative frequency of each
galaxy type, rather than changes in the intrinsic properties of each
type.

Galaxies move to regions of higher density over time, under the
influence of gravity.  The variation of galaxy properties with
environment is therefore at least partly due to
environmentally-dependent changes in individual galaxy properties over
time.  If all galaxies in a given environment were affected similarly,
we would expect to see smooth changes in the property distributions of
each individual component.  However, we find that the individual
components remain mostly unchanged with environment.  This implies
that some galaxies are transformed directly from one type to another,
in an apparently stochastic manner.  If this transformation is
sufficiently slow, we would expect to see the transitioning galaxies
appearing as a separate component in the relevant range of local
density.  If it is rapid, then the fraction of transitioning galaxies
at any time would be too low to separate from the main distribution.

\begin{figure}
\includegraphics[angle=270,width=0.45\textwidth]{fig5.ps}
\caption{\label{fig:bpt13}The BPT diagram for our $\rho_{1.3}$ sample, traditionally
  used to identify star-forming galaxies and AGN hosts.  For clarity,
  only one-fifth of our sample galaxies are plotted. The \emph{LINER},
  \emph{Seyfert 2} and \emph{SF dominated} regions are colour-coded to
  match our interpretation of their correspondence to the NMR
  components shown in the other figures (purple, blue and green,
  respectively).  Note that many galaxies cannot be placed on this
  diagram.  These are \emph{passive} galaxies, with no emission lines,
  and \emph{uncertain} galaxies, with some detected emission lines,
  but not all four of those required for inclusion in this diagram.}
\end{figure}

\section{Identifying the components}
It is easy to identify the component at zero $\EW$ with passive
galaxies, containing no star-formation or AGN activity.  The dominant
component at high $\EW$ must be associated with star-forming galaxies
(with the above caveats concerning potential AGN contamination).  We
also find two intermediate $\EW$ components.  The principle change
with environment appears to be the movement of galaxies from the
star-forming component to the others, but primarily to the passive
component.  However, interpreting either of these intermediate EW
components as a population transitioning between star-forming and
passive is inconsistent with their existence as a significant fraction
of the galaxy population even at low environmental densities.

To explore the physical interpretation of the components we have
found, we now turn to more traditional diagnostics to separate the
contributions from star formation (SF) and AGN to the emission lines.
The BPT diagram for our $\rho_{1.3}$ sample is shown in \reffig{fig:bpt13}.  In
order to appear on this plot, all four required emission lines must be
detected at $>2$~sigma significance.  The classifications we define
are as follows;
\emph{passive}: no emission lines detected,
\emph{SF dominated}: all four lines detected and below the curve of
\citet{2006MNRAS.371..972S},
\emph{AGN dominated}: above the line of \citet{2001ApJ...556..121K}
with either all four lines detected or with both lines for just one of
the ratios detected and ${\rm{[OIII]}}/{\rm{H}\beta} > 0.6$ or
${\rm{[NII]}}/{\rm{H}\alpha} > 0.05$,
\emph{AGN+SF}: all four lines detected and between the curves of
\citet{2001ApJ...556..121K} and \citet{2003MNRAS.346.1055K},
\emph{SF+AGN}: all four lines detected and between the curves of
\citet{2006MNRAS.371..972S} and \citet{2003MNRAS.346.1055K},
\emph{uncertain}: at least one of the four emission lines detected,
but none of the other classification criteria met.
Note that the majority of AGN-dominated galaxies can be robustly
identified simply from their ${\rm{[NII]}}/{\rm{H}\alpha}$ ratio
\citep{2003ApJ...597..142M,2006MNRAS.371..972S}.

Our classification method is such that galaxies classified as
\emph{AGN dominated} must contain a significant AGN component, and
will have low contribution to their emission lines from star
formation.  On the other hand \emph{SF dominated} galaxies may well
also contain up to $\sim 20$--$40$\% AGN contamination in their
emission lines \citep{2003MNRAS.346.1055K,2006MNRAS.371..972S}.  The
\emph{AGN dominated} galaxies can be further subdivided into
\emph{LINER} and \emph{Seyfert 2} sources using the BPT diagram
\citep{2003MNRAS.346.1055K}.

\begin{figure*}
\includegraphics[clip=True,trim = 3cm 22cm 9.5cm 3.3cm,width=1.0\textwidth]{fig6.ps}
\caption{\label{fig:wrho13bpt}The $\W$--$\rho_{1.3}$ distribution for
  objects in our sample colour-coded by their location in the BPT
  diagram shown in Fig.~4.  The lines indicate the median $\W$ in bins
  of $\rho_{1.3}$ for each subsample.  The left panel shows
  \emph{passive}, \emph{LINER}, \emph{Seyfert 2} and \emph{SF
    dominated} galaxies (in order of increasing $\W$), while the right
  panel shows \emph{uncertain}, \emph{AGN+SF} and \emph{SF+AGN}
  galaxies (brown, orange and cyan, respectively, and again in order
  of increasing $\W$).  A comparison with Fig.~2 reveals a
  correspondence between the NMR components and, in order of
  increasing $\W$, (1) \emph{passive} galaxies, (2) \emph{LINER} and
  \emph{uncertain} galaxies, (3) \emph{Seyfert 2} and \emph{AGN+SF}
  galaxies, and (4) \emph{SF dominated} and \emph{SF+AGN} galaxies.}
\end{figure*}

Figure \ref{fig:wrho13bpt} shows the $\W$--$\rho_{1.3}$ distributions
of galaxies classified using the BPT diagram. Comparing with
\reffig{fig:wrho13}, one can clearly identify the NMR components with
the \emph{passive}, \emph{LINER}, \emph{Seyfert 2} and \emph{SF
  dominated} BPT-classified galaxies.  The large fraction of galaxies
for which the BPT diagram gives an uncertain result may also be
identified with the components.  The galaxies with apparently mixed
star formation and AGN emission are found at similar $\W$ to the
\emph{Seyfert 2} objects, and the higher intermediate NMR component.
Galaxies with at least one emission line, but which cannot be
identified via the BPT diagram have similar $\W$ to \emph{LINER}
objects and the lower NMR component.  While not conclusive, this
strongly suggests that the components derived from the NMR technique
do represent physically distinct populations.  This is remarkable
given that the NMR components have been inferred from just a single
emission line.

\section{A new insight into the galaxy population}

By applying the newly developed NMR method to the H$\alpha$ equivalent
width distribution, a single astrophysical quantity that contains
information on both star formation and nuclear activity, we have
identified four distinct components in the galaxy population.  None of
these components vary significantly with environment, in terms of the
distribution of their H$\alpha$ equivalent widths.  However, the relative
proportions of galaxies in each component vary substantially with
environment.  This implies that any environmental processes at work do
not affect all galaxies in a gradual way, which would result in
changes in the component H$\alpha$ equivalent width distributions.
Rather, they must rapidly transform a fraction of galaxies from one
component to another, in a stochastic manner, in order to avoid
changing the properties of the individual components.

The above conclusions stand without requiring us to identify the
components with more traditional galaxy sub-populations.  However,
when we attempt such an identification, we find that the extreme
components may be associated with passive and star-forming galaxies,
while the two intermediate components display similarities to galaxies
hosting LINERs and Seyfert 2 AGN.  Galaxies with an apparent mix of
star-formation and AGN may also be identified with these components.
However, in contrast to the usual methods of classifying the
star-formation and AGN properties of galaxies, which require multiple
emission lines to be significantly detected, the technique we describe
in this paper is applicable to all galaxies.  We thereby avoid the
issue of excluding objects for which traditional methods are
uncertain, and the biases which this may introduce.

\section*{Acknowledgements}
SPB acknowledges support from an STFC postdoctoral grant.  AR
acknowledges the Qatar Foundation for Education, Science and Community
Development.  RCN holds a Marie Curie Excellence Chair from the
European Commission.  We thank the NSF for funding this
inter-disciplinary research through their KDI initiative.
Three-dimensional visualisation was conducted with the S2PLOT
programming library \citep{2006PASA...23...82B}.  We are grateful to
the referee, Dr. Nicholas Ball, for useful comments.

\bsp


Part I - A case study in post-secondary mathematics: the importance of mental health awareness

Abstract: Due to the shift of civilization from the Industrial Age to the Information Age, mathematical literacy has become a necessity in the twenty-first century. However, in order to learn and contribute to the mathematical community, one has to be in a state of good mental health. When I say 'good mental health,' I mean one has to develop a set of healthy coping strategies, be in a positive learning environment, and have a social support system. Traditionally, universities have been the venue of such higher learning. If these institutions want to remain as thriving grounds for higher education and engines of research, post-secondary institutions need to be aware of these factors and actively contribute to the well-being of its students and faulty. Unfortunately, students do not always receive the support necessary to be mentally healthy. The purpose of this paper is to examine how social awareness and sensitivity of mental health in a university setting is a key component for individuals to flourish academically and grow personally. Several voluntary surveys completed by my first and second year math students at the University of Waterloo will be presented. The surveys investigated the students' experiences at the university, particularly in the Mathematics Department. In addition, this paper explores an alternative way to structure math classes, specifically assignment scheduling, in order to help students develop healthy study habits. The scores for the alternative-scheduled assignments are statistically analyzed and compared to scores of the typical-scheduled assignments currently implemented at the University of Waterloo. Finally, this paper briefly discusses the importance of socialization, specifically for young mathematicians and scientists, and potential consequences of reduced social exposure in the Digital Age.

\section{Introduction}

According to the Mood Disorder Society of Canada (MDSC), \textit{mental health} and \textit{mental illness} is defined as follows:

\begin{itemize}
\item[] \textbf{Mental Health:}

\textit{(or well-being) is an ideal we all strive for. It is a balance of mental, emotional, physical and spiritual health. Caring relationships, a place to call home, a supportive community, and work and leisure all contribute to mental health. However, no one's life is perfect, so mental health is also about learning the coping skills to deal with life's ups and downs the best we can.}

\item[] \textbf{Mental Illness:	}

\textit{is a serious disturbance in thoughts, feelings and perceptions that is severe enough to affect day-to-day functioning} \cite{1}.

\end{itemize}
\noindent However, in this paper, I will instead be using the term \textit{poor mental health} rather than \textit{mental illness} to avoid the negative association with the term `illness.'

Some specific types of poor mental health that one typically hears about are for instance: schizophrenia, depression, bi-polar disorder, anxiety disorders, eating disorders, borderline personality disorder, autism-spectrum disorder, post-traumatic stress disorder (PTSD), etc. \cite{1}. Why such forms of poor mental health develop depends on many factors such as genetics and one's environment, but this is an over simplification of the issue \cite{2}. The issue is highly individual, and globally complex as a result. Researchers still know very little about poor mental health and the relation between one's environment and gene expression. It is a topic that is slowly starting to be taken more seriously, especially given recent media coverage which describes poor mental health as an increasing pandemic \cite{3}. Nonetheless, there is still significant and substantial stigma, and naive notions when it come to mental health.

The younger generations are especially vulnerable to poor mental health. According to the MDCS, 18$\%$ of Canadians alone between the ages of 15 - 24 are reported to have poor mental health or substance abuse problems. Sadly, suicide accounts for 24$\%$ of all deaths among Canadians aged 15 - 24, and 16$\%$ of all deaths for the age group 25 - 44. The people with the highest rate of depression are people under the age of 20; and people between the ages of 20 - 29 have the highest rate of anxiety \cite{1}. This is a major problem, especially if we want to have a functional, altruistic society in the future. This is the generation that we expect will be solving some, if not all, of the world's most difficult problems. Issues such as climate change, population growth and population sustainability, economic inequality (perhaps even the collapse of the monetary system), the impact of integration of automation technologies and robots into everyday society, and so on \cite{3}, \cite{4}, \cite{5}.

In order to solve these problems, people need to be scientifically literate and mathematically competent and proficient. In general, modern functioning society is highly dependent on technology; particularly on science, technology, engineering and mathematics (STEM) educated individuals. With the high probability of most repetitive jobs vanishing due to technological developments in automation \cite{4} -  STEM education and training will be essential and necessary for twenty-first century careers. To learn such material requires a lot of intense, focused and consistent work. To complete such work, stable and reliable mental health must exist. This implies an individual's mental health is just as important as physical health for such intense learning to take place. Hence, it is essential that we address poor mental health. This is particularly important for our young people who grew up with the Internet due to their increased physical isolation relative to previous generations. Unfortunately, very little research has been done on the impact of mental stability and the development of young adults' learning capacities and motivations.

As a young mathematics instructor at a post-secondary institution, the University of Waterloo in Canada, mental health became a major concern of mine, consuming substantial time.  The longer I teach these young people, the more I see their struggle in learning due to poor mental health \cite{6}. Unfortunately, most of the time the problem is not recognized or taken seriously. Observing this first hands motivated me to write this paper to help the well-being of my students. The focus will be to present supporting evidence that a shift in the primary and post-secondary education systems' culture is needed to be:

\begin{itemize}
\item more aware of mental health,
\item aware of how poor mental health affects the learning process, and
\item aware of how an institution's awareness affects their students.
\end{itemize}

\noindent\section{A commentary on surveys}

In order to gain more perspective and information on the mental health and opinions of my students, I created several surveys for them to complete if so desired. At the start of the Fall 2014 term at the University of Waterloo, I made the surveys available to the students I had taught in the previous terms, and the students I was teaching at the time. More specifically, I prepared these surveys using Google Forms for the mathematics students at the University of Waterloo. This included students in the undeclared math program (first and second year students), as well as students in the applied math, pure math, statistics, combinatorics and optimization, mathematical physics, and double degree math-business programs. The students who completed these surveys did so out of voluntary compliance; it was not part of any class and the students were informed the results would be summarized in a paper. Although I was skeptical that anyone would complete them, to my surprise, many students were willing and enthusiastic about completing (and helping create) the surveys. This demonstrated to me that the students wanted their opinions to be known; they just required an appropriate outlet.

For the survey results shown below, the total number of students who answered `yes' or `no' is given. The number of students who did not answer `yes' or `no' but rather left a comment, is not included in the mentioned numerical values. However, some of these comments can be found in Appendix B.

\noindent\section{Data From Surveys}

The purpose of the initial surveys was to determine the students' views on mental health and how the university setting influences their mental well-being. The following results were obtained:

\begin{itemize}
\item 82$\%$ of 238 students said: \textit{the professor I have for a class influences how much I work in that class.}

\item 88$\%$ of 236 students said: \textit{I am more motivated in class if I like the professor.}

\item 66$\%$ of 227 students said: \textit{a professor's awareness on mental health influences my enjoyment in a class.}

\item 95$\%$ of 234 students said: \textit{professor's should be aware of mental health.}

\item 84$\%$ of 221 students said: \textit{discussions on mental health are important.}

\item 77$\%$ of 215 students said: \textit{I believe I am aware of my own mental health.}

\item 77$\%$ of 212 students said: \textit{I believe that my peers and I are not aware of each other's mental health.}

\item 65$\%$ of 136 students said: \textit{we did not have discussions on mental health in high school.} 40 \% of the remaining students (35$\%$ of 136 students) who said: \textit{we did have discussion on mental health in high school}, also said: \textit{I do not think the school did a good job.}

\item 88$\%$ of 102 students said: \textit{professors should be able to take mental health into consideration.}

\item 61$\%$ of 102 students said: \textit{mental health has not been addressed thoroughly at university.}

\item 72$\%$ of 107 students said: \textit{I think professors should have training on mental health.}

\item 48$\%$ of 167 students said: \textit{I think that the professors at the University of Waterloo are not aware of mental health.}
\end{itemize}

The following questions were also asked on the student surveys:
\begin{itemize}
\item[] Q1: Are there any departments at the university that you feel are less aware of mental health?
\item[] Q2: Do you think poor mental health will become a larger issue in society as time progresses?
\item[] Q3: Has the university affected your mental health?
\item[] Q4: Do you think mental health is compromising your success in classes?
\end{itemize}

Fifty-two students responded to Q1, 58 responded to Q2, 83 answered Q3 and 102 answered Q4. The distribution of student responses to Q1 can be seen in Fig. 1. It can be seen that 52$\%$ of students said the Mathematics Department in particular was not very aware of the students' mental health. Engineering was the second highest ranked department at 17$\%$.

\begin{figure} [ht]
\centering
\includegraphics[scale=0.45]{MH_pie1.png}
\caption[]{Distribution of the 52 student responses to Q1.}
\label{fi:1}
\end{figure}

\begin{figure} [ht]
\centering
\includegraphics[scale=0.74]{MH_pie2.png}
\caption[]{Distribution of student responses to: Q2 top (58 responses), Q3 middle (83 responses), Q4 bottom (102 responses).}
\label{fi:1}
\end{figure}

The distribution of student responses to Q2, Q3 and Q4 can be seen in the top, middle and bottom plot of Fig. 2, respectively. It was found that 81$\%$ of students thought that poor mental health will become a larger issue in society as time naturally progresses. This is a hypothesis of mine, and my conjecture is that poor mental health is related to the larger amounts of information that everyone is exposed to on a regular basis. Consider that within a single minute, more information is uploaded onto the Internet than a person 200 year ago was exposed to in their entire lifetime \cite{7}. Such volumes of information will have some of influence on our minds and consciousness. Humanity has progressed this far due to natural selection and thus, how our physical bodies have evolved. Surely our minds and consciousness will evolve as our environment changes with respect to time. Therefore, I hypothesize that `mental illness' is an evolution of our minds responding to the volume of information exposed to us.

In the middle plot of Fig. 2, one can see that 44$\%$ of students surveyed think the university has influenced their mental health in a negative way; whereas 35$\%$ say the university has influenced their mental health in a positive way. Unfortunately, as depicted in the bottom plot of Fig. 2, 46$\%$ of the students surveyed think their success in the classroom is being compromised because of their mental health. Although the potential connection to gender (and the corresponding issues) will not be discussed in this paper (it will be the focus of a future paper), I hypothesize that gender equality (or a lack thereof) also has an effect on the mental health of women at the university. I speculate this is an issue because of my own experiences and from the stories women at the university have shared with me.

\noindent \section{Restructuring of Assignments}

\noindent \subsection{Introduction and Purpose}

I always tell my students that learning mathematics is like learning a language; and no one is born speaking automatically. A person needs a certain amount of exposure and practice before being able to speak any language. Therefore, it cannot be expected for individuals to form proper sentences before becoming familiar with words, sentence structure and grammar. It is natural and in fact, a part of the learning process to mispronounce words or use them incorrectly. This same argument also holds for the language of mathematics. People can find learning mathematics to be challenging and even overwhelming because of the unfamiliarity. Making mistakes is an unavoidable part of the learning process. However, this does not imply that mistakes (i.e., failing) should be treated as catastrophes unless one approaches failure in an unproductive way.

\begin{figure} [ht]
\centering
\includegraphics[scale=0.62]{learningCurve.png}
\caption[]{Ebbinghaus forgetting curve based on fading theory \cite{11}.}
\label{fi:1}
\end{figure}
\vspace{2mm}

Studies have consistently shown that review is essential in order to learn new information. The study of memory was initiated by the German psychologist named Hermann Ebbinghaus. In 1885, he published a paper describing what we now call the \textit{learning curve} and the \textit{curve of forgetting} \cite{12}, \cite{13}. Ebbinghaus' curve of forgetting is depicted in Fig. 3. It illustrates the exponential decay of learning new information from a single learning session \cite{13}. It indicates that if one does not \textit{actively} review newly learned information within the first 24 hours, only 33$\%$ of the information originally absorbed will be retained. Therefore, students should review their notes within the first 24 hours after attending class to help store the information to long-term memory. Fig. 4 depicts how reviewing on a regular basis can help ensure one does not forget new material \cite{14}. However, it is essential to review from day one. If students consistently prepare for class, learning will be a much more rewarding and enjoyable experience.

\begin{figure} [ht]
\centering
\includegraphics[scale=0.72]{reviewCurve.png}
\caption[]{Ebbinghaus forgetting curve and importance of reviewing \cite{14}.}
\label{fi:1}
\end{figure}

Due to these findings, I structured my classes for the last couple of years to have \textit{end of lecture} (EoL) assignments rather than the typical weekly assignments given in university math classes. However, this was a decision that I personally made for my own sections. I wanted to investigate if such a structure would be beneficial in a multi-section course, i.e., if all sections did EoL assignments rather than weekly assignments. In the Fall 2014 term at the University of Waterloo, I had the opportunity to carry out such an experiment. I coordinated the second year honours math linear algebra class, referred to as \textit{MATH 235}. This course is required for all math majors, but typically is not considered the most enjoyable course.

\noindent \subsection{Feedback from Students}

For the MATH 235 Fall 2014 class, all 459 students had a class every Monday, Wednesday and Friday. After each class, the students were given three to five questions that were due on the next day of the class. Upon surveying the students three weeks into the term, the following was discovered:

\begin{itemize}
\item 94$\%$ of 110 students said: \textit{I think the EoL assignments are making me keep up with the course material.}

\item 91$\%$ of 110 students said: \textit{I feel like the EoL assignments are improving my learning.}

\item 60$\%$ of 102 students said: \textit{I feel less anxious coming to class after doing an EoL assignment.}

\item 52$\%$ of 98 students said: \textit{I look forward to coming to class more after completing an EoL assignment.}
\end{itemize}

\noindent After the midterm test the students said the following:

\begin{itemize}

\item 79$\%$ of 159 students said: \textit{I think the EoL assignments helped me be more prepared for the midterm.}

18$\%$ of the 159 students said: \textit{the EoL assignments did not make a difference [in preparing them for the midterm].}

\item 54$\%$ of 142 students said: \textit{I feel like the EoL assignments influenced my confidence for the midterm in a positive way.}

25$\%$ of the 142 student said: \textit{the EoL assignment did not influence [my confidence for the midterm] at all.}

\item 37$\%$ of 159 students said: \textit{I feel like the EoL assignments helped reduce my test anxiety for the midterm.}

23$\%$ of the 159 students said: \textit{I do not have any sort of test anxiety.}
\end{itemize}

\noindent After the final exam, the students said the following:

\begin{itemize}
\item 90$\%$ of 119 students said: \textit{I enjoyed the [MATH 235] class.}

\item 90$\%$ of 124 students said: \textit{I think the EoL assignments helped me learn and keep up with the material regardless of liking the EoL assignments.}

\item 61$\%$ of 115 students said: \textit{I enjoyed the EoL assignments.}

\item 66$\%$ of 122 students said: \textit{the EoL assignments helped me structure my study habits for MATH 235.}

\item 45$\%$ of 93 students said: \textit{I have applied such study habits to other classes or am planning to do so next term.}

\item 68$\%$ of 130 students said: \textit{I understood more of the material in MATH 235 compared to my other classes.}

27$\%$ of the 130 students said: \textit{there was no difference in understanding compared to my other classes.}

\item 45$\%$ of 109 students said: \textit{the EoL assignments helped my mental well-being.}

27$\%$ of the 109 students said: \textit{ I do not know if the EoL assignments helped my mental well-being.}
\end{itemize}

\noindent \subsection{Descriptive Statistics}

At the end of term, the marks of all 459 students who took and completed MATH 235 in Fall 2014 were analyzed. I hypothesized that the Fall 2014 group would outperform other MATH 235 sections that instead completed weekly assignments rather than EoL assignments. To test against my hypothesis, I also analyzed the marks of 387 students who took and completed MATH 235 in the Fall 2013 and Spring 2014 terms. For all three of these terms I was one of the instructors and the same material was taught. In particular, I had two sections in both fall terms, and one section in the Spring 2014 term. However, one of my sections from Fall 2013 and Spring 2014 were not included in the sample of 387 students. This is due to the fact that my sections completed EoL assignments whereas the other sections had weekly assignments consisting of the same questions. For convenience, I will refer to the Fall 2014 group as the \textit{EoL assignment group} and the remaining group as the \textit{weekly assignment group}.

\begin{table}[h]
\begin{center}
\begin{tabular}{ lc c c c c}
\hline
& Assign. Average & Midterm Mark & Final Mark & Overall Grade \\ \hline
Mean &72.93005181 &70.39896373 & 59.87823834 & 66.55958549 \\
Standard Error & 1.141249365 &1.002159567 &1.015469067 & 0.915371757 \\
Median & 79 & 71 & 62 & 67\\
Mode & 96 & 75 & 64 & 61 \\
Standard Deviation & 22.42199242 & 19.68931146 & 19.95080164 &17.98420155\\
Sample Variance &502.7457439 & 387.6689859 & 398.0344862 & 323.4315053\\
Kurtosis &1.251196139 & 2.134767813 & 0.636438169 & 0.560222138\\
Skewness &-1.290170072 & -1.078715037 & -0.675483433 & -0.634431278\\
Confidence Level (95$\%$) & 2.243861437 & 1.970390762 & 1.996559165 & 1.799753366\\ \hline \\
\end{tabular}
\caption[]{Descriptive statistics of the weekly assignment group.}
\end{center}
\end{table}

For both groups of students, the \textit{marks composition data set}, i.e., the breakup of the evaluation method, consisted of the students: assignment average, midterm mark, final mark and overall grade in MATH 235. The assignments were cumulatively worth 10$\%$ of the students ' overall grade. The EoL assignment group was assigned 32 assignments throughout the term, whereas the weekly assignment group was given 10 assignments. The midterm and final exam were worth 25$\%$ and 65$\%$ of the overall grade, respectively. The descriptive statistics of the marks data set for the weekly assignment group and EoL assignment group can be seen in Tables 1 and 2, respectively.

\begin{table}[h]
\begin{center}
\begin{tabular}{ lc c c c c}
\hline
& Assign. Average & Midterm Mark & Final Mark & Overall Grade \\ \hline
Mean & 68.07259821& 74.46623094 &70.81481481& 71.31633636\\
Standard Error& 0.990440257&0.851981117&0.810359931&0.724626167\\
Median&74.64516129&77&73&73.63225806\\
Mode & 94.38709677&95& 66&74.3483871\\
Standard Deviation&21.21947462&18.25308651&17.36138234&15.52459772\\
Sample Variance&450.2661033&333.1751672&301.4175966&241.0131345\\
Kurtosis&0.62697467&3.030615182&1.631463296&0.715859109\\
Skewness&-1.082880529& -1.391674584&-0.960573545&-0.772909234\\
Confidence Level (95$\%$) &1.946370622&1.674276671&1.592484504&1.424004197\\ \hline \\
\end{tabular}
\caption[]{Descriptive statistics of the EoL assignment group.}
\end{center}
\end{table}

We can see that the weekly assignment group had a higher assignment average than the EoL assignment group. However, the EoL assignment group's mean for the midterm exam, final exam and overall grade was higher than the weekly assignment group. The standard error, median, mode, standard deviation, sample variance, Kurtosis value, skewness and confidence level is also given in the tables for both groups. To be scientific, it is important to point out the possibility that the EoL assignment group could have been a stronger group of students academically. However, such a possibility can only be investigated if more data is collected and analyzed. It is also important to mention that many students (roughly 12 - 40) from both the weekly and EoL assignment groups that were enrolled in other sections attended my lectures or office hours. Additionally, if students from the other sections requested to use my personal notes prepared for my students, they were given access to these resources. Approximately 15$\%$ - 20$\%$ of students from other sections used my personal notes.

\begin{figure}[h]
\centering
\subfigure[Weekly assignment group.]{
\includegraphics[scale=0.48]{MH_hist_weekly.png}}
\subfigure[EoL assignment group.]{
\includegraphics[scale=0.48]{MH_hist_eol.png}}
\caption[]{Histogram of the students' overall grade for weekly assignment group (left) and EoL assignment group (right).}
\end{figure}

The corresponding histogram of the weekly and EoL assignment groups' overall grades can be seen in Fig. 5a and 5b, respectively. A histogram with a normal distribution fit of both groups' overall grade can be seen in Fig. 6. These figures illustrate that there is a skew to the right for both groups of students. However, the EoL assignment group seems to be more skewed than the weekly assignment group. This comparison is depicted in Fig. 6.

\begin{figure} [ht]
\centering
\includegraphics[scale=0.45]{MH_hist_normal.png}
\caption[]{Histogram with a normal distribution fit of weekly and EoL assignment groups' overall grade data.}
\label{fi:1}
\end{figure}

One can also see a set of four box plots for the weekly assignment group and EoL assignment group in Fig. 7a and 7b, respectively. In each plot, the boxes graphed from left to right represent the students' assignment average, midterm mark, final mark and overall grade obtained in MATH 235, respectively. Notice that the assignment box plots for both groups in Fig. 7 is similar. However, the midterm, final and overall grade box plots and their corresponding relationships to each other differ. In particular, it appears that the EoL assignment group has more consistent, evenly distributed grades compared to the weekly assignment group. The width (i.e., variation of range) of the box plots for the midterm, final and overall grade match each other rather well for the EoL assignment group. Yet, for the weekly assignment group, the ranges of the box plots vary more, especially with respect to each other. Also, the assignment average box plot for the weekly assignment group consists of the highest grades relative to other box plots. However, this is not the case for the EoL assignment group. Since the assignment average box plots for both Fig. 7a and 7b are similar, yet there is more variation between the boxes in Fig. 7a rather than Fig. 7b, this may suggest that the frequency at which the students had to complete assignments in the EoL assignment group allowed the students to retain more information.

\begin{figure}[h]
\centering
\subfigure[Weekly assignment group with sample size of $n=387$.]{
\includegraphics[scale=0.48]{MHP_Weekly_Group_Boxplot.png}}
\subfigure[EoL assignment group with sample size of $n=459$.]{
\includegraphics[scale=0.48]{MHP_EoL_Group_Boxplot.png}}
\caption[]{Box plots from left to right of students' assignment average, midterm mark, final mark and overall grade.}
\end{figure}

\noindent \subsection{Further Plots}

In Fig. 8, the final mark that each student obtained was compared with their assignment marks. Recall that the weekly assignment group had to complete 10 assignments throughout the term, whereas the EoL assignment group was given 32 assignments to complete. Due to this frequency variation of assignments, Fig. 8a depicts the percentage of assignments the students completed versus the students' overall mark in MATH 235. Fig. 8b plots the students' assignment average against their overall grade in MATH 235. The line of best fit for the two groups in both Fig. 8a and 8b is also depicted.

In both the left and right plots in Fig. 8, a positive correlation can be seen. In the right plot, the line of best fit for both groups has a similar slope. However, the blue line (representing the EoL assignment group) is shifted by approximately $+10 \%$ relative to the red line (representing the weekly assignment group). In Fig. 8a, the blue line of best fit is similar to the blue line shown in Fig. 8b. However, the slope of the red line of best fit in Fig. 8a is different compared to the slope of the red line in Fig. 8b. In Fig. 8a, the red line is not as steep as the blue line. In particular, the slopes seen in Fig. 8a, are approximately $\frac{1}{3}$ and $\frac{1}{5}$ for the blue and red line, respectively. The blue line is always above the red line after the one intersection point at approximately $(35, \; 50)$. Hence, if a student from the EoL assignment group only completed 35$\%$ of the assignments, they were more likely to obtain a higher grade than a student in the same situation in the weekly assignment group. Note that 35$\%$ assignment completion for the EoL assignment group corresponds to approximately 11 assignments, which is roughly the cardinality of the entire assignment set for the weekly assignment group.

\begin{figure}[h]
\centering
\subfigure[Assignments completed vs overall mark.]{
\includegraphics[scale=0.52]{MHP_Assignment_Average_vs_Overall_Mark_Scatter_Plot_and_Line_Of_Best_Fit.png}}
\subfigure[Assignment average vs overall mark.]{
\includegraphics[scale=0.52]{MHP_Assignment_Average_vs_Overall_Mark.png}}
\caption[]{Data comparing information about students' assignments and overall mark.}
\end{figure}

Therefore, the data seems to suggest that the higher the frequency of evaluation (not necessarily dependent on marks), the better the measurement of the students' understanding. This demonstrates a decoupling of the assignment average from the overall grade for the EoL assignment group, hence giving a more accurate measure of the students' progress. More precisely, by comparing the information in Fig. 8a and 8b, it appears that the actual assignment grades for the EoL assignment group does not necessarily indicate that the students' understanding is stagnant, or that the numerical value properly characterizes the students' competency. One of the goals of the EoL assignments was for the students to make their mistakes on the assignments rather than on the exams. Thus, even if a student's assignment average is not very high, they had the opportunity to correct and learn from their mistakes at a faster rate due to the higher frequency of evaluation and feedback. This observation also seems consistent from the box plots shown in Fig. 7b.

The final observation from Fig. 8a I will discuss is potentially related to mental health. Notice there are several students who completed less than $50\%$ of their assignments, but still received a decent grade in the course. In particular, there were 38 students in the weekly assignment group and 13 students in the EoL assignment group that completed at most $50\%$ of their assignments but still obtained a final grade of over $65\%$ (which is the average students need to maintain in the Math Department at the University of Waterloo). Interestingly, it appears from Fig. 8a that it is less likely to have such `outliers' in the EoL assignment group. Perhaps this has to do with helping the students structure a study schedule, or giving them the opportunity to seek help more often. Nonetheless, I wanted to examine the number of these `outliers' because several of them (not all) are students with poor mental health. I can say this with confidence because I have spoken to several of these students (also see Appendix B for some commentary from the students).

Overall, there are a variety of mental health issues or environmental reasons (e.g., domestic violence or sexual violence) that could result in students not attending class or submitting work on time. However, this does not imply that students intend for such things to happen, nor does it imply that the students are not working on their own. Thus, it is very possible to have several outliers, i.e., students who do not submit assignments but who still do well in classes, because they work without focusing on deadlines. For instance, some of the students I spoke with suffer from depression. It is common for students with depression to spend more time in bed, especially during conventional work hours. Recent studies have revealed that ``individuals with depressed mood show as much as a 12$\%$ reduction in memory compared to individuals without depressed mood when depressive thoughts are present, but perform similarly to individuals without depressed mood when depressive thoughts are not present'' \cite{15b}. Therefore, if a student is having a terrible day, it very well may be more productive for the student to attend to their mental health and then work later on their own. That being said, if a student does struggle with mental health and finds that they continually cannot keep up with their studies, such a student should carefully reflect on their course load. It is vital for students to understand their own limitations and adapt new strategies or schedules accordingly. Even if professors do everything in their power to help or accommodate students with poor mental health, it is still up to the student to communicate with their instructors, and take hold of their own responsibilities, lifestyle and needed accommodations.

\noindent \subsection{Statistical Analysis: Analysis of Variance}

An exhaustive statistical analysis comparing the EoL assignment group and weekly assignment group can be found in Appendix A. It was strongly suggested to me by professional statisticians that considering too many statistical tests is not a good idea because eventually, an issue will arise. I am not sure I agree with such a statement. In my view, if a reasonable data set is collected, it seems sensible to think that the statistical tests would give consistent results. However, I am not a statistician, so perhaps this is a naive notion. Nonetheless, for full disclosure purposes please refer to Appendix A to review all the tests I performed on the data.

For the statistical tests presented in this section, a simple random sample of 100 individuals from both groups was taken 16 separate times and averaged. I will refer to these new populations as the \textit{random EoL assignment group} and the \textit{random weekly assignment group}. These data sets were also graphed as a set of box plots which can be seen in Fig. 9. The structure of Fig. 9 is similar to Fig. 7. Notice that for each group, the relationship between the box plots depicting the students' assignment average, midterm mark, final mark and overall grade obtained in MATH 235 exacerbates the trends seen in Fig. 7. This observation seems to further support the claims made above while analyzing Figs. 7 and 8.

\begin{figure}[h]
\centering
\subfigure[Randomly sampled weekly assignment group.]{
\includegraphics[scale=0.48]{MHP_Weekly_Group_Random_Sample_Boxplot.png}}
\subfigure[Randomly sampled EoL assignment group.]{
\includegraphics[scale=0.48]{MHP_EoL_Group_Random_Sample_Boxplot.png}}
\caption[]{Box plots from left to right of students' assignment average, midterm mark, final mark and overall grade.}
\end{figure}

The random weekly and EoL assignment groups' data was also used for further statistical analysis. An analysis of variance (ANOVA) test was used to analyze the statistically significant differences between the two groups. In particular, the one-way ANOVA test was performed on the two mutually independent groups in which the means of the two groups were compared. Thus, the null hypothesis, denoted by $H_0$, was defined as follows:
\begin{itemize}

\item[$H_0:$] there is no statistically significant difference between the randomly sampled populations, against the alternative that there is a statistically significant difference.
\end{itemize}

\noindent The one-way ANOVA test was performed separately for the four different categories of randomized data, i.e., the randomly averaged assignment average, midterm mark, final mark and overall grade data. A test decision for $H_0$ is returned as well as a table with seven columns. The first column shows the source of variability; the second depicts the sum of squares, denoted by SS, due to each source; the third shows the degrees of freedom, denoted by df, associated with each source; the fourth contains the mean squares, denoted by MS, which is the ratio $\frac{SS}{df}$; the fifth shows the F-statistic, denoted by F, which is the ratio of the mean squares; the sixth column contains the p-value, denoted by p; and finally, the seventh column shows the critical F-value, denoted by $F_{\text{crit}}$. If $F>F_{\text{crit}}$, which corresponds to a very small p-value, then the null hypothesis $H_0$, is rejected. The results can be seen below in Table 3.

\begin{table}[h]
\begin{center}
\begin{tabular}{ lc c c c c c c c}
Results for sampled assignment average data: \\
\hline
Source of Variation & SS & df & MS & F & p & $F_{crit}$ \\ \hline
Between groups & 216.082 & 1 & 216.082&59.35 & 1.36e-08 & 4.171 \\
Error (within group) &109.229 & 30 & 3.641 & & \\
Total & 325.311& 31 \\
\hline \\
Results of sampled midterm mark data. \\
\hline
Source of Variation & SS & df & MS & F & p & $F_{crit}$ \\ \hline
Between groups & 102.030 & 1 & 102.030 &42.72& 3.15e-07 & 4.171 \\
Error (within group) &71.651& 30 & 2.388 & & \\
Total & 173.682& 31 \\
\hline \\
Results of sampled final mark data. \\
\hline
Source of Variation & SS & df & MS & F & p & $F_{crit}$ \\ \hline
Between groups & 297.253 & 1 & 297.253 &136.24& 1.11e-12 & 4.171 \\
Error (within group) &65.454& 30 & 2.388 & & \\
Total & 362.707& 31 \\
\hline \\
Results of sampled overall grade data. \\
\hline
Source of Variation & SS & df & MS & F & p & $F_{crit}$ \\ \hline
Between groups & 127.932 & 1 & 127.932 & 81.27& 4.84e-12 & 4.171 \\
Error (within group) &47.226 & 30 & 1.574 & & \\
Total & 175.158 & 31 \\ \hline
\end{tabular}
\caption[]{Results of the one-way ANOVA tests.}
\end{center}
\end{table}

Since $p \leq 3.15e-07$ ($F>F_{\text{crit}}$) for all the ANOVA tests reported in Table 3, $H_0$ is rejected for the randomized assignment average ($p=1.36e-08<0.05$), midterm mark ($p=3.15e-07<0.05$), final mark ($p=1.11e-12<0.05$) and overall grade ($p=4.84e-12<0.05$) data. Hence, the differences in the means of the two populations are statistically significant. This result also agrees with the t-test results for each pair of means (which can be found in Appendix A).

A two-way ANOVA test was also performed for all the random EoL and weekly assignment group data. This test is useful when data can be classified along two different dimensions. One dimension will be the different randomized groups, i.e., the random weekly and EoL assignments groups; the other will be the different ways the students were evaluated, i.e., the assignment, midterm, final and overall marks for each group. The two-way ANOVA test will allow us to test with 95$\%$ confidence the three following hypotheses:
\\
\begin{itemize}
\item[$H_{0_1}:$] there is no statistically significant difference between the population means of the two different groups,
\item[$H_{0_2}:$] there is no statistically significant difference in the effectiveness of how the populations were taught or evaluated, and their performance in the class,
\item[$H_{0_3}:$] there is no statistically significant difference between the interaction of the groups and performance/evaluation. \\
\end{itemize}

The results of the two-way ANOVA test can be seen in Table 4. Since the p-value corresponding to the different groups (EoL and weekly) is $p=6.05e-26<0.05$, the null hypothesis $H_{0_1}$ is rejected. This implies that the different groups of students are statistically different. The p-value corresponding to the different types of course evaluation (i.e., assignments, midterm, final and overall grade) is $p=5.26e-12<0.05$, and so the null hypothesis $H_{0_2}$ is also rejected. Hence, there is a statistically significant difference between the effectiveness of how the different groups where taught or evaluated (or how the course was structured), and the outcome (i.e., marks) of the two groups. The two-way ANOVA test cannot say what the difference is, but whatever the difference may be, it significantly effected the performance or understanding of the two different groups of students.

\begin{table}[h]
\begin{center}
\begin{tabular}{ lc c c c c c c c}
\hline
Source of Variation & SS & df & MS & F & p & $F_{crit}$ \\ \hline
Group & 143.437 & 1 & 143.437 & 58.633& 5.36e-12 &3.920\\
Evaluation types &504.637 & 3&168.212 & 68.761 & 6.05e-26 & 2.680\\
Interaction & 599.861 & 3 &199.954 & 81.736 &7.21e-29 &2.680\\
Error & 293.561 & 120 &2.446 &\\
Total & 1541.495 &127 \\ \hline
\end{tabular}
\caption[]{Results of the two-way ANOVA test for random overall grade data.}
\end{center}
\end{table}

The p-value of the interactions is $p=7.21e-29<0.05$ and so, the null hypothesis $H_{0_3}$ is rejected. We conclude there are statistically significant differences in the interaction between the different types of course evaluation (i.e., assignment, midterm, final and overall marks) and the groups of students. We can plot the mean interactions to examine the interactions more carefully. These plots can be seen below in Fig. 10. Lines that are roughly parallel are indications of the lack of interaction, while lines that are not roughly parallel indicate interaction.

\begin{figure}[h]
\centering
\subfigure[Interaction factor 1.]{
\includegraphics[scale=0.44]{MH_Inter_1.png}}
\subfigure[Interaction factor 2.]{
\includegraphics[scale=0.41]{MH_Inter_2.png}}
\caption[]{The mean interactions between the levels of the two factors}
\end{figure}

Therefore, in summary, these are incredibly small p-values indicating that the null hypotheses do not seem appropriate. This strongly indicates that the groups of students were not comparable. It appears there was a fundamental difference between how the students from the weekly assignment group and EoL assignment group absorbed and understood the MATH 235 material. However, in order to make any further statements, more data would have to be collected and analyzed.

\noindent \subsection{Conclusions for the restructuring of assignments}

In conclusion, by comparing Figs. 5 to 10 and the obtained p-values for the ANOVA tests, it appears that there is a statistically significant difference between the performance of the weekly assignment group and EoL assignment group. In addition, the results depicted in Figs. 7 to 9 seem to support the concept illustrated in Ebbinghaus' forgetting curve. That is, if students are able to review and test their knowledge on a more frequent basis, they will better understand the material and thus, perform better in classes. Perhaps if the math assignments, at least for certain courses, were structured in a similar fashion as the EoL assignments rather the weekly assignments, the students may obtain higher grades. Furthermore, from the survey results, we can see that at least for a portion of the students, such a structure may increase confidence, reduce anxiety and stress.

\noindent \section{Other Considerations}

After reviewing the survey results discussed above, I spent much time talking with students individually and in small groups. All the feedback inspired the suggestions that I will now present. As a note to the reader, one suggestion is implementing better ways to integrate technology into the education system. I think this is an important aspects that needs serious consideration, and so it will be the focus of the proceeding paper.

\noindent \subsection{Basic Mental Health Training and Fall Reading Week}

During orientation week training at the University of Waterloo, orientation leaders are told to expect one out of two students to have some sort of mental health issue. Orientation leaders are also told to be alert for suicidal students. This high value suggests that mental health training for professors would be beneficial. Recall that from the student surveys, out of 107 students, 72$\%$ think professors should have such training. Desperate, mentally vulnerable people have the potential to be a danger to themselves or others; the faculty needs to know how to deal with such situations.

A wide variety of Canadian universities and colleges have conducted surveys showing that poor mental health, including suicidal thoughts, attempts and overall rates, is becoming a larger issue \cite{8}. One proposal to counteract this pandemic is to introduce a fall reading week. This is a solution that could be easily implemented; and in fact, many colleges and universities have already done so. For instance, Flemming College, Brock University, Wilfrid Laurier University, and Queen's University all now have fall reading weeks \cite{9}.

A recent survey at the University of Waterloo showed that 74$\%$ of over 6,000 students polled want a fall reading week. However, according to university spokesman Nick Manning: ``implementing a reading week at the University of Waterloo is complicated by the demands of our co-op education program. So it's not as straight forward as the other 14 universities in Ontario who have been able to implement a reading week" \cite{10}. Nonetheless, the University of Waterloo senate will be voting hopefully this year on the perspective of having a fall reading week \cite{11}. Although, this was done approximately ten years ago, and the referendum never passed. A fall reading week at the University of Waterloo would only require two additional school days to be added to the fall term \cite{10}.

\noindent \subsection{Mental Health Days}

Another suggestion would be for the university to consider students' mental health in the same fashion as physical health, i.e., universities and colleges could employ mental health days like we do sick days. For instance, if a student is physically ill, especially with the influenza virus, alternative arranges or exceptions for courses are usually made for the student. Until recently, in the case of influenza, the University of Waterloo created a special form when H1N1 was at its peak. Students could submit the `declaration of absence for influenza-like illness' form online excusing them from class and the corresponding missed work without being penalized \cite{16}. Consider if a student for instance was going through a suicidal episode, a similar form could be equally beneficial, humane and perhaps even reduce drop-out rates \cite{17}. Upon surveying my students I found that:

\begin{itemize}
\item 71$\%$ of 102 students said: \textit{there should be mental health days like there are sick days.}

\item 67$\%$ of 95 students said: \textit{I think mental health days would be useful in a university setting.}
\end{itemize}

I personally hypothesize that mental illness is like the seasonal flu - anyone can have it or contract it. If this conjecture turns out to have enough supporting evidence, implementing mental health days especially seems appropriate. Interestingly, in a survey, students were asked if they agreed with such a hypothesis and 76$\%$ of 45 students agreed.

\begin{figure} [ht]
\centering
\includegraphics[scale=0.54]{MH_pie3.png}
\caption[]{Distribution of student responses (111 answered this question).}
\label{fi:1}
\end{figure}

A typical concern regarding the implementation mental health days is that they will be abused far more than sick days. Thus, in a survey, I asked my students what their view were on the potential abuse of mental health days. Fig. 11 shows the distribution of the 111 students who responded. One can see that 44$\%$ of the students said mental health days are more likely to be abused than sick days; whereas 42$\%$ did not think there would be a difference. My hypothesis is that mental health days would be abused the same amount as sick days. As some students have pointed out in the survey:
\begin{itemize}
\item \textit{They would be abused the same amount as sick days, in some cases, due to the physical effects of certain mental illnesses, mental health days $==$ sick days.}
\item \textit{I think mental health is the same as ``sick", both can make it difficult to perform at your best.'}
\end{itemize}

\noindent\subsection{A Hypothesis of Mine - The Importance of Socialization}

One complaint that I often hear from students (and agree with) is that ``many mathematicians are socially-inept." The comparison of STEM individuals to those of autism-spectrum disorder (e.g., Asperger syndrome) is not uncommon, especially by the student body. Typically, students label and associate `social-ineptness' when their peers or professors seem to lack empathy, have abnormal social interaction(s), minimal social interactions(s), or scant communication(s) \cite{20}. Although this may initially sound offensive, it is supported by modern research.

For example, one study was conducted in 2001 with a group of 840 students from Cambridge University (454 males, 386 females) where each student was given an autism-spectrum quotient test (AQ) \cite{AQ}. This test attempts to measure autistic traits for `individuals of normal intelligence.' The researchers found that scientists obtained higher AQ scores compared to non-scientists. In particular, it was found that within scientific fields, mathematicians had the highest AQ scores. Physical scientists, computer scientists and engineers also had higher scores compared to students in `human or life-centered' sciences (specifically medicine, veterinary science and biology) \cite{AQ}. These results also agreed with an earlier study carried out in 1998 again at Cambridge University. It was found that autism occurred ``significantly more often in families of students in the fields of physics, engineering, and mathematics'' \cite{AM}.  Overall, ``it seems that people who have autistic traits may find it easier or more natural to work in these fields," said Simon Baron-Cohen; the co-director of the Autism Research Centre at Cambridge University \cite{21}.

I found these studies to be quite intriguing since it has also been revealed that on average, millennials spend around seven to eleven hours a day in front of some sort of computer screen (desktop, tablet, smart phone etc.) \cite{22}. This indicates that the younger generations are not socializing (face-to-face) nearly as much as previous generations. Hence, children now do not gain as much experience physically socializing with others, which can be quite problematic. Learning how to read body language, especially facial expressions, is an imperative part of physiological development \cite{23}. Facial expressions for instance are a form of nonverbal communication that conveys emotion.  In fact, facial expressions are one of the primary ways that emotion is communicated between individuals \cite{24}. Being able to witness a person's reaction activates mirror neurons, allowing one to simulate what another is experiencing \cite{25}.

A part of the brain that plays a key role in facial recognition is the amygdala \cite{28}. As we currently understand, the amygdala which is part of the limbic system in the brain, processes emotional reactions, memory formation and decision making \cite{26}. It plays an essential role in storing and forming memories associated with emotional experiences \cite{27}, and researchers speculate that the amygdala may be highly involved in recognizing negative emotions \cite{29}. Individuals with a larger amygdala are better able to recognize the emotional state of others through facial expressions. Therefore, researchers hypothesize that larger amygdalae may be an indication of greater emotional intelligence. The development of the amygdala is correlated with involvement in social networks. More specifically, the size of the amygdala increases according to the number of social interactions a person has. The variations and complexity of different social groups also influences development \cite{30}, \cite{31}.

In addition, several studies have indicated that the amygdala has a significant role in mental states and hence, poor mental health \cite{32}. For instance, some studies have shown that children with anxiety disorders usually have a smaller left amygdale. An overactive left amygdala seems to be associated with borderline personality disorder and depression; whereas an underactive left amygdala seems to be correlated with psychopathy. In general, the amygdala seems to be connected with social anxiety, obsessive compulsive disorders, social phobias, bipolar disorder and not surprisingly, autism-spectrum disorders \cite{32}, \cite{33}, \cite{34}, \cite{35}. Studies have demonstrated that people with autism-spectrum disorders have a difficult time identifying facial expressions. In particular, some research has suggested that such disorders escalate due to a lack of physical social experience: ``the development of face perception and social cognitive skills are supported by the amygdala–fusiform system, and that deficits in this network are instrumental in causing autism'' \cite{36}.

Due to these studies and findings, I hypothesize that all mathematicians are ``a little bit autistic" and potentially, more susceptible to poor mental health. Consider that to be a successful mathematics student, the student has to spend more time alone doing problems compared to individuals studying less technical fields. To complete an undergraduate degree in mathematics, at least four years has to be dedicated to learning. Typically such an endeavour starts around the age of 17 - 19, sometimes even younger. If a student then enrolls into a masters or doctorate degree, that is another two to six years of intensive work. In total (not including post-doctoral fellowships), that is approximately four to ten years during the period in one's life where individuals grow into the people they want to be \cite{Jay}. This is a lot of time being spent alone. This may imply that the parts in the brain responsible for empathy and socialization, such as the amygdala, will not sufficiently develop; especially compared to people with more complex social networks and interactions.

If this hypothesis is true, we may have a serious problem. If a STEM background is necessary in the twenty-first century, a time where computers and the Internet are an integral part of our lives, will people start becoming less socially aware and less empathetic (for some more information see \cite{37})? Could the people who are needed to solve the world's biggest problems be more likely to suffer from psychological disorders? These are questions that we cannot afford to ignore or avoid.

This hypothesis also seems to be connected with a major complaint people, and especially women have regarding mathematics or generally STEM education: it is lonely. Social interaction is more important biologically for women than it is for men \cite{38}, and women are more likely to change academic and professional paths than men due to their lack of confidence \cite{CBC}. Perhaps these are reasons why there are more men in fields like maths, physics and engineering. A potential solution may sound simple: make mathematics and the education system more social with the help of technology (which will be the focus of the next paper). In my opinion, we need to create an open, positive, social learning environment to nurture healthy educated people.

\noindent\section{Conclusions}

Mathematics is the only subject where proofs exist. As Chris Rock said: ``there is mathematics and everything else is debatable.'' It is a beautiful subject that is based in truth without ambiguity. However, to learn mathematics is no trivial task. It requires an excellent teacher, and a good state of mental health so the mind's computational power can be allotted to learning. How best to learn and teach mathematics is debatable and generally is gaining more attention by researchers. Nonetheless, it is essential to listen to the students' needs and suggestions. Providing a positive, open and safe learning environment is necessary. In the post-secondary education system, this requires a lot of time and effort from the university and in particular, its faculty. However, that is the purpose of post-secondary institutions: to be a place of higher learning, open-mindedness and innovation. We need to actively prepare individuals to create and maintain a dynamic, functional, altruistic society that has many challenging problems to solve within the next 50 - 100 years \cite{4}. Without awareness of mental health, this goal is very unlikely to be achieved.

\noindent\section{Acknowledgments}

I would like to thank all of the students who participated in these surveys. Thanks for your openness and honesty. Also, I would like to thank the following undergraduate math students for their discussions, ideas, questions, interpretations and suggestions: Amar Bhalla, Andre Jodat-Danbrani, Andrew Park, Ashna Singh, Glen Chalatov, Ivan Li, Jessie Huang, Jessie Yeung, Joshua Thompson, Cheyenne Guo, Varmanaa, Mohammad Ali, Monica Pierce and all others who shared their stories.