In this paper we show how to simulate and estimate a COGARCH(p, q) model in the R package yuima. Several routines for simulation and estimation are introduced. In particular, for the generation of a COGARCH(p, q) trajectory, the user can choose between two alternative schemes. The first is based on the Euler discretization of the stochastic differential equations that identify a COGARCH(p, q) model while the second considers the explicit solution of the equations defining the variance process.Estimation is based on the matching of the empirical with the theoretical autocorrelation function. Three different approaches are implemented: minimization of the mean squared error, minimization of the absolute mean error and the generalized method of moments where the weighting matrix is continuously updated. Numerical examples are given in order to explain methods and classes used in the yuima package.