Rationale

The elaboration likelihood model states that people are cognitive misers and have plenty of energy to spend processing information, but are reluctant to do so. There are two routes to attitude change within the model, central and peripheral. The central route contains high-quality arguments that are information-rich and can be difficult to process, which creates a change in cognitive structure, an enduring attitude change. However, recipients must be motivated and able to process the information. The peripheral route features peripheral cues that require less cognitive effort and lead to short-term attitude changes. This typically happens because the hurdles of processing through the central route are too taxing for the recipient, or they are not persuaded and have not changed their cognitive structure; they default to this type of processing.

Considering the model, eye-tracking technology measuring the average number of minutes spent looking at “high-quality” persuasive information arguing for legalizing recreational marijuana use in Tennessee may produce long-term attitude change.

Hypothesis

The amount of time spent looking at “high quality” persuasive information arguing for legalizing recreational marijuana use in Tennessee is associated with the likelihood of whether the participant favored (1) or opposed (0) legalizing recreational marijuana in Tennessee.

Variables & method

The continuous independent variable in the analysis was the amount of time, in minutes, the respondent spent looking at (and, presumably, reading) “high-quality” persuasive information arguing for legalizing recreational marijuana use in Tennessee.

The dependent variable was categorical with two categories as an indication of whether the participant favored (1) or opposed (0) legalizing recreational marijuana in Tennessee. Participants who were undecided or who declined to answer were also coded as 0.

A logistical regression analysis was conducted to test the odds that the longer the respondent spends looking at “high quality” persuasive information arguing for the legalization of recreational marijuana use in Tennessee, the more likely they are to favor (1) legalization.

Results & discussion

The results supported the hypothesis. The probability of being in favor of the legalization of marijuana in Tennessee increases for every minute spent looking at “high-quality” persuasive information arguing for the legalization of marijuana in Tennessee. In particular, after about 15 minutes of exposure to the persuasive information, respondents were more likely than not to favor legalizing marijuana in Tennessee. The dataset below is statistically significant with a p-value of 0.6869.





Logistic Regression Results
Odds Ratios with 95% Confidence Intervals
term Odds_Ratio CI_Lower CI_Upper P_Value
(Intercept) 0.099 0.044 0.205 0.0000
IV 1.167 1.116 1.227 0.0000



Linearity of the Logit Test (Box-Tidwell)
Interaction term indicates violation if significant
term Estimate Std_Error P_Value
(Intercept) −2.566 1.100 0.0196
IV 0.262 0.290 0.3657
IV_log −0.032 0.078 0.6869



Inflection Point of Logistic Curve
Value of IV where predicted probability = 0.50
Probability Inflection_Point
0.5 14.965



CODE

# ------------------------------
# Install and load required packages
# ------------------------------
if (!require("tidyverse")) install.packages("tidyverse")
if (!require("gt")) install.packages("gt")
if (!require("gtExtras")) install.packages("gtExtras")
if (!require("plotly")) install.packages("plotly")

library(ggplot2)
library(dplyr)
library(gt)
library(gtExtras)
library(plotly)


# ------------------------------
# Read the data
# ------------------------------
mydata <- read.csv("ELM.csv") # <-- EDIT filename

# ################################################
# # (Optional) Remove specific case(es)s by row number
# ################################################
# # Example: remove rows 10 and 25
# rows_to_remove <- c(10, 25) # Edit and uncomment this line
# mydata <- mydata[-rows_to_remove, ] # Uncomment this line

# Specify dependent (DV) and independent (IV) variables
mydata$DV <- mydata$Favor_1   # <-- EDIT DV column
mydata$IV <- mydata$Minutes   # <-- EDIT IV column

# Ensure DV is binary numeric (0/1)
mydata$DV <- as.numeric(as.character(mydata$DV))


# ------------------------------
# Logistic regression plot 
# ------------------------------
logit_plot <- ggplot(mydata, aes(x = IV, y = DV)) +
  geom_point(alpha = 0.5) +   # scatterplot of observed data
  geom_smooth(method = "glm",
              method.args = list(family = "binomial"),
              se = FALSE,
              color = "#1f78b4") +
  labs(title = "Logistic Regression Curve",
       x = "Independent Variable (IV)",
       y = "Dependent Variable (DV)")

logit_plotly <- ggplotly(logit_plot)


# ------------------------------
# Run logistic regression
# ------------------------------
options(scipen = 999)
log.ed <- glm(DV ~ IV, data = mydata, family = "binomial")

# Extract coefficients and odds ratios
results <- broom::tidy(log.ed, conf.int = TRUE, exponentiate = TRUE) %>%
  select(term, estimate, conf.low, conf.high, p.value) %>%
  rename(Odds_Ratio = estimate,
         CI_Lower = conf.low,
         CI_Upper = conf.high,
         P_Value = p.value)

# Display results as a nice gt table
results_table <- results %>%
  gt() %>%
  fmt_number(columns = c(Odds_Ratio, CI_Lower, CI_Upper), decimals = 3) %>%
  fmt_number(columns = P_Value, decimals = 4) %>%
  tab_header(
    title = "Logistic Regression Results",
    subtitle = "Odds Ratios with 95% Confidence Intervals"
  )


# ------------------------------
# Check linearity of the logit (Box-Tidwell test)
# ------------------------------
# (Assumes IV > 0; shift IV if needed)
mydata$IV_log <- mydata$IV * log(mydata$IV)
linearity_test <- glm(DV ~ IV + IV_log, data = mydata, family = "binomial")

linearity_results <- broom::tidy(linearity_test) %>%
  select(term, estimate, std.error, p.value) %>%
  rename(Estimate = estimate,
         Std_Error = std.error,
         P_Value = p.value)

linearity_table <- linearity_results %>%
  gt() %>%
  fmt_number(columns = c(Estimate, Std_Error), decimals = 3) %>%
  fmt_number(columns = P_Value, decimals = 4) %>%
  tab_header(
    title = "Linearity of the Logit Test (Box-Tidwell)",
    subtitle = "Interaction term indicates violation if significant"
  )


# ------------------------------
# Calculate the inflection point (p = .50)
# ------------------------------
p <- 0.50
Inflection_point <- (log(p/(1-p)) - coef(log.ed)[1]) / coef(log.ed)[2]

inflection_table <- tibble(
  Probability = 0.5,
  Inflection_Point = Inflection_point
) %>%
  gt() %>%
  fmt_number(columns = Inflection_Point, decimals = 3) %>%
  tab_header(
    title = "Inflection Point of Logistic Curve",
    subtitle = "Value of IV where predicted probability = 0.50"
  )


# ------------------------------
# Outputs
# ------------------------------
# Interactive plot
logit_plotly

# Tables
results_table
linearity_table
inflection_table