You seem to be double-counting the elements of $S_4$ with cycle type $(2, 2)$: There are only three, namely, $(12)(34)$, $(13)(24)$, $(14)(23)$.
More generally, a permutation $\sigma \in S_n$ satisfies $\sigma^2 = e$, i.e., is an involution, if its cycle type contains only $1$'s and $2$'s, that is, if it is a product of (possibly zero) disjoint transpositions.
The number of permutations with $k$ transpositions and $n - 2 k$ fixed points is, in multinomial notation,
$$\frac{1}{k! (n - 2 k)!}{{n}\choose{1, \ldots, 1, 2, \ldots, 2}} = \frac{n!}{2^k k! (n - 2 k)!} :$$ The multifactorial counts the number of ways to put $n$ labeled elements in $n - k$ bins so that $n - 2 k$ bins have $1$ element and $k$ bins have $2$ elements. The labeling on the bins of a given size are irrelevant, so we divide the multinomial by $k! (n - 2 k)!$ as indicated.
For example, for the given case $n = 4$, we can have
- $k = 0$ transpositions (there are $\frac{1}{0! 4!} {4 \choose{1, 1, 1, 1}} = 1$ of these, namely the identity),
- $k = 1$ transpositions (there are $\frac{1}{1! 2!} {4 \choose{1, 1, 2}} = 6$ of these, the transpositions themselves), or
- $k = 2$ transpositions (there are $\frac{1}{2! 0!} {4 \choose{2, 2}} = 3$ of these, the double transpositions).
Summing over the possible values of $k$ gives a general formula number $T_n$ of involutions in $S_n$:
$$T_n = \sum_{k = 0}^{\lfloor\frac{n}{2}\rfloor} \frac{n!}{2^k k! (n - 2 k)!} .$$
The first few values are $T_1 = 1$, $T_2 = 2$, $T_3 = 4$, $T_4 = 10$, $T_5 = 26$.
This turns out to be a well-known sequence, A000085, Number of self-inverse permutations on $n$ letters, also known as involutions; number of standard Young tableaux with $n$ cells. See that link for other interpretations and various formulas for and identities of this sequence.
For more references, see the second paragraph of the introduction of this 1957 article by Miksa, Moser, and Wyman.