The β function in higher covariant derivative regularisation