####################
# R-harjoitukset 1 #
# viikko 36, 2007  #
####################

# Ohjelman kotisivut:  
# http://r-project.org 

######################
# esimerkkejä aluksi #
######################

# Lukujonojen erilaisia määrityksiä ja nimeäminen:

x <- 1:10
x

1:10-1

1:(10-1)

seq(1,10,0.5)

seq(length=9,from=1, to=10)

rep(1,10)

# R:llä on hyvä help-toiminto! Esim. rep()-funktiosta voi lukea tarkemmin:
help(rep)
#tai
?rep 

y <- sequence(1:4) 
y # huomaa muodostuminen: 1, (1,2) , (1,2,3) , (1,2,3,4)

sequence(c(2,5,7))

z <- c(1,seq(1,10,0.5),rep(1,10))
z

# Viittaaminen muuttujan havaintoihin: hakasulut

x[1] # ensimmäinen alkio x-vektorista
y[1:5] # viisi ensimmäistä alkiota y-vektorista
z[c(1,5,10)] # ensimmäinen, viides ja kymmenes alkio z-vektorista

# Eräiden perusfunktioiden käyttöä

# keskiarvo
m1 <- mean(x); m1
m2 <- sum(x)/length(x); m2

# keskihajonta
s1 <- sd(x); s1
s2 <-sqrt(var(x)); s2
s3<-sqrt(sum((x-(sum(x)/length(x)))^2)/(length(x)-1)); s3

median(x)
min(x); max(x)

# Aineiston syöttäminen:
data <- data.frame(nimi1=x,nimi2=y)
data

# lasketaan matriisin data rivisummat:
apply(data,1,sum)

# lasketaan matriisin data sarakesummat:
apply(data,2,sum)

# lasketaan sarakkeittain matriisista data keskiarvot ja keskihajonnat:
apply(data,2,mean)
apply(data,2,sd)


# Muuttujan arvojen järjestäminen suuruusjärjestykseen: sort-funktio,
# joka järjestää vektorin alkiot oletusarvoisesti pienimmästä isoimpaan.
# kirjoittamalla "?sort" saat taas help-ikkunan avautumaan

koe <- c(4,5,2,3,1); koe
koe_jarj <- sort(koe) ; koe_jarj
# muutetaan järjestys laskevaksi:
koe_jarj_uusiksi <- sort(koe_jarj, decreasing=TRUE) ; koe_jarj_uusiksi 


# Työskentelyhakemiston määrittäminen ja muuttaminen
# Huom! Käytä hakemistopolussa joko \\ tai / merkkiä.

# näyttää nykyisen työskentelyhakemiston:
getwd()

# voidaan muuttaa työskentelyhakemistoa:
setwd("X:\\")
# Polku voidaan määritellä myös valikosta File/Change Dir...

# oman R-koodin tallentaminen?
# yksinkertainen tapa on tallentaa se notepadiin, ja sieltä sitä voi kopioida suoraan R:ään

# myös koko R-istunnon voi tallentaa haluttuun paikkaan: valikosta File/Save to File