library(tidyverse)

library(agricolae)

library(semPlot)

library(lavaanPlot)

data(haynes)

X <- haynes %>% dplyr::select(FL,   MI,   ME)

Y <- haynes$WI

# correlation analysis ---------------------------------------

corr.x <- correlation(X, X)

corr.y <- correlation(Y, X)

# Path Analysis -----------------------------------------------

result <- path.analysis(corr.x$correlation, corr.y$correlation)

## Direct effect path coefficients

DEffect <- diag(result$Coeff)

## Indirect effect path coefficients

diag(result$Coeff) <- 0

IEffect <- rowSums(result$Coeff)

# print results of path analysis -----------------------------------------------

result1 <- result$Coeff %>% as.data.frame %>%

  mutate(

    '总间接作用' = IEffect,

    '直接作用' = DEffect,

    '总作用' = corr.y$correlation %>% as.numeric()

  ) 

rownames(result1) <- DEffect %>% names

cat("\n\nDirect and indirect effect:\n========================\n" )

print(result1)

sprintf("\nResidual Effect^2 =  %2.3f\n========================\n",

        result$Residual) %>% cat