"create"<- function(filename) { www <- read.table(filename, sep = "\t", header = T) tnames <- names(www) w2 <- as.matrix(www) dimnames(w2) <- list(NULL, NULL) total <- w2[, grep("P0030*", tnames)] all <- w2[, grep("P0130*", tnames)] w.male <- w2[, grep("P014A0*", tnames)] w.female <- w2[, grep("P014B0*", tnames)] b.male <- w2[, grep("P014C0*", tnames)] b.female <- w2[, grep("P014D0*", tnames)] h.male <- w2[, grep("P015A0*", tnames)] h.female <- w2[, grep("P015B0*", tnames)] Ages <- cbind(w.male, w.female, b.male, b.female, h.male, h.female) data.frame(all, Ages) }