extract rse from lm in r

# Fit a linear regression model
model <- lm(y ~ x, data = your_data)

# Extract residuals
residuals <- residuals(model)

# Calculate the residual standard error (RSE)
rse <- sqrt(sum(residuals^2) / (length(residuals) - 2))