How to best rewrite and fit a chunk of R code to lapply()?

Question

I have some SAS programming experience, but not so much in R.

I have been using the following code successfully to carry out regressions on a large twin sample: