A method is presented for computing maximum likelihood, or Gaussian, estimators of the structural parameters in a continuous time system of higherorder stochastic differential equations. It is argued that it is computationally efficient in the standard case of exact observations made at equally spaced intervals. Furthermore it can be applied in situations where the observations are at unequally spaced intervals, some observations are missing and/or the endogenous variables are subject to measurement error. The method is based on a state space representation and the use of the Kalman–Bucy filter. It is shown how the Kalman-Bucy filter can be modified to deal with flows as well as stocks.