rm(list=ls())
library(DAAG)
par(mfrow=c(1,1),ask=F)
pause()
#-----------------------------------------------------------------------------
# Charakterisierung bodenphysikalischer Messungen durch deskriptive Statistik
#-----------------------------------------------------------------------------
#
# Bereitstellung der Datensätze: 
#	asw - Abscherwiderstand (N/m^2)
#	ld  - Trockenraumdichte (g/cm^3)
#	prs - Eindringwiderstand (N/m^2)
#	wsv - Vol. Wassergehalt (% vol)
#
asw <- c(23,31,39,47,55,64,72,80,88,96,104,113,121,129,137,145,153,162,170,178)
ld  <- c(1.15,0.94,1.21,0.99,1.14,1.18,1.04,0.97,1.04,1.11,1.07,1.11,1.10,1.14,1.27,1.37,1.11,1.39,1.25,1.34)
prs <- c(24,19,21,35,14,35,49,40,50,50,56,42,55,54,86,67,88,50,33,119)
wsv <- c(15.6,13.9,13.8,15.5,15.7,14.6,17.3,17.3,17.0,18.8,14.6,15.2,16.7,18.3,18.4,19.8,19.4,18.2,20.2,19.6)
#
# erstellen eines Dataframe (Datenobjekt mit Daten in Spalten und ggf. Überschrift)
#
dfr <- data.frame(asw,ld,prs,wsv)
attach(dfr)
smr <- summary(dfr)
#
par(mfrow=c(1,4))
hist(wsv,prob=T,breaks=5, col="lightblue", border="blue")
hist(ld,prob=T,breaks=5, col="lightblue", border="blue")
hist(prs,prob=T,breaks=5, col="lightblue", border="blue")
hist(asw,prob=T,breaks=5, col="lightblue", border="blue")
pause()
#
boxplot(wsv, notches=F,col="wheat2",border="red",lwd=2)
boxplot(ld, notches=F,col="wheat2",border="red",lwd=2)
boxplot(prs, notches=F,col="wheat2",border="red",lwd=2)
boxplot(asw, notches=F,col="wheat2",border="red",lwd=2)
#
plot(dfr)
crl <- cor(dfr)
#
# Anzeigen der Daten:
#
cat("\nSummary:\n")
cat("========\n")
print(smr)
cat("\nKorrelationsmatrix:\n")
cat("===================\n")
print(crl)

par(mfrow=c(1,1))
detach(dfr)
#
# Aufgaben in der R-Konsole:
# --------------------------
# 1) Lassen Sie sich die im workspace vorhandenen Objekte anzeigen. 
# 2) Überprüfen Sie die Größe der Datenvektoren mit dem Befehl "length()"
# 3) Geben Sie den Spalten des dataframes andere Spaltenüberschriften mit dem
#    Befehl names(df) <- c("name1","name2",...)
# 4) Entfernen Sie aus dem dataframe die prs-Spalte
# 5) Zeichnen Sie die boxplots mit notches
# 6) Beschriften Sie die plots und die Achsen
# 7) Analysieren Sie den pause()-Befehl
# 8) Schreiben Sie den dataframe in eine Datei
#                     
#---Ende---