e1 <- as.data.frame(table(subset(base, EDUC1 == 1)$MARRIED)/nrow(base))
e2 <- as.data.frame(table(subset(base, EDUC2 == 1)$MARRIED)/nrow(base))
e3 <- as.data.frame(table(subset(base, EDUC3 == 1)$MARRIED)/nrow(base))
e4 <- as.data.frame(table(subset(base, EDUC4 == 1)$MARRIED)/nrow(base))

table(base)
df <- data.frame(cbind(e1,e2,e3,e4))
df <- df[c(2,4,6,8)]
df <- rbind(df, colSums(df))
df <- cbind(df, rowSums(df))
colnames(df) <- c("EDUC1","EDUC2","EDUC3","EDUC4", "Total")
rownames(df) <- c('0', '1', 'Total')
View(df)