################## ERSTE EINFACHE BEISPIELE ZUM UMGANG MIT R ##################

# Das Doppelkreuz # fungiert als Kommentar-Zeichen; Text dahinter wird von R
# ignoriert; gut zum Kommentieren von eigenen Programmen geeignet

########### R als Taschenrechner ##############################################

(5.2-3.4)/2
log(2.718)

########### Zuweisung #########################################################

### Wert rechts vom "=" wird dem Element links vom "=" zugewiesen

a=5.2
b=-3.4
c=2
(a+b)/c

########### Erzeugen einfacher Beispiel-Daten #################################

# Folge von 1 bis 10 (ganze Zahlen)
x=1:10
x

# Folge von 2 bis 18 im Abstand 4
y=seq(2,18,4) 
y

# Folge von 1.2 bis 0.2 im Abstand -0.1
z=seq(1.2,0.2,-0.1)
z

# Laenge von z
length(z)

# Folge von 1.2 bis 0.2 im Abstand -0.01
x1=seq(1.2,0.2,-0.01)
x1

# Laenge von x1
length(x1)

# irgendeine beliebige Folge (z.B. Klausurnoten)
noten=c(3,5,3,2,4,4,3,2,2,1)

# Zugriff auf einzelne Elemente
noten[1] # 1.Element
noten[5] # 5.Element

# Vektor noten ordnen: ergibt die Positionen, wo die Elemente
# der Reihe nach auftauchen
order(noten)   

### Ergebnis:  [1] 10  4  8  9  1  3  7  5  6  2
 
# geordneter Vektor
notenORD=noten[order(noten)]
notenORD

### Ergebnis:  [1] 1 2 2 2 3 3 3 4 4 5

########### Einfaches erstes Histogramm #######################################

hist(noten)  ### ACHTUNG: die linke Saeule repraesentiert die 1er UND 2er! 

hist(noten,col=2)  ### dasselbe Histogramm mit etwas Farbe 

hist(noten,col=2,right=FALSE)  ### nun repraesentiert die rechte Saeule die 
### 4er UND 5er!

hist(noten,col=2,breaks=seq(0.5,5.5,1)) # schoenes Histogramm (eine Saeule fuer
### jede Note); mit 'breaks' die Bruchpunkte fuer die Klassen angeben

hist(noten,col=2,breaks=c(0.5,1.5,2.5,3.5,5.5)) ### Zusammenlegung 4er und 5er
### ACHTUNG: jetzt erscheint 'density' auf der y-Achse: Die Flaechen geben
### jetzt die relativen Haeufigkeiten wieder

########### Einlesen von Daten aus Dateien ####################################

### wir gehen davon aus, dass im Arbeitsverzeichnis eine Datei mit dem Namen 
### income.txt liegt; in der einzigen Spalte (ohne Beschriftung in der ersten
### Zeile) stehen Einkommensdaten; am besten den genauen Dateipfad eruieren und 
### dann in der Form z.B. 
### "C:\\Users\\ich\\Desktop\\income.txt" 
### an read.table uebergeben 

mydata=read.table(file="income.txt",h=F) ### h=F, weil keine Beschriftung in
### der ersten Zeile (andernfalls h=T verwenden)

income=mydata[,1]  ### Zuweisung der ersten Spalte der Daten an 'income'

length(income)     ### Anzahl der vorhandenen Daten
min(income)        ### Minimum der Daten
max(income)        ### Maximum der Daten

hist(income,col=2,breaks=seq(600,2800,200))  ### Histogramm zeichnen


########### Beenden von R #####################################################

q()

### beim Beenden i.d.R. mit 'y' die Sicherung des Workspace vornehmen

###############################################################################