Úvodné príkazy

Pomoc

help()

Priradenie premennej

assign("x", c(1,2,3))

Výpis premennej

Práca s vektormi

Priradenie stĺpcového vektora premennej

y<-c(1,2,3,4,5,6)

Transpozícia vektora (príklad vytvorí riadkový vektor)

t(x)

Sčítanie vektorov

x + y
y + x

Vytvorenie vektora z vektorov: vytvorí sa jeden vektor

w<-c(x,0,3,x)

Prístup k prvkom vektora cez []

x[2]

Výpis definovaných premenných

ls()

Odstránenie definovanej premennej rm()

rm(x)

Vektor postupnosti od 1 po 10 vo vektore

x<-c(1:10)

Vektor postupnosti od 2 do 30 s krokom 2

x<-c(2*1:15)

Postupnosť od 1 do 15

1:15

Postupnosť od 1 do 29, krok 2

2*1:15-1

Všeobecné postupnosti cez seq()

seq(from=1, to=10, by=1, length=10)
seq(from=1, to=10, by=1)
seq(from=1, by=1, length = 10)
seq(from=1, by=1, length = 20)

Postupnosť s dĺžkou vektora v parametri along

seq(along=x)

Postupnosť od 3 do 17 s krokom 2

seq(3,17,by=2)
seq(3,17,2)

Postupnosť od 3 do 2 dĺžky 10

seq(3,2,length=10)

Vytvorenie vektora opakovaním iného vektora.

xx<-rep(x, times=5)

Vypísanie minimálneho a maximálneho prvku vektora

range(x5)

Minimum vektora

min(x5)
range(x5)[1]

Maximum vektora

max(x5)
range(x5)[2]

Usporiadanie vektora

sort(x5)

Odmocnina (druhá)

sqrt(pi)

Odmocnina komplexného čísla

sqrt(-17 + 0i)

Operátory

5 >= 3
5 <= 3
5 == 5
5 != 5

Vektor B bude vektor pravdivostných hodnôt, podľa splnenia podmienky

b<-c(1,3,5,7,9)
B <- b>5

Negácia

!B

Špeciálne hodnoty: NaN (Not a Number), Inf (nekonečno)

u=c(2:5, NA, 0, 1, NA, Inf)

Vektor p bude obsahovať tie hodnoty, ktoré nie sú NA

p = u[!is.na(u)]

Vektor p bude obsahovať všetky hodnoty, ktoré nie sú NA a sú > 2

p = u[!is.na(u) & u>2]

Matice

Vytvorenie matice z postupnosti (matica bude vyplnená po stĺpcoch)

m = array(1:20, dim=c(4,5))

Alternatívne vytvorenie matice

matrix(1:20, 5, 4)

Vytvorenie matice po riadkoch

matrix(1:20, 5, 4, TRUE)

Operácie z maticami: skalárne násobenie

2*m1

Sčítanie matíc

m + n

Násobenie matíc po prvkoch (nie štandardné!)

m * n

Násobenie matíc štandardné

m %*% m

Vektorový súčin vektorov

crossprod(x,y)

Počet riadkov vektora

nrow(m)

Počet stĺpcov

ncol(m)

Vytvorenie diagonálnej matice alebo vektora z diagonály matice

diag(v)

Riešenie maticovej sústavy Ax = B:

solve(A, B)

Hľadanie inverznej matice

solve(A)

Hľadanie pseudoinverznej matice pomocou knižnice MASS

library(MASS)
ginv(A)

Vlastné čísla a vlastné vektory matice

eigen(A)

Vlastné čísla matice (A musí byť štvorcová matica)

eigen(A)\$values

Vlastné vektory matice (A musí byť štvorcová matica)

eigen(A)\$vec

Prirodzený logaritmus (log = ln)

log(3)

Dekadický logaritmus

log10(3)

Logaritmus o ľubovoľnom základe

logb(27,3)

Data frames

Načítanie tabuľky zo súboru do štruktúry data.frame

Longley <- read.table("longley1.txt", header=TRUE);

Výpis stĺpca tabuľky (ako vektor)

Longley$NPE

,,Primontovanie" tabuľky -- k stĺpcom bude možno pristupovať pomocou premenných, ktoré prislúchajú názvom stĺpca

attach(Longley)

Výpis stĺpca pomocou premennej

NPE

,,Odmontovanie" tabuľky

detach(Longley)

Po odmontovaní už nie je premenná prístupná (Object "NPE" not found)

NPE

,,Konkatenácia matíc" -- `cbind()` a #rbind()

Príkaz vytvorí trojstĺpcovú maticu s prvým stĺpcom (1,2,3,4), druhým stĺpcom z 2 a tretím stĺpcom z 3

z<-cbind(c(1,2,3,4), 2, 3)

Matica a bude tvorená z 3 stĺpcov: prvý stĺpec bude vektor z postupnosti 1..2, zvyšné stĺpce tvoria maticu m.

m<-matrix(c(1,2,3,4), 2, 2)
a<-cbind(seq(1,l=2, by=0), m)

Výsledok cbind() je stále matica. ,,Typecasting" vektora na maticu (v bude matica)

v <- cbind(c(1,2,3))a

Varovanie, ak ,,skonkatenované výsledné matice" nemajú rovnaký počet riadkov

cbind(1:7, diag(3))

Prvý vektor sa opakuje, aby mal počet riadkov ako druhá matica

cbind(1:3, diag(6))

Použitie rbind() po konkatenovanie ,,po riadkoch"

rbind(c(5, 9, 12), c(1,3,5));

Funkcie

Vlastná funkcia pre výpočet mocniny

nmocnina <- function(n, x) {x^n;}

nmocnica(3,3)

Funkcia pre vytvorenie jednotkovej matice

jednotkovamatica <- function(n) {  v<-seq(1, 1, l = n);   diag(v); }

jednotkovamatica(3)

Lineárne modely

Príkaz lm()

lm(NPE~GNPPD + GNP + UE + AF + NIP14 + YEAR)

Vytiahnutie koeficientov modelu coef() -- a je zoznam

a <- coef(lm(NPE~GNPPD + GNP + UE + AF + NIP14 + YEAR))

Prístup k hodnotánm

a["GNPPD"]

Štatistické znaky súboru

Štatistické veličiny súboru summary()

summary(a)
summary(Longley)
summary(lm(NPE~GNPPD + GNP + UE + AF + NIP14 + YEAR))

Načítanie data setu -- data()

data(faithful)

Výpis

faithful

Primontovanie

attach(faithful)

Výpis stĺpcov cez premenné

eruptions
waiting

Vykreslenie jednoduchého histogramu. | udáva pozíciu desatinnej čiarky.

stem(eruptions)
stem(waiting)