This paper presents an implementation in R of the Cluster Estimated Standard Errors (CESE). The method estimates the covariance matrix of the estimated coefficients of linear models in grouped data sets with correlation among observations within groups. Cluster Estimated Standard Errors (CESE) is an alternative solution for the classical Cluster Robust Standard Errors (CRSE), which underestimates the standard errors in most of the situations encountered in practice.