Box.test(residuals(model.tslm), type=“Ljung-Box”, lag=20)

Box-Ljung test

data: residuals(model.tslm) X-squared = 48.55, df = 20, p-value = 0.000356