ssfit: Fitting of parametric models using summary statistics. Fits complex parametric models using the method proposed by Cox and Kartsonaki (2012) without likelihoods: The fitting of complex parametric models. Consider parametric models that are too complicated to allow calculation of a likelihood but from which observations can be simulated. We examine parameter estimators that are linear functions of a possibly large set of candidate features. A combination of simulations based on a fractional design and sets of discriminant analyses is then used to find an optimal estimator of the vector parameter and its covariance matrix. The procedure is an alternative to the approximate Bayesian computation scheme.