Hi,

I am trying to do association analysis between soil toxic elements and epigenetic age based on road buffer zone by using below code; but could you please suggest any other modelling or data transformation if needed for this analysis. I am just using raw ICP-MS value(mg/kg) toxic element value and epigenetic age acceleration calculated by regressing again chronological age.

raw <- read.table("clean.raw.elements.txt", header=T, sep="\t")
colnames(raw)


model1 <- lm (EpigeneticAgeAcceleration_EEAA ~ 0 + buffer.zone+ buffer.zone:( Zn+ Cu +Hg + Mo +Pb + Smoking_Status + Sex + Age), data = raw)
summary_model1 <- summary(model1)