#erstmal alle variablen löschen
remove(list=ls())
#aktuelles Arbeitsverzeichnis zuweisen
setwd("/Users/michael/Dropbox/Angewandte Statistik SoSe11/Uebungen/RIntro/")
#Zufallszahlengenerator initialisieren
set.seed(23514)

#Werte zuweisen
x <- 4

#Vektoren zuweisen/erzeugen
x <- 1:6
y <- 2*x - 3
#Nullvektor der Länge 5
v <- rep(0,5)
#von 0 nach 1 in Schritten von 0.1
v2 <- seq(0,1,0.1)
#Teilvektor
v3 <- v2[1:5]
#Vektoren aneinander hängen
z <- c(x,v)

t(x) %*% y
x*y

#Matrix
M <- matrix(1:10, nrow=2, ncol=5)
M2 <- matrix(1:10, nrow=2, ncol=5,byrow=TRUE)
#Matrix aus Vektoren (spaltenweise)
M3 <- cbind(x,y,x+y)
#auf Eintrag zugreifen
M3[2,3]
#zweite Zeile von M2
M[2,]

#Daten einlesen (aus der Datei dataIntro.dat im Verzeichnis wie in setwd() oben)
data <- matrix(scan("./dataIntro.dat"),ncol=2,byrow=TRUE)

#Vektoren aus den beiden Spalten der eingelesenen Daten
u1 <- data[,1]
u2 <- data[,2]
#Einträge gegeneinander plotten
plot(u1,u2)
x <- seq(0,1,0.2)
#Die Punkte mit den Koordinaten (x[i] , x[i]*x[i]) mit Linien verbinden und zum Plot hinzufügen
lines(x,x*x)

#Zufallszahlen erzeugen
#Gleichverteilung auf [0,1]
x <- runif(100,min=0,max=1)

#Gleichverteilung auf [-10,10]
x <- runif(100,min=-10,max=10)
#wieviele Elemente enthaelt x?
length(x)
#gerundet auf ganze Zahlen
z <- round(x)
#wo im Vektor z taucht der Wert 3 auf?
z3 <- z==3
#welche Werte wurden auf 3 gerundet?
x[z==3]
#Summe der Eintraege
sum(z)
#Stichprobenmittel
mean(x)
#Zusammenfassung
summary(x)


#if-Abfrage
if (runif(1,-1,1)>0) { 
	print("Die Zufallszahl ist groesser als Null.") 
} else { 
	print("Doch nicht.") 
}

#for-Schleife
res <- 0
for (n in 1:10) {
	print(n)
	res <- res +n
}
#jetzt sollte res=55 gelten.