################################################################################
#
# PhD. European convergence 20-07-2015
# Function. Extract estimates of simple lm models
# NOTE: the function designed for just a simple linear regression with 1 regressor
# Ilya Kashnitsky, ilya.kashnitsky@gmail.com
#
################################################################################
ik_an_extract.lm <- function(list, round=T){
extract <- data.frame(matrix(0,length(list),7))
names(extract) <- c('model_y.x','coef','conf2.5','conf97.5','pValue','Rsq','cor')
for (i in 1:length(list)){
modeli <- list[[i]]
namei <- names(list)[i]
coef <- coef(modeli)[[2]]
ci1 <- confint(modeli)[2,1]
ci2 <- confint(modeli)[2,2]
p.value <- summary(modeli)$coefficients[2,4]
Rsq <- summary(modeli)$r.squared
cor <- cor(modeli$model)[2,1]
extract[i,1] <- namei
extract[i,2:7] <- c(coef,ci1,ci2,p.value,Rsq,cor)
}
if (round==T) {extract[,2:7] <- apply(extract[,2:7],2,round,3)} else {}
return(extract)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.