This report analyzes a 2³ full factorial experiment designed to understand how water–cement ratio, curing time, and additive type affect concrete compressive strength.
The dataset is provided separately as
concrete_strength_sim.csv.
df <- read.csv("concrete_strength_sim.csv")
head(df)
df$A <- df$A_coded
df$B <- df$B_coded
df$C <- df$C_coded
model <- lm(Strength_MPa ~ A * B * C, data=df)
summary(model)
library(car)
Anova(model, type=2)
par(mfrow=c(2,2))
plot(model)
shapiro.test(resid(model))
newdat <- expand.grid(A=c(-1,1), B=c(-1,1), C=c(-1,1))
newdat$pred <- predict(model, newdat)
transform(newdat,
A_wcr = ifelse(A==-1, 0.40, 0.55),
B_days = ifelse(B==-1, 7, 28),
C_additive = ifelse(C==-1, "None", "Plasticizer")
)