x = c(10, 8.65, 8.90, 8.01) #havaintojen x-koordinaatit vektorissa x
y = c(8.65, 8.90, 8.01, 7.81) # havaintojen y-koordinaatit vektorissa y

fit = lm(y~x) #luo objektin fit, jossa on lineaarisen mallin sovite kun vakio on mukana mallissa
fit = lm(y~x-1) # ... vakio ei mukana mallissa
summary(fit) # antaa oleellisimmat tiedot mallista: selitysaste, parametrien arvot, p-arvot, jne

plot(x,y, xlim=c(-5,5), ylim=c(-5,5)) # piirtŠŠ havaintopisteet x,y -koordinaatistoon. ylim ja xlim ovat ylimŠŠrŠisiŠ parametreja, jotka kertovat minkŠ alueen koordinaatistosta r piirtŠŠ.
abline(fit) #piirtŠŠ sovitteen (auki olevaan) kuvaan
points(k,l, pch=15, col="red") #piirtŠŠ (auki olevaan) kuvaan pisteet havainnot (k,l) (mŠŠritellŠŠn samalla tavalla kuin x ja y), chp ja col ylimŠŠrŠisiŠ parametreja, jotka kertovat
					ĘĘhavaintopisteen muodon ja vŠrin

t(x) #x:n transpoosi
X%*%Y #matriisikertolasku
a*b #tavallinen reaalilukujen kertolasku
inv(t(X)$*$X) # = (X«X)^(-1)
sum(x) #vektorin x alkioiden summa
mean(x) # vektorin x alkioiden kekiarvo