A function to display
stats::predict's standard errors
ggally_nostic_se_fit( data, mapping, ..., lineColor = brew_colors("grey"), linePosition = NULL )
|data, mapping, ..., lineColor||
parameters supplied to
base comparison for a perfect fit
ggplot2 plot object
As stated in
If the logical 'se.fit' is 'TRUE', standard errors of the predictions are calculated. If the numeric argument 'scale' is set (with optional ''df'), it is used as the residual standard deviation in the computation of the standard errors, otherwise this is extracted from the model fit.
Since the se.fit is
TRUE and scale is unset by default, the standard errors are extracted from the model fit.
A base line of 0 is added to give reference to a perfect fit.
# Small function to display plots only if it's interactive p_ <- GGally::print_if_interactive dt <- broomify(stats::lm(mpg ~ wt + qsec + am, data = mtcars)) p_(ggally_nostic_se_fit(dt, ggplot2::aes(wt, .se.fit)))