(5), however, exhibit heteroscedasticity: because we are averaging over all time origins within the trajectory, data points at longer times have higher variance. In this case, we can use weighted least squares to fit the data, where the weights are the reciprocal of the variance of the quantity in angular brackets. By providing higher weights for the short-time data (where we have better statistics), we observe consistently lower error bars in our transport coefficient measurements.