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 |
---|---|

linePosition | base comparison for a perfect fit |

ggplot2 plot object

As stated in `stats::predict`

documentation:

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)))