Let $a_n$ denote the answer. I claim that
$\sum_{n \ge 0} \frac{a_n}{n!} z^n = \frac{1}{2 - e^z}.$
To see this, write the RHS as $\frac{1}{1 - (e^z - 1)}$, or
$\sum_{k \ge 0} (e^z - 1)^k.$
Then I claim that the coefficient of $\frac{z^n}{n!}$ in $(e^z - 1)^k$ is precisely the number of ways to weakly order $n$ objects such that there are $k$ equivalence classes. This is equivalent to the number of ways to partition $n$ objects into $k$ non-empty subsets in order, so you can prove this by using identities for Stirling numbers, but it actually follows directly from basic properties of exponential generating functions: $e^z - 1$ is the generating function for non-empty sets, so $(e^z - 1)^k$ is the generating function for $k$-tuples of non-empty sets.
The generating function above shows that one can't expect, say, a formula in terms of a bounded number of binomial coefficients or anything especially nice like that. However, $\frac{1}{2 - e^z}$ is a meromorphic function, and the pole closest to the origin is at $z = \log 2$. This turns out to imply, after some computations, that
$\frac{a_n}{n!} \sim \frac{1}{2 (\log 2)^{n+1}}$
and one can even precisely write down the rest of the terms in the asymptotic expansion. For example, the next poles are at $z = \log 2 \pm 2 \pi i$, so the next term in the asymptotic expansion is $\frac{a_n}{n!} \sim \frac{1}{2 (\log 2)^{n+1}} + \frac{\cos (n+1) \theta}{r^{n+1}}$
where $r = \sqrt{(\log 2)^2 + 4 \pi^2}, \theta = \arctan \frac{2\pi}{\log 2}$. It actually follows from the fact that all of the other terms are much smaller than this that the first term in the above estimate has exponentially small error. For example, the largest value in the OEIS is $\frac{a_{18}}{18!} = 528.793...$
whereas $\frac{1}{2 (\log 2)^{19}} = 528.794...$
For more on these kinds of techniques, see Flajolet and Sedgewick's Analytic Combinatorics.