Find the systematic component in the linear form for fitted values in across each simulation (note largely for internal use in qi_builder)

Description

Find the systematic component in the linear form for fitted values in across each simulation (note largely for internal use in qi_builder)

Usage

1
linear_systematic(b_sims, newdata, inc_intercept = TRUE)

Arguments

b_sims

a data frame created by b_sim of simulated coefficients

newdata

a data frame of fitted values with column names corresponding to variable names in b_sims. Variables in b_sim not present in newdata will be treated as fitted at 0. Interactions will automatically be found if they were entered into to the model using the * operator.

inc_intercept

logical whether to include the intercept in the lineary systematic component.

Value

A data frame fitted values supplied in newdata and associated linear systematic component estimates for all simulationed coefficient estimates. The linear systematic components are included in a column named ls_.

Source

King, Gary, Michael Tomz, and Jason Wittenberg. 2000. "Making the Most of Statistical Analyses: Improving Interpretation and Presentation." American Journal of Political Science 44(2): 341-55.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
library(car)

# Estimate model
m1 <- lm(prestige ~ education + type, data = Prestige)

# Create fitted values
fitted_df <- expand.grid(education = 6:16, typewc = 1)

# Simulate coefficients
m1_sims <- b_sim(m1, nsim = 1000)

# Find linear systematic component for fitted values
ls <- linear_systematic(b_sims = m1_sims, newdata = fitted_df)