made up example list wise PCA with AIC

pca_example.R

# just an implementation example that does not make a whole lot of sense... 

library(MASS)
y <- mtcars$mpg
ll <- list(cars1 = mtcars[,-1], cars2 = mtcars[,-1]*100)
out_put_list <- lapply(ll,function(x){
  pca <- prcomp(x,scale.=T)
  pca_scores <- as.data.frame(pca$x)
  # maybe give back the result and do it 
  # multiple dfs
  fit1 <- lm(y ~ .,data =pca_scores)
  best_model <- stepAIC(fit1)
  best_model
})