The performance of (i) the reduced multireference (RMR) coupled-cluster (CC) method with singles and doubles (RMR CCSD) that employs a modest-size MR CISD wave function as an external source for the most important (primary) triples and quadruples in order to account for the nondynamic correlation effects in the presence of quasidegeneracy, (ii) the RMR CCSD(T) method that adds a perturbative correction for the remaining (secondary) triples to the RMR CCSD energy, and (iii) the recently developed partially linearized MR CCSD method, which determines primary triples and quadruples using a subset of linear CC equations projected onto the corresponding higher-than-doubly excited configurations, are tested by considering the singlet-triplet splitting for several diradicals, ranging from a prototypical methylene radical to trimethylenemethane, and benzyne and pyridynium cation isomers. Both RHF and multiconfigurational self-consistent field molecular orbitals are employed. The equilibrium geometries for the lowest-lying singlet and triplet states are determined using both the density functional theory (DFT) and various CC approaches, and a comparison with both the experiment and other theoretical results, wherever available, is made. The RMR CCSD(T) results provide the most satisfactory description in all cases. The dependence of the MR diradical character on a spatial separation of radical centers, as well as the artifactual DFT geometry in the case of benzyne and pyridynium meta-isomers, is also pointed out.