# 
# Joonte lisamine - tüüpolukordi
#
# Viimati muudetud - 28. august 2005
# Koostanud: Märt Möls
#


# Millega png-pilte vaadata
vaataja="C:\\Program Files\\Mozilla Firefox\\firefox.exe"

# Kuhu pildfailid salvestada
tee="F:\\kursus\\R\\temp\\lisamine\\"

algus=function(pilt, xs, ys){
assign("teenimi", paste(tee,pilt,sep=""), env=.GlobalEnv)
assign("teenimi2", paste('"',vaataja,'" ',tee,pilt,sep=""), env=.GlobalEnv)
png(filename = teenimi, width = xs, height = ys, pointsize = 12, bg = "white")
}

# ---------------------------------------------------
#
# 	Esimene pilt -- Teoreetiline tihedus & histogramm
#
# ---------------------------------------------------

algus("pilt50.png", 400, 300)

# Algandmete genereerimine
x=rnorm(60)

hist(x, probability=TRUE, xlim=c(-4,4), ylim=c(0,0.4),
	main="Tihedusfunktsioon ja histogramm", 
	ylab="tihedus",
	col="gold", border="darkgoldenrod2")

xkeskmine=mean(x)
xshalve=sd(x)

x=seq(-4,4,0.1)
y=dnorm(xabi, mean=xkeskmine, sd=xshalve)

lines(x,y, lwd=3, col="darkgoldenrod3")

dev.off()
system(teenimi2)



# ---------------------------------------------------
#
# 	Teine pilt -- Hinnatud tihedus & histogramm
#
# ---------------------------------------------------

algus("pilt51.png", 400, 300)

# Algandmete genereerimine
x=runif(60)

hist(x, probability=TRUE, xlim=c(-0.5,1.5), ylim=c(0,1.5),
	main="Hinnatud tihedusfunktsioon", 
	ylab="tihedus",
	col="gold", border="darkgoldenrod2")

lines(density(x), lwd=3, col="darkgoldenrod3")

dev.off()
system(teenimi2)





# ---------------------------------------------------
#
# 	Kolmas pilt -- usaldatavad ja mitteusaldatavad andmed
#
# ---------------------------------------------------

algus("pilt52.png", 400, 300)

# Algandmete genereerimine
x=runif(60)*2
y=2+3*x+rnorm(60)

plot(x,y, col="gray80", xlim=c(-1,3), ylim=c(-5,15))

m1=lm(y~x)
xuus=seq(0,2,0.01)
yuus=predict(m1, data.frame(x=xuus))
lines(xuus, yuus, lwd=2)

xuus=seq(-3,5,0.01)
yuus=predict(m1, data.frame(x=xuus))
lines(xuus, yuus, lty=2, lwd=2)

dev.off()
system(teenimi2)





# ---------------------------------------------------
#
# 	Neljas pilt -- Silumine
#
# ---------------------------------------------------

algus("pilt53.png", 400, 300)

x=runif(200)*2
y=2+3*x+4*sin(2*x)+1.5*sqrt(x)+rnorm(200)

plot(x,y, col="gray75", xlim=c(0,2), ylim=c(0,15), main="Silujaid")

# Esimene siluja
lines(lowess(x,y), lwd=2)

# Veidi tundlikum joon (mida väiksem f-väärtus, seda hüplikum on siluja joon)
lines(lowess(x,y, f=0.25), lty=2, lwd=2)

# Teisel tööpõhimõttel töötav siluja
lines(ksmooth(x,y), col="orange", lwd=2)

dev.off()
system(teenimi2)




# ---------------------------------------------------
#
# 	Viies pilt -- qq-plot e. tõenäosuspaber
#
# ---------------------------------------------------

algus("pilt54.png", 400, 300)

x=rnorm(100)
qqnorm(x, xlab="teoreetilised kvantiilid", ylab="valimikvantiilid", main="Tõenäosuspaber (QQ-plot)")
qqline(x, lwd=2)

dev.off()
system(teenimi2)