How many groups of pentagonal flower bouquets can be formed? 
A florist has three types of flowers: tulips, roses, and daisies. There are 4 tulips, 5 roses, and 6 daisies. These 15 flowers are to be arranged into three bouquets of 5 flowers each. Assume that

*

*the order of the three bouquets is irrelevant,

*flowers of the same type are indistinguishable.

How many groups of pentagonal bouquets can the florist bundle?

Attempt
Let us denote tulips, roses, and daisies with T, R, and D, respectively. If we form all 15-letter strings and add dashes after every five letters, we can obtain all possible groups of bouquets. For instance, one possibility would be
$$\mathrm{TRRTR-TRDDD-DDTRD}.\tag{ex. 1}$$
There are $\dfrac{15!}{4!\ 5!\ 6!}$ such strings. While, of course, all groups of bouquets can be obtained this way, we are overcounting. For strings, $\mathrm{TRDDD-TRRTR-DDTRD}$ is different from the example above, yet it makes no difference for the group of bouquets since order was assumed to be irrelevant. It might be tempting to divide the number of strings by $3!$ but this would also be incorrect. As an example, $\mathrm{TDDDT-TDDDT-RRRRR}\tag{ex. 2}$ is a valid group of three bouquets that should instead be divided by $\dfrac{3!}{2!} = 3$.
So, one way to proceed is to divide all groups of bouquets into two non-intersecting classes. First those for which all groups of three bouquets are pairwise different, and then those with exactly two matching bouquets out of three. Note that forming groups with three identical bouquets is impossible because 4 tulips cannot be shared equally among three bouquets. Once partitioned in this manner, we may divide the first kind of partition with $3!$, and the second with $3$.
However, such partitioning seems overly tedious, and is further complicated by the following aspect. We still have to consider that whenever there are at least two different kinds of flowers in a single bouquet, there is a further overcounting with the string-approach. Namely, e.g., the bouquets
$$\mathrm{TRDDD\equiv DTRDD\equiv DDTRD\equiv DDDTR\equiv RDDDT}\tag{ex. 3}$$
are all equivalent since they can be transformed into one another by a rotation in space. (So a division with $5$ might additionally be in order for such bouquets). The 'further complication' is then the fact that groups of bouquets which initially seem to warrant division by $3!$, actually require division with $3$, as is the case for our first example. Indeed, by ex. 3 we have $\mathrm{TRDDD\equiv DDTRD}$ and so
$$\mathrm{TRRTR-TRDDD-DDTRD\equiv TRRTR-TRDDD-TRDDD}$$ which should by divided by $3$.
Clarification from comments: bouquets which can be transformed into one another by reflection are not equivalent, and should be counted as different bouquets.
Question
The above discussion seems to lead into various subcases where mistakes could be easy to make, and is tedious to generalise. Is there cleaner approach? Regardless, an answer which carefully carries the above scheme to completion has value, too. For the record, the answer I obtain with the above method is $898$.
Edit: I have now also "confirmed" the answer $898$ with an independent Python program.
Trying to find partitions of the multiset $\{\mathrm{T}:4, \mathrm{N}:5, \mathrm{D}:6\}$ into classes of size five is something that I admittedly have not thought much about, but on the face of it would lead to undercounting, since, for instance, the multiset $\{\mathrm{D, D, R, R, T}\}$ would not differentiate between non-equivalent bouquets $\mathrm{DDRRT}$ and $\mathrm{DTDRR}$.
(This question is from the context of introductory combinatorics without recurrences, generating functions and so on).
 A: Not an answer, but just a long comment, regarding the partition of the multiset you mention, thus when we don't consider the order of flowers in bouquets. However, the order of bouquets is taken into account here, so that e.g. TTDDD-TTDDD-RRRRR is different from TTDDD-RRRRR-TTDDD.
I don't think this calculation can help to solve your problem, anyway.
We can define $t_1, t_2, t_3$ the number of tulips in bouquet $1,2,3$ respectively, $r_1, r_2, r_3$ the number of roses in bouquet $1,2,3$ respectively, $d_1, d_2, d_3$ the number of daisies in bouquet $1,2,3$ respectively. We need then to find the number of non-negative integer solutions to the following system of equations:
$$\begin{cases}
t_1+r_1+d_1=5 \\
t_2+r_2+d_2=5 \\
t_3+r_3+d_3=5 \\
t_1+t_2+t_3=4 \\
r_1+r_2+r_3=5 \\
d_1+d_2+d_3=6
\end{cases}$$
Then we can use generating functions, and assign the following variables to each equation: $x$ to equation $1$, $y$ to equation $2$, $z$ to equation $3$, $t$ to equation $4$, $u$ to equation $5$, $v$ to equation $6$.
The generating function is built with one factor for each variable of the linear system:
$$f(x,y,z,t,u,v)=\frac{x^6t^6-1}{xt-1}\frac{x^6u^6-1}{xu-1}\frac{x^6v^6-1}{xv-1}\frac{y^6t^6-1}{yt-1}\frac{y^6u^6-1}{yu-1}\frac{y^6v^6-1}{yv-1}\frac{z^6t^6-1}{zt-1}\frac{z^6u^6-1}{zu-1}\frac{z^6v^6-1}{zv-1}
$$
For example, $\frac{x^6t^6-1}{xt-1}$ takes into account equation $1$ and $4$ for the variable $t_1$. We limit the terms to $x^5$ because this is the total of the first equation.
We use then WolframAlpha to compute the product of the first three factors, see here. We are interested only on the term with $x^5$ (the other factors don't contribute to $x$), which is:
$$x^5 (t^5 + t^4 (u + v) + t^3 (u^2 + u v + v^2) + t^2 (u^3 + u^2 v + u v^2 + v^3) + t (u^4 + u^3 v + u^2 v^2 + u v^3 + v^4) + u^5 + u^4 v + u^3 v^2 + u^2 v^3 + u v^4 + v^5)$$
and we can add the other two factors to get:
$$x^5 y^5 z^5 (t^5 + t^4 (u + v) + t^3 (u^2 + u v + v^2) + t^2 (u^3 + u^2 v + u v^2 + v^3) + t (u^4 + u^3 v + u^2 v^2 + u v^3 + v^4) + u^5 + u^4 v + u^3 v^2 + u^2 v^3 + u v^4 + v^5)^3$$
Then, we need to find the coefficient of $x^5y^5z^5t^4u^5v^6$. To do that, we still use WolframAlpha, but note that I had to replace $t,u,v$ with $x,y,z$ there only to make the engine understand the query, so this is a little confusing, but think of having $t,u,v$ there.
The result is $210$, I hope to have done all correctly.
I don't know if there is a way to do this by hand in an elegant way.
I noticed that $\frac{15!}{4!5!6!}=630630=210 \times 3003$ but I don't know if this might mean something.
