\section{Introduction}

Let $G$ be a finite group and $H$ be a subgroup of $G$. We denote by $o_H(x)$ the order of an element $x\in G$ relative to $H$, i.e. the smallest positive integer $m$ such that $x^m\in H$. Also, $C_n$ denotes the cyclic group with $n$ elements, where $n$ is a positive integer.

The main objective of \cite{9} was to study the quantity
$$\psi_H(G)=\sum\limits_{x\in G}o_H(x),$$
which is a generalization of the sum of element orders of $G$, that is
$$\psi(G)=\sum\limits_{x\in G}o(x),$$
where $o(x)$ denotes the usual order of an element $x\in G$. The second quantity appeared in \cite{1} and was further studied in several papers such as \cite{2,3,4,8}. The main result of \cite{1} states that the maximum value of the sum of element orders among all finite groups of order $n$ is $\psi(C_n)$. In other words, we have:\\

\textbf{Theorem 1.1.} \textit{Let $G$ be a group of order $n$. Then $\psi(G)\leq \psi(C_n)$ and the equality holds if and only if $G\cong C_n$.}\\

Up to the class of finite nilpotent groups, in \cite{9}, it was proved that a similar property occurs if one works with element orders relative to a subgroup. More exactly, the following result holds:\\

\textbf{Theorem 1.2.} \textit{Let $G$ be a nilpotent group of order $n$ and $H$ be a subgroup of order $m$ of $G$. Then
\begin{align}\label{r1}
\psi_H(G)\leq \psi_{H_m}(C_n),
\end{align}
where $H_m$ is the unique subgroup of order $m$ of $C_n$.}\\

Moreover, the author conjectured that Theorem 1.2 is valid for any finite group. However, in this paper we show that this statement is false by outlining a way to construct a finite group $G$ such that inequality (\ref{r1}) is not satisfied for at least one subgroup $H$ of $G$. For the ease of writing, we introduce the following ratio
$$\psi'_H(G)=\frac{\psi_H(G)}{\psi_{C_{|H|}}(C_{|G|})}.$$
Using this notation, our main result may be stated as follows:\\

\textbf{Theorem 1.3.} \textit{There are infinitely many finite groups $G$ having a subgroup $H$ such that $\psi'_{H}(G)>1$.}\\

A proof of Theorem 1.3 along with other several properties related to this generalization of the sum of element orders of a finite group are highlighted in the following section.

\section{Proof of the main result and new properties of $\psi_H(G)$}

As a starting point, we outline a preliminary result which puts together several known facts concerning the sum of element orders of a finite group $G$ and its generalization. For an element $x\in G\setminus H$, a novelty would be an upper bound for $o_H(x)$. As a consequence, we also determine a new upper bound for $\psi_H(G)$.\\

\textbf{Lemma 2.1.} \textit{Let $G$ be a group of order $n=p_1^{\alpha_1}p_2^{\alpha_2}\ldots p_k^{\alpha_k}$, where $p_1<p_2<\ldots<p_k$ are primes and $\alpha_i$ is a positive integer, for all $i\in\lbrace 1, 2,\ldots, k\rbrace$, $H$ be a subgroup of $G$ of index $q$, $K$ be a normal subgroup of $G$ and $(G_i)_{i=\overline{1,s}}$ be a family of finite groups of coprime orders. Then:
\begin{itemize}
\item[i)] $\psi(C_n)=\prod\limits_{i=1}^k\frac{p_i^{2\alpha_i+1}+1}{p_i+1};$ (see Lemma 2.9 (1) of \cite{4})
\item[ii)] $\psi(C_n)\geq\frac{p_1}{p_k+1}n^2;$ (see Lemma 2.9 (2) of \cite{4})
\item[iii)] $\psi_K(G)=|K|\psi(\frac{G}{K})$; (see p. 628 of \cite{9})
\item[iv)] $\psi_{H_1\times H_2\times\ldots\times H_s}(G_1\times G_2\times\ldots\times G_s)=\prod\limits_{i=1}^s\psi_{H_i}(G_i)$, where $H_i$ is a subgroup of $G_i$, for all $i\in\lbrace 1,2,\ldots, s\rbrace$; (see Lemma 2.1 a) of \cite{9})
\item[v)] $o_H(x)\leq q, \ \forall \ x\in G\setminus H;$
\item[vi)] $\psi_H(G)\leq |H|(q^2-q+1).$
\end{itemize}}
\textbf{Proof.} For item \textit{v)}, we have
$$o_H(x)=\frac{|\langle x\rangle|}{|\langle x\rangle\cap H|}=\frac{|H\langle x\rangle|}{|H|}\leq \frac{|G|}{|H|}=q, \ \forall \ x\in G\setminus H.$$

The above inequality is further used to prove the result given by item \textit{vi)}. More exactly, one can write
$$\psi_H(G)=|H|+\sum\limits_{x\in G\setminus H}o_H(x)\leq |H|+q(|G|-|H|)=|H|(q^2-q+1).$$
\hfill\rule{1,5mm}{1,5mm}\\

Before proving our main result, we recall that a subgroup $H$ of a finite group $G$ is isolated in $G$ if, for every element $x\in G$, we have either $x\in H$ or $\langle x\rangle\cap H=1$. For more details about these subgroups, we refer the reader to \cite{5}. Also, in \cite{10}, the isolatedness property was characterized as follows: $H$ is isolated in $G$ if and only if $\psi_H(G)=|H|+\psi(G)-\psi(H).$ \\

\textbf{Proof of Theorem 1.3.} Let $p$ be a prime and $r\geq 3$ be an  integer. By following Example 7 of \cite{7}, starting with the field $\mathbb{F}_q$ with $q=p^r$ elements, we can construct a Frobenius group $\overline{G}=\overline{K}\rtimes \overline{H}$, whose kernel $\overline{K}$ is isomorphic to the additive group of $\mathbb{F}_q$, while $\overline{H}$ is isomorphic to the multiplicative group of $\mathbb{F}_q$. Hence, $\overline{K}\cong C_p^r$ and $\overline{H}\cong C_{p^r-1}$. Moreover, since $\overline{H}$ is isolated in $\overline{G}$ and its conjugates trivially intersect, we obtain
$$o_{\overline{H}}(\overline{G})=|\overline{H}|+\psi(\overline{G})-\psi(\overline{H})=(p^r-1)(\psi(\overline{H})+p)=(p^r-1)(\psi(C_{p^r-1})+p).$$
Further, suppose that $p=2$ and $r$ is chosen such that $2^r-1$ is a prime (i.e. $2^r-1$ is a Mersenne prime). Using Lemma 2.1 \textit{i)}, \textit{iii)}, we get
$$\psi'_{\overline{H}}(\overline{G})=\frac{o_{\overline{H}}(\overline{G})}{o_{C_{|\overline{H}|}}(C_{|\overline{G}|})}=\frac{(2^r-1)(\psi(C_{2^r-1})+2)}{(2^r-1)\psi(C_{2^r})}
=\frac{3\cdot 2^{2r}-9\cdot 2^r+15}{2^{2r+1}+1}>1.$$
Consequently, if $r\geq 3$ is chosen such that $2^r-1$ is a Mersenne prime, the Frobenius group $\overline{G}$ is an example of a finite group for which inequality (\ref{r1}) does not hold. Finally, to obtain infinitely many such examples, take $G=\overline{G}\times C_q$, where $q$ is any odd prime such that $q\nmid 2^r-1$. Then, $H=\overline{H}\times C_q$ is a subgroup of $G$. By Lemma 2.1 \textit{iv)}, we deduce that
$$\psi'_H(G)=\psi'_{\overline{H}}(\overline{G})\psi'_{C_q}(C_q)=\psi'_{\overline{H}}(\overline{G})>1,$$
and this concludes the proof.
\hfill\rule{1,5mm}{1,5mm}\\

Using the notations from our previous proof, note that $\psi'_{\overline{H}}(\overline{G})=\frac{3\cdot 2^{2r}-9\cdot 2^r+15}{2^{2r+1}+1}$ is an increasing function with respect to $r$ and $\displaystyle \lim_{r \to\infty}\psi'_{\overline{H}}(\overline{G})=\frac{3}{2}$. Hence, the smallest of our examples of finite groups for which inequality (\ref{r1}) does not hold is obtained for $p=2$ and $r=3$. In this case, we have $\overline{G}\cong C_2^3\rtimes C_7$ (SmallGroup(56,11)) and $\psi'_{\overline{H}}(\overline{G})=\frac{45}{43}.$ Also, note that all the groups that were outlined in our previous proof are solvable so, in general, Theorem 1.2 does not hold for this class of finite groups.

An unsolved problem in group theory is the following one: \textit{Given a group $G$ of order $n$, does there necessarily exist a bijection $f$ from $G$ onto $C_n$ such that for each element $x\in G$, we have $o(x)|o(f(x))?$} (see Problem 18.1 posed by I. M. Isaacs in \cite{6}). The answer is true if $G$ is a solvable group. We can formulate a similar question concerning element orders relative to subgroups: \textit{Let $G$ be a group of order $n$, $H$ be a subgroup of order $m$ of $G$ and $H_m$ be the unique subgroup of order $m$ of $C_n$. Does there necessarily exist a bijection $f$ from $G$ onto $C_n$ such that for each element $x\in G$, we have $o_H(x)|o_{H_m}(f(x))?$} As a consequence of the proof of Theorem 1.3, we can state that the answer is negative even for solvable groups.

In \cite{9}, it was proved that inequality (\ref{r1}) holds whenever $H$ is a normal subgroup of a finite group $G$. In this paper, we show that this property also occurs for any subgroup $H$ of $G$ of prime index.\\

\textbf{Proposition 2.2.} \textit{Let $G$ be a group of order $n=mq$ and $H$ be a subgroup of $G$ of prime index $q$. Then $\psi_{H}(G)\leq \psi_{H_m}(C_n)$, where $H_m$ is the unique subgroup of order $m$ of $C_n$.}

\textbf{Proof.} By Lemma 2.1 \textit{i)}, \textit{iii)}, we have
$$\psi_{H_m}(C_n)=|H_m|\psi(C_q)=m(q^2-q+1).$$
Then, the result easily follows by Lemma 2.1 \textit{vi)}.
\hfill\rule{1,5mm}{1,5mm}\\

Finally, we obtain some upper bounds for $\psi'_H(G)$ in terms of the prime divisors of the index of $H$ in $G$.\\

\textbf{Proposition 2.3.} \textit{Let $G$ be a finite group and $H$ be a subgroup of $G$ of index $q=p_1^{\alpha_1}p_2^{\alpha_2}\ldots p_k^{\alpha_k}$, where $p_1<p_2<\ldots<p_k$ are primes and $\alpha_i$ is a positive integer, for all $i\in\lbrace 1, 2,\ldots, k\rbrace$. Then,
\begin{itemize}
\item[i)] $\psi'_H(G)<\prod\limits_{i=1}^k\frac{p_i+1}{p_i}\leq (\frac{3}{2})^k$;
\item[ii)] $\psi'_H(G)<\frac{p_k+1}{p_k}$; in particular, if $k=1$, then $\psi'_H(G)<\frac{p_1+1}{p_1}\leq\frac{3}{2}.$
\end{itemize}}
\textbf{Proof.} By Lemma 2.1 \textit{iii)}, \textit{vi)}, we get
$$\psi'_H(G)=\frac{\psi_H(G)}{\psi_{C_{|H|}}(C_{|G|})}\leq \frac{q^2-q+1}{\psi(C_q)}.$$
Consequently, for item \textit{i)}, we apply Lemma 2.1 \textit{i)} to obtain
$$\psi'_H(G)\leq(q^2-q+1)\prod\limits_{i=1}^k\frac{p_i+1}{p_i^{2\alpha_i+1}+1}<\frac{q^2-q+1}{q^2}\prod\limits_{i=1}^k\frac{p_i+1}{p_i}<\prod\limits_{i=1}^k\frac{p_i+1}{p_i}\leq \bigg (\frac{3}{2}\bigg )^k,$$
while, for item \textit{ii)}, we use Lemma 2.1 \textit{ii)} to get
$$\psi'_H(G)\leq \frac{q^2-q+1}{q^2}\cdot \frac{p_k+1}{p_1}<\frac{p_k+1}{p_1}.$$
\hfill\rule{1,5mm}{1,5mm}\\

In what concerns the last result and its proof, note that it is not always true that $f(q)=\frac{q^2-q+1}{\psi(C_q)}<\frac{3}{2}$. For instance, if we take $q=2^a\cdot 3$, where $a\geq 1$ is an integer, we have
$$\displaystyle \lim_{a \to\infty}f(q)=\displaystyle \lim_{a \to\infty}\frac{9\cdot 2^{2a}-3\cdot 2^a+1}{\frac{7}{3}(2^{2a+1}+1)}=\frac{27}{14}>\frac{3}{2}.$$

\bigskip\noindent {\bf Acknowledgements.} This work was supported by a grant of the "Alexandru Ioan Cuza" University of Iasi, within the Research Grants program, Grant UAIC, code GI-UAIC-2021-01.


\section{Introduction}

\begin{quote}
    ``You have to simulate in your mind how the program's going to work, and you have to have a complete grasp of how the various pieces of the program work together.'' ---Bill Gates \cite[p. 73]{lammers1986}
\end{quote}

There is a considerable body of research that has studied software engineers as ``problem solvers,'' that seeks to understand how developers make sense of problems and how they solve design puzzles \citep{petre1997}.
The sense of accomplishment from solving problems and programming activities is often seen as a critical characteristic of developers \citep{almstrum2003attraction}.
Similarly, a review article on motivation factors of software engineers concludes that:
\begin{quote}
``learning, exploring new techniques and problem solving appear to be the motivating aspects of SE. However little work has focused on the specific nature of Software Engineering itself [...] neither did any of the models focus on the nature of the SE's job itself such as [...] the logical nature of problem solving, use of creativity, complex problem solving, and so on.''
\cite[p. 874]{beecham2008motivation}
\end{quote}

Beecham et al. suggest that a software engineer's job involves problem solving and use of creativity---in other words, cognitive tasks.
The role of cognition and, in particular, the ability to solve problems in programming has been acknowledged as an important factor in developer productivity since the early days of computing \citep{curtis1984fifteen}.
One key psychological attribute for software engineers is the joy of effortful cognitive activities, or what has been termed \textit{need for cognition} (NfC) \citep{cacioppo1982need}. The level of need for cognition indicates the extent to which a developer enjoys cognitive effort, which is directly related to problem solving skills \citep{heppner1983cognitive}.
Understanding what drives developers' need for cognition can be
useful to software development organizations, for example, in hiring decisions and onboarding new staff \citep{sharma2020}.
Despite decades of research and an increasing level of attention for human aspects in software engineering \citep{cruz2015forty,curtis1984fifteen,lenberg2015behavioral,storey2020}, need for cognition of software engineers has, as of yet, not been studied (cf. \cite{capretz2003}).
Given the importance of the role that need for cognition plays in software engineers' problem solving skills, it is surprising that there is a paucity of research that has focused on this construct.
Prior studies that focused on software development tasks and developer preferences have sought explanations in developers' personality traits.
For example, studies have linked  personality traits to developers' ability in code review \citep{da2007does}, exploratory testing \citep{shoaib2009empirical}, object-oriented programming \citep{cegielski2006makes}, and programming proficiency \citep{evans1989best}.
Similarly, other studies highlight the relation between developer behavior and attitudes,  and personality traits. For example, decision-making \citep{feldt2010links}, working attitude \citep{smith1989personality}, or preferences for using offline batch versus online time-sharing equipment \citep{lee1978personality}. More recently, \cite{mellblom2019} found a strong link between neuroticism (one of the personality traits in the Five Factor Model that is also used in this study) and developer burnout.
These studies suggest that personality traits theory can be an insightful theoretical perspective to study software developers.
Hence, in this study we adopt this perspective as we seek to find empirical evidence, rather than relying on anecdote, to establish developers' need for cognition. Specifically, we address the following research question:

\begin{quote}
    \textbf{Research Question}: What is the relationship between personality traits of software engineers and their need for cognition?
\end{quote}

To answer this question, we conducted a large-scale sample study of almost 500 software professionals. We analyzed this rich data set using Bayesian analysis procedures.
Based on the findings, we offer recommendations for software organizations for attracting and retaining software professionals.

This article is structured as follows.
Section~\ref{sec:related} provides an overview of prior studies of personality traits of software developers.
Sec.~\ref{sec:design} presents the research design.
Afterward, in Sec.~\ref{sec:analysis}, we present the results of our analysis.
Findings with recommendations are discussed in Sec.~\ref{sec:discussion}. Finally, Sec.~\ref{sec:conclusion} concludes the paper and offers suggestions for future research.

\section{Background and Related Work}
\label{sec:related}

Software engineering is a socio-technical field of research \citep{perry1994people}, which seeks to advance the state of practice of software development.
While historically the software engineering field has had a strong focus on the technical aspects, there has been a growing interest in studying human factors in software engineering \citep{boehm1988understanding,glass2002facts,lenberg2015behavioral,feldt2010links,storey2020}.

To explain developer preferences and behavior, several studies have considered personality traits \citep{weinberg1971psychology,lee1978personality,shneiderman1980software}.
For example, \cite{li2020distinguishes} provides ``18 attributes pertained to engineers’ personalities'' by asking developers to rank 54 attributes.
\cite{cruz2015forty} identified around 90 articles on personality research in software engineering, published between 1970 and 2010.
However, the interest for personality studies has surged dramatically in the last 20 years or so.  Over 70\% of studies reviewed by \cite{cruz2015forty} were published after 2002.

Need for cognition \citep{cacioppo1982need} is a critical psychological attribute for software professionals, and  related to problem-solving skills \citep{heppner1983cognitive},
which are essential for productive software engineers \citep{graziotin2014happy}.
We identified only three studies that discussed the need for cognition related to computer applications, though none of them studied software professionals specifically.
The most recent study concludes through a small-scale laboratory experiment involving 18 undergraduate students that using a web development learning platform improves students' need for cognition \citep{bergande2020codetripping}.
However, since several previous studies in social psychology  through investigations involving a large number of subjects suggest that need for cognition is a stable trait \citep{furnham2013stable,roberts2008stable,bruinsma2018longitudinal}, we suggest treating the findings of that laboratory experiment with care.
Another study by \cite{amichai2007effects} tested several hypotheses related to the level of need for cognition and internet use.
Findings from a laboratory experiment with 182 participants suggest that people use the same number of hyperlinks regardless of their level of need for cognition; that individuals scoring high in NfC are less influenced by interactive or aesthetic components; and that the time spent on the internet is the same for both groups.
Finally, \cite{van2006paradox} conducted a laboratory experiment involving 43 participants and demonstrated that the usability of Graphical User Interfaces does not depend on the users' NfC.
We identified only one study that linked NfC to personality traits \citep{sadowski1997need}, but like the studies discussed above, this study also relied on a sample of participants from the general population rather than software engineers.
In this study, we seek to investigate the relationship between personality traits and need for cognition that specifically focuses on software developers.
Studying software developers  will allow us to draw comparisons with studies using samples of a general population, and to identify specific implications for a software engineering context.

\subsection{Personality Traits}

Personality can be considered as a set of patterns of thinking, feeling, and behaving, which are based on a set of traits that are predictors of a person's behavior and action \citep{corr2009personality}.
Karl Gustav Jung first introduced personality traits-based theories in the early 1930s \citep{jung1931basic}.
These theories were later operationalized employing measurement instruments, such as the Meyers-Briggs type indicators (MBTI) \citep{myers1976introduction}.
At least until 2010, the MBTI was the most widely used personality test within software engineering research \citep{cruz2015forty}.

However, the Myers-Briggs classification has been criticized for its lack of validity and utility, such as unstable test-retest reliability and inaccurate predictive validity \citep{boyle1995myers}.
Moreover, the MBTI is not correlated with other personality scales \citep{furnham1996big}, and has not been validated by empirical research \citep{druckman1991mind}.
Psychology researchers observed that the Meyers-Briggs type indicators do not measure categorically distinct types, generating quasi-random traits assignments~\citep{mccrae1989reinterpreting}.
The reason is that the MBTI is based on a dichotomous preference score and not a continuous one, which limits the extent to which personality traits can be measured~\citep{devito1985review}.

To address these limitations, McCrae and colleagues developed the Five-Factor Model (FFM), which is grounded in a series of empirical evaluations \citep{mccrae1989reinterpreting,mccrae1992introduction}.
The FFM is also known as the OCEAN model, an acronym for the five  personality traits it measures: Openness (O), Conscientiousness (C), Extraversion (E), Agreeableness (A), and Neuroticism\footnote{The term `Neuroticism' should not be used as such since it might suggest a mental health issue. Emotionality or Emotional stability are better terms.} (N).
\cite{ashton2004six} have extended the OCEAN model by including an additional trait: Honesty-Humility.
This trait emerged from large-scale studies that identified it as a distinct personality trait \citep{saucier2009recurrent}.
The HEXACO model has become the new standard in social psychology \citep{ashton2014hexaco}.
Replication studies of the OCEAN model (which used more advanced computation techniques that were not available when the OCEAN model was initially proposed) confirmed the validity of this new trait \citep{ashton2013individual}.

The characteristics (or \textit{facets}) of the six personality traits of the HEXACO model are:

\begin{enumerate}
\item Honesty-Humility (H-H): sincerity, fairness, greed avoidance, and modesty.
\item Emotionality (Emo): fearfulness, anxiety, dependence, and sentimentality.
\item Extraversion (Ext): social self-esteem, social boldness, sociability, and liveliness.
\item Agreeableness (Agr): forgivingness, gentleness, flexibility, and patience.
\item Conscientiousness (Con): organization, diligence, perfectionism, and prudence.
\item Openness to Experience (OtE): aesthetic appreciation, inquisitiveness, creativity, and unconventionality.
\end{enumerate}

Personality traits can be classified as either `bright' or `dark' \citep{judge2009bright}.
The HEXACO (as also the OCEAN) model focuses on bright traits, i.e. those traits that are considered ``desirable.''
High scores (low scores only for Emotional stability, since it measures the tendency to experience danger) indicate bright personality traits.

There are also dark personality traits; that is, high scores on these traits suggest less desirable characteristics.
Commonly used in personality studies is the so-called ``Dark Triad'' comprising
Narcissism, Psychopathy, and Machiavellianism, which are defined as follows \citep{paulhus2002dark}:

\begin{enumerate}
\item Narcissism (Nar): is related to the individual gratification through self-centered, selfish, and egoistic behaviors.
\item Psychopathy (Psy): anti-social behaviors such as lack of empathy or remorse characterize this trait.
\item Machiavellianism (Mac): is named after Niccolo Machiavelli (author of \textit{The Prince} \citep{machiavelli1532}), focusing on the manipulation to gain self-interested goals and readiness to exploit other people to achieve them.
\end{enumerate}

\subsection{Need for Cognition}

Need for cognition is a psychological attribute advanced in the early 1980s by \cite{cacioppo1982need}, and refers to an individual's tendency to engage and enjoy effortful thinking.
This construct was developed to assess differences among individuals of cognitive processing \citep{sadowski1997need}.
Individuals with high scores on this construct have positive attitudes toward problem-solving, reasoning, and abstraction \citep{cacioppo1996dispositional}.
People scoring low rely on cognitive heuristics and stereotypes to interpret reality \citep{leary2009handbook}.

Individual psychological attributes are measurable proxies of behavioral characteristics \citep{salzberger2013attempting}, such as intelligence and satisfaction with life \citep{kosinski2013private}.
Personality traits are also defined as a set of psychological attributes, which are able to predict other attributes, for example cognitive ability  \citep{roberts2007power}.
In this study we assess the link between personality traits and need for cognition.

It is important to emphasize that need for cognition does not represent a person's intelligence, although some studies found a positive correlation between these two constructs \citep{hill2013need}. A key difference is that need for cognition  only reflects information processing motivations, and not cognitive abilities \citep{cacioppo1996dispositional}.
Notwithstanding, need for cognition is related to other attributes, such as cognitive and sensory innovativeness, which refers to the joy to engage in new experiences to ``stimulate the mind'' \citep{venkatraman1990differentiating}, attributional complexity (i.e., the level of interest or motivation about reality, a preference for complex over simple explanations, meta-cognition capabilities, abstraction skills, and a tendency to infer abstract or complex attributions) \citep{fletcher1986attributional}, and a high level of tolerance for ambiguity and cognitive flexibility \citep{petty1996individual}.

\cite{tuten2001understanding} noted that need for cognition is usually correlated with certain personality traits, in particular with openness to experience (defined above).
However, scholars have pointed out that need for cognition is a theoretically distinct construct: an individual might be open to new experiences \textit{without} reflecting on the consequences of such new experiences \citep{kearney2009and}.
On the other hand, introverts who are reluctant to new experiences might enjoy thinking in a reassuring environment, such as reading an interesting book in the comfort of their home.
In fact, extraversion (the tendency to enjoy human interactions) is not related to need for cognition \citep{sadowski1997need}.
Therefore, regardless of the fact whether personality traits are correlated or not with need of cognition, all those personality attributes have unique meanings and are thus, considered separately \citep{mccrae1996social}.

Need for cognition is also a relevant predictor of positive behaviors in working environments.
Research suggests that a high level of need for cognition is a robust negative moderator of social loafing (the tendency to exert less effort when working in groups), showing significantly higher team efforts \citep{smith2001individual}.
Further, people with high need for cognition values tend to exhibit higher job satisfaction if they face complex engineering tasks, and are generally more satisfied with their jobs than employees scoring low \citep{park2008need}.
Finally, the need for cognition reduces depressed moods at work \citep{gallagher2012managing}.

Anecdotal evidence suggests that software engineers would score high on need for cognition, but there are no studies that have empirically investigated this. Given the socio-technical nature of software development and the importance of human aspects in software engineering, we set out a large-scale study to investigate need for cognition among software engineers.

\section{Research Design}
\label{sec:design}
To investigate the need for cognition attribute of software developers, we conducted a large-scale sample study.
Sample survey research is highly suitable to seek answers that can be generalized to a larger population \citep{stol2018abc}.
In this study we seek to understand the relationship between personality traits and need for cognition of software engineers, considering both bright and dark personality traits.

\subsection{Measurement Instrument}

All constructs in the study were defined as latent variables which we measured using existing instruments that we adopted from the psychology literature.
To measure the bright personality traits, we adopted the HEXACO model of personality \citep{ashton2004six}.
While different survey instruments exist for measuring traits, we adopted the Brief HEXACO Inventory (BHI), which consists of 24 questions \citep{de2013BHI}.
This instrument has shown a high level of stability and accuracy; moreover, it correlates well with other similar and longer instruments (e.g., the HEXACO-PI-R) \citep{de2013BHI}.
To measure the dark traits, we adopted the so-called ``Dirty Dozen'' \citep{jonason2010dirty}, which is also a widely used survey instrument in psychology research \citep{rauthmann2012dark}.
Need for cognition was measured using the NCS-6 instrument, which is a parsimonious instrument with six items \citep{lins2018very}.
We followed the instrument guidelines (which indicate which items are reverse-coded).
We calculated factor scores (which represent a respondent's score on a particular trait) using the mean function; which is a standard procedure when using frequently used instruments such as this \citep{weiner2017handbook}.

\subsection{Data Collection and Screening}

To collect data we used the Prolific data collection platform, which is a dedicated platform for academic research (www.prolific.co) with over 88,000 active members worldwide at the time of the survey.
The use of Prolific for data collection has several advantages over other data collection procedures (for example, invitations through mailing lists), such as reliability, replicability, and data quality \citep{peer2017beyond,palan2018prolific}.
Prolific is widely used in other research fields such as economics \citep{marreiros2017now}, psychology \citep{callan2017interrelations}, and food science \citep{simmonds2018show}.
More recently, studies within computer science and software engineering have adopted the Prolific as a data collection platform \citep{Hosio2020CrowdsourcingDiets,van2021effect,russo2020predictors,russo2021agile_success_model,russo2021daily,cucolas2021impact,russo2021understanding}.
To implement the survey questionnaire, we used Qualtrics (www.qualtrics.com).
To minimize response bias, we randomized questions within their blocks \citep{paulhus1991,gravetter2018research}.

The choice of sampling strategy is key for sample research.
\cite{baltes2020sampling} recently suggested that most studies in software engineering research use non-probability sampling, such as purposive or convenience sampling, and less than 10\% use random sampling.
We selected a cluster-randomized sampling strategy \citep{gravetter2018research}: a number of potential respondents were selected at random from the Prolific population.
In other words, we did not selected the entire world population, but just a cluster of it (i.e., the Prolific community).
From our cluster we proceeded to identify software professionals through a multi-stage selection process described below.
Compared to random sampling, this is less precise (since \textit{only} a portion of the population is considered, in this case, that portion which has a presence on Prolific), but it is very cost-effective compared to the other strategies.
To be sure to include only software developers in our sample, we ran a three-step screening process as detailed below.

\paragraph{Step 1. Pre-Screening.} The first step was a pre-screening of the active members on Prolific based on some general criteria related to software development, such as knowledge of software development techniques, doing computer programming for a living, and the use technology at work. Further, we required a Prolific platform approval rate of 100\%, suggesting that all included respondents successfully completed prior surveys (and did not fail any screening criteria in previous surveys).
After pre-screening, our sample included 2,897 people.

\paragraph{Step 2. Competence Screening.} In the second step we invited candidates who self-identified as software engineers to take part in a competence screening to ensure that they were sufficiently knowledgeable in the subject matter.
We asked three multiple-choice questions,\footnote{Questions are included in the supplementary material.} one about software design and two about programming, and also time-boxed the responses within three minutes to avoid suspicious behavior and cheating.
At the end of this stage, we had a sample of 514 participants.

\paragraph{Step 3. Quality Screening.} At this point we administered the survey to the remaining participants, but included a quality screening to ensure high-quality data. The quality screening comprised three attention checks; participants who did not pass these checks were excluded from the sample.
As a final result of the entire screening process, we received 483 valid and complete responses that could be used for analysis.

Ethical issues were considered according to the recommendations in the Declaration of Helsinki \citep{general2014Helsinki}.
Participants were identified only through a user ID provided by the Prolific data collection platform.
We did not collect any potentially sensitive information of participants.
Further, participants were free to withdraw their participation at any time up to the point of submitting, and some (\textit{n}=306) decided not to participate to this investigation although invited.
Finally, several of the authors have had research ethics and integrity training, which was specifically tailored to behavioral sciences.

\subsection{Sample Description}\label{sub:sample}
The final sample included 483 complete responses, representing the largest sample for a personality research study in software engineering \citep{cruz2015forty}.
In addition to data on personality traits, we also collected demographic data.
Respondents were between 18 and 68 years of age, with an average of 33 years (see Fig.~\ref{fig:age}).
Table~\ref{tab:gender} shows that almost 19\% of  respondents were women.
Most of the respondents (ca. 93\%) were born in Western countries.
Table~\ref{tab:Country} lists the 15 most frequent countries of origin (representing approximately 90\% of all subjects).
Almost three-quarters of the selected professionals have at least one university degree (Table~\ref{tab:Education}).

\begin{figure}[!ht]

\includegraphics[width=2.5in,trim={0 2cm 0 0cm},clip]{agebplot.pdf}
\caption{Distribution of participants' age}
\label{fig:age}
\end{figure}

\begin{table}[!t]
\sisetup{
    group-digits=true,
    group-minimum-digits=4,
    mode=text,
    detect-weight=true,
    detect-family=true
}
\caption{Gender distribution of the sample}
\label{tab:gender}
\begin{tabular}{m{3cm}SS}
\toprule
Gender & {Frequency} & {Percent} \\ \midrule
Men & 393 & 81.4 \\
Women & 90 & 18.6 \\
Total & 483 & 100 \\ \bottomrule
\end{tabular}
\end{table}

\begin{table}[!ht]
\sisetup{
    group-digits=true,
    group-minimum-digits=4,
    mode=text,
    detect-weight=true,
    detect-family=true
}
\caption{Respondents' country of origin}
\label{tab:Country}
\begin{tabular}{m{4cm}SS}
\toprule
Country & {Frequency} & {Percent} \\ \midrule
United Kingdom & 141 & 29.2 \\
USA & 135 & 28.0 \\
Portugal & 33 & 6.8 \\
Poland & 22 & 4.6 \\
Italy & 18 & 3.7 \\
Canada & 15 & 3.1 \\
Germany & 12 & 2.5 \\
Spain & 9 & 1.9 \\
Ireland & 9 & 1.9 \\
Greece & 8 & 1.7 \\
Mexico & 8 & 1.7 \\
Australia & 7 & 1.4 \\
France & 6 & 1.2 \\
Hungary & 5 & 1.0 \\
Estonia & 4 & 0.8 \\
Other & 51 & 10.5 \\ \bottomrule
\end{tabular}
\end{table}

\begin{table}[!ht]
\sisetup{
    group-digits=true,
    group-minimum-digits=4,
    mode=text,
    detect-weight=true,
    detect-family=true
}
\caption{Highest degree of education of respondents}
\label{tab:Education}
\begin{tabular}{m{4cm}SS}
\toprule
Education & {Frequency} & {Percent} \\ \midrule
Bachelor's degree & 241 & 49.9 \\
Master's degree & 105 & 21.7 \\
Some college but no degree & 77 & 15.9 \\
High school graduate & 37 & 7.7 \\
Doctoral degree & 17 & 3.5 \\
Less than high school degree & 2 & 0.4 \\
Other & 4 & 0.8 \\ \bottomrule
\end{tabular}
\end{table}

\begin{table}[!ht]

\sisetup{
    group-digits=true,
    group-minimum-digits=4,
    mode=text,
    detect-weight=true,
    detect-family=true
}
\caption{Respondents' roles within their organizations}
\label{tab:Role}
\begin{tabular}{m{5.4cm}SS}
\toprule
Role &  {Freq.} & {Percent} \\ \midrule
Software developer, programmer & 252 & 52.2 \\
Data analyst, engineer, scientist & 44 & 9.1 \\
Technical support & 32 & 6.6 \\
Team Lead & 31 & 6.4 \\
DevOps engineer, infrastructure developer & 22 & 4.6 \\
Product manager & 21 & 4.3 \\
Tester, QA engineer & 16 & 3.3 \\
Architect & 12 & 2.5 \\
CIO, CEO, CTO & 12 & 2.5 \\
Systems analyst & 12 & 2.5 \\
UX, UI designer & 9 & 1.9 \\
Other & 20 & 4.1 \\ \bottomrule
\end{tabular}
\end{table}

\section{Bayesian Data Analysis}
\label{sec:analysis}

We conducted the analyses using several software packages. After collecting the data, we performed several sanity checks using IBM SPSS (version 26).
The primary analyses were conducted using the JASP (version 0.11.1) statistical package \citep{jasp2018jasp}. JASP is a relatively new open source package that is actively being developed, and seeks to offer a user-friendly interface to conduct Bayesian analyses.
We modeled personality traits as exogenous variables (or predictors) and need for cognition as the endogenous (or criterion) variable.

The interest in Bayesian statistics is rising in the empirical software engineering research community \citep{furia2017,furia2019bayesian,furia2021}.
One reason for this is an increasing awareness of the limitations of `traditional' frequentist statistics, which relies on \textit{p}-values and significance testing of null-hypotheses \citep{cohen1994earth,baker2016statisticians}.
By convention, \textit{p}-values are considered as a measure of evidence in favor (or against) of a null hypothesis. Small \textit{p}-values indicate strong evidence against the null hypothesis, and \textit{p}-values lower than 0.05 are widely considered as the critical threshold that determines when we have sufficiently strong evidence against the null hypothesis. Nevertheless, \textit{p}-values cannot be interpreted as the probability that the null hypothesis is true. What they provide is the probability of having data at least as extreme as the one observed, if the null hypothesis was indeed true.\footnote{Also, the threshold lower than 0.05 for the \textit{p}-values is a convention that has to be corrected based on the number of predictors to avoid Type I errors \citep{ACMStd,russo2020predictors}.} Consequently, the use of \textit{p}-values as a measure of strength of the evidence to support any claims must be done with caution \citep{nuzzo2014scientific}. In recent years, scholars of statistics have stressed these arguments and advocated for the use of Bayesian approaches instead of frequentist ones, since the former can offer a more comprehensible, nuanced understanding of the statistical data analysis \citep{wagenmakers2018bayesian,dienes2014Bayes,lee2018determining,furia2019bayesian,JASP2019guidelines}.

Bayesian statistics are not based on a dichotomous evaluation mode (i.e., hypotheses are not accepted or rejected). Instead, the Bayesian  paradigm accepts that different explanations may be possible and that evidence favoring one or another hypothesis can be properly quantified. Bayesian methods allow the computing of the probability that a given hypothesis is true.\footnote{When the set of considered hypothesis does not contain a hypothesis that perfectly represents the underlying data generation process, this probability should be interpreted as the probability that a given hypothesis is the best among all considered hypotheses \citep{masegosa2020learning}.} This probability is the so-called \textit{Bayesian posterior probability}, which is denoted as P(H\textpipe{}data). This probability  quantifies the strength of our evidence in favor of \textit{H}. Bayesian posterior probabilities are computed by applying the Bayes rule over the following two quantities: the observed data's prior probability and the likelihood. Prior probabilities allow us to integrate any prior evidence that we have in favor of a given hypothesis. However, non-informative (usually uniform) probabilities can also be used when little or no prior information is available. The likelihood of the data can be easily understood in terms of Bayes factors. A Bayes factor, defined for two competing hypotheses and noted as BF\textsubscript{10}, represents the observed evidence in favor of \textit{H}\textsubscript{1} over \textit{H}\textsubscript{0}.\footnote{An alternative notation is BF\textsubscript{01}, representing the evidence in favor of \textit{H}\textsubscript{0} over \textit{H}\textsubscript{1}.} Table~\ref{tab:bayesian} presents heuristics to interpret Bayes factors. Rather than delineating explicit thresholds or cut-off values, beyond which a conclusion is supported, as is typical in frequentist statistics, Bayes factors should be interpreted more qualitatively. For example, an analysis resulting in a Bayes factor BF\textsubscript{10} of, say, 5, suggests moderate evidence for \textit{H}\textsubscript{1} over \textit{H}\textsubscript{0}. If the Bayes factor is 50, on the other hand, the heuristics suggest that there is robust evidence for \textit{H}\textsubscript{1} over \textit{H}\textsubscript{0}. Consequently, the `thresholds' in Table~\ref{tab:bayesian} should not be interpreted strictly but rather  heuristically, as a continuum representing strength of evidence. Finally, the following equation described how Bayesian posterior probabilities, prior probabilities, and Bayes factors relate:

\[\frac{P(H_1|D)}{P(H_0|D)} = \frac{P(H_1)}{P(H_0)} \cdot BF_{10} \]

\noindent For example, if \textit{H}\textsubscript{1} and \textit{H}\textsubscript{0} are, a priori, equally probable (i.e., $\frac{P(H_1)}{P(H_0)}=1$), but the observed data suggests strong evidence for \textit{H}\textsubscript{1} with a BF\textsubscript{10} equal to 20, then we have that $\frac{P(H_1|D)}{P(H_0|D)} = 20$. And, solving the equation, we have that \textit{P(H\textsubscript{1}|D)=$\frac{20}{1+20}\approx0.95$} and \textit{P(H\textsubscript{0}|D)=$\frac{1}{1+20}\approx0.05$}.

One of the critical advantages of Bayesian statistics is that the above reasoning can be naturally extended to multiple hypotheses or to more complex settings like multi-model linear regression (as we do in the next section). Because, in Bayesian statistics, all models and observations are equally treated as random variables. So, all computations are always done following standard probability operations, like the ones described above.

\begin{table}[!t]
\sisetup{
    group-digits=true,
    group-minimum-digits=4,
    table-format=0.3,
    mode=text,
    detect-weight=true,
    detect-family=true
}
\caption{Heuristics for interpretation of Bayes factors BF$_{10}$ \cite[p. 105]{lee2013}}
\label{tab:bayesian}
\begin{tabular}{p{3cm}p{5cm}}
\toprule
Bayes factor & Evidence category \\
\midrule
$>$ 100 & Extreme evidence for H$_1$\\
30 -- 100& Very strong evidence for H$_1$\\
10 -- 30& Strong evidence for H$_1$\\
3 -- 10& Moderate evidence for H$_1$\\
1 -- 3& Anecdotal evidence for H$_1$\\
1& No evidence\\
1/3 -- 1& Anecdotal evidence for H$_0$\\
1/10 -- 1/3& Moderate evidence for H$_0$\\
1/30 -- 1/10& Strong evidence for H$_0$\\
1/100 -- 1/30  & Very strong evidence for H$_0$\\
$<$ 1/100& Extreme evidence for H$_0$\\
\bottomrule
\end{tabular}
\end{table}

\subsection{Bayesian multi-model linear regression}\label{sec:BMMLR}

In this study we use a Bayesian multi-model linear regression. Under this model, we assume that the nine personality traits studied here (a.k.a. the independent variables) have a linear relationship with need for cognition (a.k.a. the dependent variable). But, crucially, this analysis assumes that not necessarily all the personality traits are needed to predict developers' need for cognition. In consequence, we explore all regression models defined by all possible combinations of the different personality traits. As we are considering nine traits, this analysis explores 2\textsuperscript{9}=512 alternative regressions models. Each model is treated as a different hypothesis.

We followed the standard guidelines reported by \cite{JASP2019guidelines} and \cite{van2021tutorial} for  conducting the Bayesian data analysis. We first performed several checks to verify that the use of a linear regression model to study the relationship between the considered personality traits and need for cognition is a reasonable choice. Following \cite{van2021tutorial}, we assessed the assumptions of normality, linearity, homoscedasticity, and the absence of multicollinearity. We discuss the procedures for assessing these next.

\textbf{Data Normality.} We first evaluated data normality.
Residuals can be seen as the differences between the observed value of need for cognition (i.e., the endogenous variable) and the personality trait values (i.e., the exogenous variables).
Figure~\ref{fig:PPplot} presents a normal predicted probability plot (P-P plot), which shows that the residuals are normally distributed since they are reasonably aligned to the diagonal normality line.

\begin{figure}[!b]
\includegraphics[trim={4cm 0cm 0cm 0cm},clip,width=3.5in]{PPplot.png}
\caption{Normal P-P plot of regression standardized residuals for need for cognition.}
\label{fig:PPplot}
\end{figure}

\textbf{Homoscedasticity.} Second, we evaluated homoscedasticity by checking whether the residuals are equally distributed.
Heteroscedasticity means that residuals are unevenly spread, leading to clustering on some values. This is undesirable as it would mean that the variance of the observations is not equal.
The scatterplot shown in Figure~\ref{fig:ResidualsScatterplot} shows that the residuals are randomly distributed.

\begin{figure}[!t]
\includegraphics[width=4in]{ResidualsScatterplot.png}
\caption{Scatterplot of the residuals to assess homoscedasticity for need for cognition}
\label{fig:ResidualsScatterplot}
\end{figure}

Since our residuals are normally distributed and homoscedastic, we can also assume that the exogenous variables (i.e., personality traits) in the regression have a straight-line relationship with the endogenous variable (i.e., need for cognition), concluding that the data have a linear relationship.
It is linear because the residuals appear to be linear since we do not see any curved pattern in Figure~\ref{fig:ResidualsScatterplot}.

\textbf{Multicollinearity.} Finally, we checked for multicollinearity; that is, we evaluated whether the exogenous variables are correlated with each other.
This step is essential because our primary goal is to isolate the relationship between the endogenous variable and each of the exogenous variables.
If the exogenous variables are correlated, changes in one of them will also influence another, leading to unreliable results.
A standard measure for assessing multicollinearity is the variance inflation factor (VIF). A commonly accepted cut-off for the VIF is 5, but a more conservative cut-off is 3 \citep{Miles2014VIF,Obrien2007VIF}.
Table~\ref{tab:VIF} shows that all VIF values are below two, well below the cut-off value.
Hence, we are confident that our regression model will associate the variance of the endogenous variable with the exogenous ones.

\begin{table}[!t]
	\sisetup{
    group-digits=true,
    group-minimum-digits=4,
    table-format=0.3,
    mode=text,
    detect-weight=true,
    detect-family=true
}
	\caption{Coefficients for need for cognition to assess multicollinearity}
	\label{tab:VIF}
	{
		\begin{tabular}{p{4.5cm}SSSSS}
			\toprule
			Variables & {Unstandardized}	 & {Standardized} & {\textit{t}} & {\textit{p}} & {VIF}  \\
			\midrule
			Machiavellianism & 0.058 & 0.069 & 1.422 & 0.156 & 1.662  \\
			Psychopathy & 0.073 & 0.079 & 1.608 & 0.109 & 1.718  \\
			Narcissism & 0.012 & 0.014 & 0.318 & 0.751 & 1.400  \\
			Honesty-Humility & 0.258 & 0.260 & 5.226 & 0.000 & 1.747  \\
			Emotionality & -0.129 & -0.128 & -3.088 & 0.002 & 1.210  \\
			Extraversion & 0.060 & 0.060 & 1.384 & 0.167 & 1.316  \\
		    Agreeableness  & -0.097 & -0.082 & -1.958 & 0.051 & 1.235  \\
			Conscientiousness  & 0.237 & 0.217 & 5.268 & 0.000 & 1.202\\
		    Openness to Experiences & 0.426 & 0.360 & 8.976 & 0.000 & 1.135  \\
			\bottomrule
		\end{tabular}
	}
\end{table}

Another key step in a Bayesian analysis is the specification of a prior \citep{JASP2019guidelines,van2021tutorial,furia2021}. In our case, the prior of a Bayesian multi-model regression is composed of two parts. A first part refers to the prior over the 512 different linear regression models explored in this analysis. In our case, we decided to use a uniform prior (i.e. all models have the same prior probability). Eventually, a non-uniform prior could be employed encoding some previously available information about the relevance of some of our personality traits for predicting need for cognition. However, we did not find a sensible way to define a precise probability value for all (or some of) the personality traits reflecting the information present in the literature.

The second part of the prior in a Bayesian multi-model regression refers to the prior over the linear coefficients associated to each of the independent variables (i.e. the personality traits). For example, a negative linear coefficient for a given personality trait indicates an negative correlation between this personality trait and need for cognition. Finding a prior for these continuous values is not as easy as in the case of the different models. Following the guidelines by \cite{JASP2019guidelines} and \cite{van2021tutorial}, we considered four different priors which are commonly used in the Bayesian literature for the coefficients of a linear regression model \citep{liang2008mixtures,rouder2012default}.
By using different priors, we can assess the sensitivity of our analysis to the choice of the prior, which is a standard robustness check in Bayesian data analysis \citep{JASP2019guidelines,van2021tutorial}.
In the following section, we  report the average plus/minus the standard deviation of the different computations, such as the Bayes Factors and the posterior probabilities, obtained by repeating the whole data analysis using a different prior at a time. The standard deviations  will help us to identify the sensitivity of the Bayesian analysis to the choice of the prior. Small standard deviation values, in comparison to the mean value,
suggests that the results do not vary much across different priors, and thus that the prior of choice does not have an impact on the results that would raise any concerns as to the sensitivity of the prior of choice. On the other hand, large standard deviation values, in comparison to the mean value, suggest that the choice of prior \textit{does} have an impact on the results, and this should be reason for concern, because the results depend on a specific prior.

The results of the various computations are included in the supplementary material at the end of this article.
In Appendix \ref{app:BLR}, we provide full details about these priors used and how to reproduced the analysis shown in this paper using JASP.

\subsection{Results}

We now present the results of our analysis. We establish (i) which personality traits best predict developers' need for cognition; (ii) how much developers' personality traits contribute in explaining their need for cognition; and (iii) whether or not software engineers have a higher need for cognition compared to a general population.

The result of the ten best regression models out the 512 models considered in this Bayesian analysis are summarized in Table~\ref{tab:modelComparison-NfC}.
The table shows P(H\textpipe{}data)
which is the posterior probability of each model/hypothesis after observing the data, the Bayes factor (BF\textsubscript{10}),
and the explained variance \textit{R}\textsuperscript{2}.
As noted, Table~\ref{tab:bayesian} offers heuristics to interpret Bayes factors. We  note that JASP fixes as the reference model (H\textsubscript{0} following previous notation) the model on top of the list (this is because BF\textsubscript{10} is equal to one in this case). The other models have Bayes factors smaller than one, showing that they are all less supported by the data compared to the reference model. Nevertheless, as can be seen, no model has an overwhelming advantage over the others. This suggests that any of these subsets of personality traits explains NfC similarly well.
Based on the \textit{R}\textsuperscript{2} values, we conclude that personality traits explain approximately one-third of the variation in need for cognition of software engineers in our sample.

As discussed in the previous section, Table~\ref{tab:modelComparison-NfC} reports the average value plus/minus the standard deviation of the different quantities over four analysis, each of them performed with a different prior. As previously discussed, the standard deviations help us to identify the sensitivity of the Bayesian analysis to the choice of the prior. As can be seen, standard deviations values are very small and, in consequence, the main conclusions of this analysis are not affected by the particular choice of the prior.

\begin{table*}[!t]
\caption{Model Comparison --- need for cognition (machiavellianism (Mac), psychopathy (Psy), narcissism (Nar), honesty-humility (H-H), emotionality (Emo), extraversion (Ext), Agreeableness (Agr), conscientiousness (Con), openness to experiences (OtE)). We report the average value $\pm$ the standard deviation of the different Bayesian quantities over four analysis, each of them performed with a different prior, as detailed in Section \ref{sec:BMMLR}.}
\label{tab:modelComparison-NfC}

	\robustify\bfseries
	\sisetup{
    group-digits=true,
    group-minimum-digits=4,
    table-format=0.3,
    mode=text,
    detect-weight=true,
    detect-family=true
}

	{
 		\begin{tabular}{clccc
			\toprule
 			No. & Models & {P(H\textpipe{}data) } & {BF\textsubscript{10}} & {\textit{R}\textsuperscript{2}}  \\
			\cmidrule[0.4pt]{1-5}
1 & Mac + H-H + Emo + Agr + Con + OtE & 0.13$\pm$ 0.00 & 1.00$\pm$ 0.00 & 0.33 \\
2 & H-H + Emo + Agr + Con + OtE & 0.10$\pm$ 0.01 & 0.77$\pm$ 0.08 & 0.32 \\
3 & Psy + H-H + Emo + Con + OtE & 0.09$\pm$ 0.01 & 0.68$\pm$ 0.07 & 0.32 \\
4 & Psy + H-H + Emo + Agr + Con + OtE & 0.08$\pm$ 0.00 & 0.56$\pm$ 0.00 & 0.32 \\
5 & Mac + Psy + H-H + Emo + Con + OtE & 0.07$\pm$ 0.00 & 0.50$\pm$ 0.00 & 0.32 \\
6 & Mac + H-H + Emo + Con + OtE & 0.06$\pm$ 0.01 & 0.45$\pm$ 0.05 & 0.32 \\
7 & Mac + Psy + H-H + Emo + Agr + Con + OtE & 0.05$\pm$ 0.00 & 0.36$\pm$ 0.03 & 0.33 \\
8 & Psy + H-H + Emo + Ext + Agr + Con + OtE & 0.05$\pm$ 0.00 & 0.34$\pm$ 0.03 & 0.33 \\
9 & Mac + H-H + Emo + Ext + Agr + Con + OtE & 0.04$\pm$ 0.00 & 0.28$\pm$ 0.03 & 0.33 \\
10 & Psy + H-H + Emo + Ext + Con + OtE & 0.03$\pm$ 0.00 & 0.25$\pm$ 0.00 & 0.32 \\
			\bottomrule

		\end{tabular}
	}
\end{table*}

If we jointly consider the different models shown in Table~\ref{tab:modelComparison-NfC}, several traits appear in most of them. Table~\ref{tab:posteriorSummariesOfCoefficients}  lists the individual probabilities of each trait to be relevant for predicting need for cognition. The inclusion probability of a given trait is computed by summing up (i.e., marginalizing) the posterior probabilities of all of the models where this trait is included (all the 512 evaluated models are considered here, not only the ten ones shown in Table~\ref{tab:modelComparison-NfC}). Whereas Table~\ref{tab:modelComparison-NfC} considers traits only in \textit{combinations},
Table~\ref{tab:posteriorSummariesOfCoefficients} presents information on the likelihood of individual traits to be needed to predict the developer's need for cognition. Note that, as in the previous table, we report average values plus/minus the standard deviation in order to evaluate the sensitivity of the choice of the prior.

\begin{table}[!ht]
	\sisetup{
    group-digits=true,
    group-minimum-digits=4,
    table-format=0.3,
    mode=text,
    detect-weight=true,
    detect-family=true
}
	\caption{Posterior summaries of coefficients. {P(incl\textpipe{}data)} denotes the posterior inclusion probability. The last two columns denote the 95\% central credible interval (CI) for the value of the associated linear coefficients. We report the average value $\pm$ the standard deviation of the different quantities over four analyses, each of them performed with a different prior, as detailed in Section~\ref{sec:BMMLR}.}
	\label{tab:posteriorSummariesOfCoefficients}
	{

 		\begin{tabular}{lccc}
			\toprule
 			\multicolumn{1}{c}{}  & \multicolumn{1}{c}{} & \multicolumn{2}{c}{95\% CI} \\
 			\cmidrule{3-4}
			Coefficient & {P(incl\textpipe{}data)} & {Lower} & {Upper}  \\
			\cmidrule[0.3pt]{1-4}
Intercept & 1.00$\pm$ 0.0 & 3.92$\pm$ 0.0 & 4.02$\pm$ 0.0 \\
Mac & 0.49$\pm$ 0.0  & -0.00$\pm$ 0.0 & 0.13$\pm$ 0.0 \\
Psy & 0.52$\pm$ 0.0  & -0.00$\pm$ 0.0 & 0.16$\pm$ 0.0 \\
Nar & 0.16$\pm$ 0.0  & -0.02$\pm$ 0.0 & 0.06$\pm$ 0.0 \\
H-H & 1.00$\pm$ 0.0  & 0.13$\pm$ 0.0 & 0.32$\pm$ 0.0 \\
Emo & 0.96$\pm$ 0.0  & -0.22$\pm$ 0.0 & -0.02$\pm$ 0.0 \\
Ext & 0.26$\pm$ 0.0  & -0.01$\pm$ 0.0 & 0.10$\pm$ 0.0 \\
Agr & 0.61$\pm$ 0.0  & -0.18$\pm$ 0.0 & 0.00$\pm$ 0.0 \\
Con & 1.00$\pm$ 0.0  & 0.15$\pm$ 0.0 & 0.32$\pm$ 0.0 \\
OtE & 1.00$\pm$ 0.0  & 0.34$\pm$ 0.0 & 0.52$\pm$ 0.0 \\
            \bottomrule
		\end{tabular}
	}
\end{table}

As a result of this Bayesian analysis, four traits in particular have an extremely high chance of almost 100\% to be needed to predict need for cognition, namely: openness to experience (OtE), conscientiousness (Con), honesty-humility (H-H), and emotionality (Emo). This is reflected in the probability of inclusion in the regression model after observing the data (i.e., P(incl\textpipe{}data)), with values of  100\% (with the only exception being emotionality which has a probability of 96.3\%). Note that in this analysis, we assumed that \textit{a priori} (i.e., before considering the data from the survey), all traits had the same chance of 50\% to be relevant for predicting need for cognition (i.e., P(incl)=0.5).

The next question is whether these personality traits are positively, negatively, or not correlated with the need for cognition. Table~\ref{tab:bayesianPearsonCorrelations} presents the Pearson coefficients that represent the relations between need for cognition and personality traits. Based on the correlation using Pearson's \textit{r}, we can see the sign (negative or positive) of the individual relation but not the strength of the evidence supporting this relationship. To address this issue, we also compute the Bayes factor for a model where no relationship is established between the trait and need for cognition. In this case, the specification of the prior is different from the case of the Bayesian multi-model regression because we now use a simple correlation model. Again, we consider four different priors and report average plus/minus standard deviations for the quantities depending on the prior. In this case, the standard deviations of some of the Bayes factors are, at first sight, relatively high, but we should note that in these cases, Bayes factors always vary within a range of extremely supportive evidence according to the ranges provided in Table~\ref{tab:bayesian}. This effect can be appreciated when looking at the associated posterior probabilities, which can be computed from the Bayes factors, because, in those cases, the associated standard deviations become negligible. Again, we can see that the main conclusions of this analysis are not affected by the particular choice of the prior. In Appendix \ref{app:BLR}, we provide more details about these priors and about how to reproduce this analysis using JASP.

In particular, we can see how dark traits (i.e., psychopathy, machiavellianism, and narcissism) are negatively associated with need for cognition, although those associations are barely significant (i.e., they have negative Pearson's r coefficients but very small Bayes factors). Honesty-humility, extraversion, conscientiousness, and openness to experience are significantly positively associated with need for cognition (i.e., they have positive Pearson's r coefficients and large Bayes factors). Also, emotionality is strongly negatively associated to NfC. However, according to this analysis, we can not infer whether agreeableness is associated to NfC in a significant way.

\begin{table}[!ht]
	\caption{Bayesian correlation analysis between need for cognition and personality traits. We report the average value $\pm$ the standard deviation of the different Bayesian quantities over four analysis, each of them performed with a different prior, as detailed in Appendix \ref{app:BLR}.}
	\label{tab:bayesianPearsonCorrelations}
		\sisetup{
    group-digits=true,
    group-separator={,},
    group-minimum-digits=4,
    table-format=-1.3,
    mode=text,
    detect-weight=true,
    detect-family=true
}
		\begin{tabular}{lScc
			\toprule
			Trait & {Pearson's \textit{r}} & {P(correlation\textpipe{}data)} & {BF\textsubscript{10}}\\
			\midrule
Mac & -0.02 & 0.06$\pm$ 0.02 & 0.06$\pm$ 0.02 \\
Psy & -0.04 & 0.07$\pm$ 0.02 & 0.08$\pm$ 0.03 \\
Nar & -0.07 & 0.17$\pm$ 0.05 & 0.21$\pm$ 0.07 \\
H-H & 0.24 & 1.00$\pm$ 0.00 & 4.3e+04$\pm$ 1.3e+04 \\
Emo & -0.20 & 1.00$\pm$ 0.00 & 1.1e+03$\pm$ 3.3e+02 \\
Ext & 0.20 & 1.00$\pm$ 0.00 & 1.2e+03$\pm$ 3.6e+02 \\
Agr & -0.11 & 0.44$\pm$ 0.08 & 0.81$\pm$ 0.27 \\
Con & 0.35 & 1.00$\pm$ 0.00 & 3.3e+12$\pm$ 8.3e+11 \\
OtE & 0.42 & 1.00$\pm$ 0.00 & 2.2e+19$\pm$ 4.6e+18 \\
 			\bottomrule
 			\end{tabular}

\end{table}

Finally, we assess the level of need for cognition for software engineers through an analysis of variance (ANOVA).
To do so, we consider the descriptive statistics in Table~\ref{tab:NfC_DescriptiveStatistics} and the boxplot in Figure~\ref{fig:NfC}.
To interpret these values, we compare these to findings by  \cite{popoviciu2011students}, who observed  a mean value of 3.97 (on a scale of 1-5), suggesting that software professionals are ``very likely to engage in and enjoy effortful cognitive activities.''

\begin{table}[!ht]
	\sisetup{
    group-digits=true,
    group-minimum-digits=4,
    table-format=0.3,
    mode=text,
    detect-weight=true,
    detect-family=true
}
	\caption{Descriptive statistics of need for cognition}
	\label{tab:NfC_DescriptiveStatistics}
	{
		\begin{tabular}{lS}
			\toprule
			Statistic & {Value} \\
			\cmidrule[0.4pt]{1-2}
			N & 483  \\
        	Minimum value & 1  \\
			Maximum value & 5  \\
			Mean & 3.970  \\
			Standard Error of Mean & 0.032  \\
			Standard Deviation & 0.694  \\
			Skewness & -0.926  \\
			Standard Error of Skewness & 0.111  \\
			Kurtosis & 1.643  \\
			Standard Error of Kurtosis & 0.222  \\
			\bottomrule
		\end{tabular}
	}
\end{table}

\begin{figure}[!ht]
\includegraphics[trim={0cm 21mm 0cm 0cm },clip,width=3in]{NfC.pdf}
\caption{Distribution of software professionals' need for cognition.}
\label{fig:NfC}
\end{figure}

\begin{table}[!ht]
	\sisetup{
    group-digits=true,
    group-minimum-digits=4,
    table-format=0.3,
    mode=text,
    detect-weight=true,
    detect-family=true
}
	\caption{Need for cognition scores in this and other studies}
	\label{tab:NfC_genpop}
	{
		\begin{tabular}{llSS}
			\toprule
			Study & {Sample size} & {Mean} & {SD} \\
			\midrule
			\cite{kearney2009and} & 83 teams (549 professionals) & 3.33 & 0.22  \\
        	\cite{popoviciu2011students} & 30 graduate students & 3.56 & {NA} \\
        	\cite{madrid2016creativity} & 220 professionals & 3.79 & 0.60 \\
	This study  & 483 software engineers & 3.97 & 0.69 \\
			\bottomrule
		\end{tabular}
	}
\end{table}

We now compare those values with previous studies reporting scores of need for cognition, summarized in Table~\ref{tab:NfC_genpop}.
Based on the reported mean and standard deviation (SD) values of need for cognition, as reported by \cite{kearney2009and} and  \cite{madrid2016creativity}, we replicate the prior research data using two normal distributions (we could not replicate \cite{popoviciu2011students}, since the standard deviation was not reported). We subsequently use these distributions to compare prior work to our results.
First, we ran a one-way ANOVA comparing need for cognition between these two studies and our results. This shows a significant effect of the study sample on need for cognition (F (2) = 38.16, \textit{p} \textless{} 0.01).
The ANOVA confirms that there is an overall significant difference among groups; however, it does not provide any information about which group scores significantly higher in need for cognition.
Therefore, we performed a two-way Tukey posthoc comparison between the studies.
This test is a widespread posthoc analysis based on the assumption of normally distributed data.

This analysis shows that need for cognition in our study is significantly higher, compared to \cite{kearney2009and} (\textit{p} \textless{} 0.001) and \cite{madrid2016creativity} (\textit{p} = 0.001).
Figure~\ref{fig:Boxplot} shows the boxplot distribution for need for cognition scores.
We conclude, thus, that software engineers report a higher need for cognition compared to other professionals.

\begin{figure}[!ht]
\includegraphics[width=3.8in]{Boxplot.pdf}
\caption{Boxplot distribution of need for cognition studies}
\label{fig:Boxplot}
\end{figure}

Upon closer inspection of need for cognition,
Table~\ref{tab:NfC_DescriptiveStatistics} shows that the mean value is high, close to four (3.97) on a scale of one to five.
Also, the descriptive statistics of the variables suggest that need for cognition is relatively skewed towards high values.
The boxplot shown in Figure~\ref{fig:NfC} confirms this, where both upper and lower quartiles are close to the median value of four, with very few outliers below the lower whisker.
Accordingly, we conclude that software engineers have a significantly higher need for cognition than other knowledge workers.

\section{Discussion}
\label{sec:discussion}

\subsection{Discussion of Findings}
Based on the results of our analysis presented in the previous section, we can now address our central research question: What is the relationship between the personality traits of software engineers and their need for cognition?

Following increasing attention to human factors in general and personality traits in particular, this paper presents an investigation that links personality traits and the need for cognition among software professionals.
Psychology literature suggests that personality is stable over time; i.e., it does not overall change along the maturation process of people \citep{cobb2012stability,gustavsson1997stability,blonigen2008stability}.
This stability of personality has also been confirmed for software developers \citep{calefato2019large}.
Need for cognition has also been identified as a stable psychological attribute \citep{furnham2013stable,roberts2008stable}.
A longitudinal study confirmed this stability over time\footnote{For the sake of completeness, need for cognition slightly increases over time for young individuals and slightly decreases for older ones (in a statistically non-significant way). However, Bruinsma and Crutzen did not consider confounding factors that might explain these trends, such as dementia which could mostly explain such decrease.} \citep{bruinsma2018longitudinal}.
This evidence allows us to make reliable predictions because all traits included in this study are stable and not time-dependent.
Therefore, we can conclude that software engineers have a higher need for cognition for the general population, which is predicted by their defining characteristics, i.e., personality traits.
Although we did not investigate causality, we found that personality traits are good predictors of developers' high need for cognition.
Moreover, since both bright and dark personality traits and the need for cognition are stable over time, it suggests that the predictive relationship is also stable over time

Anecdotal evidence has suggested that software engineers must enjoy effortful thinking to be good developers \citep{petre1997}.
This investigation aims to move from anecdote to evidence through a large-scale empirical study.
Operationalizing the joy of thinking with the need for cognition (as shown by  \cite{strobel2018predicting}) brings several advantages to the community's body of knowledge because it allows us to infer previous evidence on the need for cognition with tailored recommendations for software professionals.

We now focus our discussion on the regression model and the role of each trait.
From the posteriors of Table~\ref{tab:posteriorSummariesOfCoefficients}, we can see that four traits have a probability of inclusion that is close to 100\%, namely openness to experience, conscientiousness, honesty-humility, and emotionality.
These are the only four traits that have a significant BF\textsubscript{inclusion} (i.e., higher than ten \citep{lee2014bayesian}), and are present in all ten of the most likely regression models listed in Table~\ref{tab:modelComparison-NfC}.
To recall, we reported only the ten best models of the Bayesian regression that showed the greater likelihood to predict our dependent variable (i.e., Need for Cognition).
In particular, based on the correlation signs of Table~\ref{tab:bayesianPearsonCorrelations}, \textit{we conclude that high openness to experience, high conscientiousness, high honesty-humility, and low emotionality predict a high level of need for cognition.}
Similar patterns have also been observed in a previous study \citep{sadowski1997need}.

Prior studies of personality traits in software engineering are limited to only bright personality traits \citep{cruz2015forty}.
This study is among the first to include  dark personality traits.\footnote{In a previous study, we focused on gender differences in the same set of bright and dark personality traits of software engineers \citep{russo2020gender}.}
The inclusion of dark traits offers a more comprehensive and holistic understanding of personality traits of software engineers and, in this study, in particular, how these relate to developers' need for cognition.
However, Table~\ref{tab:bayesianPearsonCorrelations} shows that none of the three dark traits, machiavellianism, psychopathy, and narcissism, are significantly correlated to the need for cognition (with a BF\textsubscript{10} close to 0).
We therefore conclude that bright traits play a much greater role in predicting the need for cognition than dark traits.
Furthermore, the BF\textsubscript{inclusion} values of these dark traits are also relatively low.
This suggests that dark traits do not play a significant role in the need for the cognition of software professionals.

Additionally, our ANOVA analysis confirmed that \textit{software professionals score significantly higher in their need for cognition compared to the general population.}
Need for cognition seems to be one unique characteristic of developers.

Finally, we also highlight that personality traits predict one-third of the  need for cognition software engineers.
This evidence has important implications for both practice and research discussed next.

\subsection{Implications for Research and Practice}

Table~\ref{tab:implications} summarizes the key findings of our study and provides recommendations for software organizations.
First, our study finds that software engineers score, on average, high on the construct need for cognition.
Second, in seeking to predict need for cognition with personality traits, we find that the latter account for ca. 33\% of the variation in need for cognition. This is a crucial finding of this study and has several implications for research and practice.

\hyphenation{experience}
\begin{table*}[!ht]
\caption{Summary of findings and propositions}
\label{tab:implications}
\begin{tabular}{p{1.7cm}p{3cm}p{6.3cm}}
\toprule
Theme & Findings & Propositions \\ \midrule
Prediction & Personality traits of software professionals are good predictors for their need for cognition. & Personality traits of software engineers are able to explain 33\% of the variation in need for cognition. \\ \addlinespace
Recruitment & Personality traits of software engineers predict a high need for cognition. & To attract software professionals, organizations should provide cognitive recruitment messages and not emotional ones. \\
\addlinespace
Teaming & Personality traits of software engineers predict a high need for cognition. & Software team composition with an individual high in need for cognition enhances task-relevant information and collaborative team identification, potentially leading to better team performance. \\
\addlinespace
Working behavior & Personality traits of software engineers predict a high need for cognition. & Software engineers have important characteristics that have to be embraced and reflected by a software organization. In particular, professionals tend to generate new ideas and are open to them. They are sensitive to organizational fairness and have a high work ethic. Also, individual innovation behavior is likely to be high. \\
\bottomrule
\end{tabular}
\end{table*}

Recruitment strategies can be better tailored for software engineers.
A study showed that affective (i.e., feeling-based) recruitment messages on company websites are the most attractive for a general population \citep{kraichy2014tailoring}.
However, for an individual with a high need for cognition, such messages are rather uninteresting.
The same research showed that software professionals are more receptive to cognitive (i.e., thought-based) recruiting messages.
From a practical perspective, software organizations should not devote resources to create an emotional link between themselves and candidates.
Instead, it should focus on communicating fairly upon reasons why they are attractive to future employees.
One example is to describe in recruitment advertisements work‐life benefits since they positively influence applicant perceptions of the software company
\citep{firfiray2017lure}.

Moreover, regarding software teaming, the need for cognition predicts team performance \citep{kearney2009and}.
According to Kearney et al., a team with a high need for cognition can better process task-relevant information.
For example, such a team will better understand the customer's needs of the technology stack they are supposed to work on.
Also, the collaborative team identification will be higher since the team cohesion will not depend on emotions but cognitive motivations.
These two points are significant predictors of how a team can perform, i.e., delivering quality software.

Further, a high need for cognition of software engineers likely influences also working behaviors.
The working ethic of an individual with a high need for cognition is above-average \citep{singer1998consideration}.
Research by Singer et al. highlighted employees' high level of empathy with a high need for cognition with potential victims (i.e., colleagues who have been unjustly treated) of their organizations, which predicted the likelihood of whistleblowing behavior and harmful consequences for the company.
Similarly, \cite{pohling2019moral} suggests that individuals' high need for cognition increases active cooperation and pro-social behaviors within organizations.
Thus, if developers' working environment is perceived as unfair, unwelcome, or toxic by developers, they will likely start looking for another workplace and increase their whistleblowing intention \citep{valentine2019moral}.
Consequently, an organizational culture based on fairness, transparency, and accountability is considered very important for subjects high in need of cognition, such as software professionals.
Accordingly, software companies that embrace such organizational behaviors will more likely experience a higher sense of belonging and community of their employees.

Likewise, although individuals such as software engineers are open to change, they do not follow it without questioning it \citep{haugtvedt1992personality}.
They tend to have a critical approach, especially towards contradictory or ill-defined problems.
Therefore, whatever changes a software organization might be willing to introduce, it must ensure that such changes have been communicated sufficiently precisely so that software engineers understand the need and benefits of such change; otherwise, they will be reluctant to change stick to old habits.

Another appealing implication to bear in mind is a high degree of creativity at work \citep{madrid2016creativity}.
Cognitive-laden individuals are very likely to generate new ideas.
Interestingly, as \cite{madrid2016creativity} found, a fair working environment leads to novel idea implementation.
These findings are essential in creative jobs, where employees have to develop new solutions, such as software development.
A fair and transparent organizational culture plays an essential role in delivering customers with novel solutions.
Similarly, Wu et al. reported that the need for cognition is an antecedent of individual innovation behaviors.
In other words, the authors suggest that the high need for cognition explains innovation behaviors.
This results in a higher degree of job autonomy and the ability to better handle time pressure \citep{wu2014need}.
It is also worth mentioning that people scoring high on need for cognition are less prone to the so-called primacy effect: they recall presented information comprehensively and uniformly and are not biased by the first items of a list, showing a higher mental flexibility~\citep{ahlering1989need}.

\subsection{Threats to Validity}
We now discuss threats to validity following \cite{gren2018standards}'s validity concerns for behavioral software engineering studies.

\textit{Reliability}.
We only used stable and reliable measurement instruments that have been validated by social psychology scholars and used in a large number of prior studies. For example, HEXACO has been used in almost 1,000 studies.\footnote{www.hexaco.org/reference}

\textit{Construct validity}.
Our three measurement instruments reflect the constructs for which they have been developed.
Constructs are well-established in the scholarly debate.
Instruments have been refined and improved over time.
We acknowledge having used short versions of the measurement instruments, which were slightly less accurate.
Nonetheless, all three instruments can estimate with high accuracy more extended inventories.
This permitted us to be parsimonious and to use multiple psychometric models (i.e., bright, dark traits, and the need for the cognition of software developers), with a high engagement rate and a low drop-out ratio in the study.

\textit{Conclusion validity}.
Bayesian statistics help us to gather a more comprehensive understanding of the subject matter, overcoming the typical shortcomings of \textit{p}-values-based findings of frequentist null-hypothesis significance testing.
Conclusions are drawn based on the odds related to the best predictive model and not just by answering yes or no to a specific hypothesis.
We also stress the fact that we did not observe the discussed implications directly.
Instead, we inferred our findings with relevant literature on the subject matter.
Accordingly, we developed four propositions that need further investigation.
We artificially generated normal distributions of previous studies to compute the ANOVA since we did not have access to such studies'  data for replication.

\textit{External validity}.
External validity is a crucial concern when conducting sample studies, which seek to maximize the potential for generalizability \citep{stol2018abc}.
In order to establish a sample size that is representative of the software developer population, we relied on \cite{yamane1973statistics}'s formula.
While there are several estimates of the size of the global software developer pool \citep{russo2020gender}, we used a conservative (i.e., the highest) estimate, which was approximately 36.5 million at the time of our study. Following Yamane's formula, a sample of ca. 400 was suggested. Our sample is well above that minimum.
We put the highest care in our multi-stage data collection process using a cluster-randomized sampling strategy, through which we gathered 483 validated responses.
Such a multi-stage process aimed to identify software professionals in our cluster.
As a cluster, we used the academic crowdsourcing platform Prolific.
Although cluster-randomized sampling is considered less efficient than standard-randomized sampling, it provides a fair representation of a population and is very cost-effective \citep{gravetter2018research}.
This approach also follows the most recent recommendations for sampling strategies in software engineering research \citep{baltes2020sampling}.

\section{Conclusion}
\label{sec:conclusion}

The individual motivation to pursue a software engineering career is a puzzling issue to address since there are so many confounding factors that are very hard to control, e.g., the socio-economic background, individual experiences, and aspirations.
Nevertheless, some specific psychological characteristics are intuitively linked to any software professional, such as the joy for problem-solving and, more generally, thinking.
Therefore, we used in this study need for cognition (i.e., a tendency to engage and enjoy effortful thinking \citep{cacioppo1982need}) as a regressor of individual defining or innate characteristics, namely personality traits.
Our analysis shows that the personality traits of software engineers have a considerable explanatory power with a coefficient of determination (\textit{R}\textsuperscript{2}) of ca. 0.33 of the need for cognition.
Thus, individuals who have a high need for cognition are more likely to start a career in software engineering than those who do not.
This finding also has several managerial propositions for software organizations since software professionals are analytical subjects, likely prefer cognitive communication with a transparent, accountable, and fair organizational culture.

Most of the insights we have reported might be seen as self-evident or already apparent, at least to some practitioners.
Nevertheless, we wish to stress that so far, those insights had no scholarly grounding.
Therefore, we put considerable effort into recruiting a large sample of almost 500 software professionals to claim generalization of results, who have been screened throughout a multi-stage selection process, using a cluster-randomized sampling strategy.
Data have been analyzed through Bayesian linear regression to provide a solid understanding of the studied phenomenon.

This is regression analysis, and it does not assess any causal relationship between independent variables and dependent opens.
Nevertheless, our research design implicitly assumes a causal-predictive assumption of our theoretical hypothesis that has been validated.
Also, we observe that the general population scores significantly lower than software engineers in need of cognition and that the personality traits of developers can predict the individual need for cognition.
Accordingly, we speculate that people with a high need for cognition are more likely to start a career in computer science than those who score low.
In other words, there are some individual characteristic that will lead someone more naturally to a computer science career compared to others.
This insight is, for now, an educated guess and needs to be carefully validated with a proper research design to assess causation.

Although we grounded our recommendations in well-established social psychology and management literature, we have not reported them directly with this study.
Therefore, further research should investigate the suggested relations in a software engineering population.

\section*{Supplementary Material}
Raw data, survey questions, and additional material are openly available under a CC-BY 4.0 license at the following URL:
\url{https://osf.io/n57as/}.

\section*{Acknowledgments}
This work was supported by the Science Foundation Ireland grant 13/RC/2094\_P2 and 15/SIRG/3293 and co-funded under the European Regional Development Fund through the Southern \& Eastern Regional Operational Programme to Lero---the Irish Software Research Centre (www.lero.ie). For the purpose of Open Access, the authors have applied a CC-BY public copyright licence to any Author Accepted Manuscript version arising from this submission.

\section*{Conflict of Interest}
All authors declare that they have no conflicts of interest.

\bibliographystyle{spbasic}