source("Clinical.R") descriptive.stats <- function(d) { cat(ncol(d), "columns by ", nrow(d), "rows\n") for (i in 1:ncol(d)) { cat("===============================================================\n") cat(i, " \"", names(d)[i], "\" ", typeof(d[,i]), "\n", sep="") if (typeof(d[,i]) == "double") { print(summary(d[,i])) hist(summary(d[,i]), main=paste(i, names(d)[i]), xlab="", breaks=50) } if (typeof(d[,i]) == "character") { print( table(d[,i]) ) } } } sink("BloodEvaluation.txt") pdf("BloodEvaluation.pdf") par(mfrow=c(3,2)) descriptive.stats(blood) sink() dev.off() sink("Illness.txt") pdf("Illness.pdf") par(mfrow=c(3,2)) descriptive.stats(illness) sink() dev.off()