Code: "ueb" <- function () {
werner_bcd <<- read.table("werner_bcd.txt") werner_bcd <<- matrix(scan("werner_bcd.txt",na.strings="*"), ncol=9,byrow=T) dimnames(werner_bcd) <<- list(werner_bcd[,1], c("Id","Age","Height","Weight","Brthpill","Cholstrl","Albumin","Calcium","Uricacid"))
werner_bcd <- as.data.frame(werner_bcd)
tfvec <- !is.na(werner_bcd[,8]) tfvec <- !is.na(werner_bcd[,9]) & tfvec
Ca <- werner_bcd[tfvec,8] Ur <- werner_bcd[tfvec,9]
par(mfrow=c(3,1))
# 1. Histogramm mit "eindimensionalen Streuungsdiagramm"
hist(Ur, col="darkgreen", breaks=seq(min(Ur),max(Ur)+3,by=3), main="URICACID Verteilung", xlab = "URICACID", ylab="Anzahl" ) rug(Ur, side=1)
# 2. Schätzungen
print(paste("Lokation:")) print(summary(Ur))
print(paste("Streuung:", diff(range(Ur)))) print(paste("Schiefe:", skew(Ur)))
print(paste("Kurtosis:", kurt(Ur)))
#3 Korrelationskoeffizient
print(paste("Korrelationskoeffizient von Calcium und Uricacid: p=",(cov(Ur,Ca)/(var(Ur)*var(Ca)))))
plot(Ca,Ur, xlab="Calcium", ylab="Uricacid") # Streuungsdiagramm
abline(lsfit(Ca,Ur),col="brown") legend(86,100,c("Calcium-Uracid", "Regression"))
#4 Boxplot
boxplot(werner_bcd[,c(2:length(werner_bcd))], col="darkgreen", outline=T) title("Werner Blood Chemistry Data")
} |