...
| Code Block |
|---|
prs1phen=read.csv("test_clump_1_prs_5.csv",header=T)
head(prs1phen)
lmprs<-lm(data=prs1phen, dummybmi~SCORE + sex) #adding sex as a covariate here; consider adding population structure PCs as covariates
summary(lmprs) |
6.2- PLOTTING
...