Background: While many phenomenological models for nuclear fission have been developed, a microscopic understanding of fission has remained one of the most challenging problems in nuclear physics. Purpose: We investigate the applicability of the generator coordinate method (GCM) as a microscopic theory for cluster radioactivities of heavy nuclei, which can be regarded as a fission with large mass asymmetry, that is, a phenomenon in between fission and α-decays. Methods: Based on the Gamow theory, we evaluate the preformation probability of a cluster with GCM while the penetrability of the Coulomb barrier is estimated with a potential model. To this end, we employ Skyrme interactions and solve the one-dimensional Hill-Wheeler equation with the mass octupole field. We also take into account the dynamical effects of the pairing correlation using BCS wavefunctions constructed with an increased strength of the pairing interaction. Results: We apply this scheme to the cluster decay of 222 Ra, i.e., 222 Ra→ 14 C+ 208 Pb, to show that the experimental decay rate can be reproduced within about one order of magnitude. We also briefly discuss the cluster radioactivities of the 228 Th and 232 U nuclei. For these actinide nuclei, we find that the present calculations overestimate the decay rates by several orders of magnitudes. Conclusions: The method presented in this paper provides a promising way to describe microscopically cluster decays of heavy nuclei.