# 
# Joonte lisamine - erinevad joonetüübid
#
# 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 - kolm firmat
#
# ---------------------------------------------------


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

par(mar=c(3,3,4,2))

#Algandmed:

aasta=1999:2005
Firma1=c(12,14,18, 24,25,26,24)
Firma2=c(20, 21, 22, 20, 21, 25, 27)
Firma3=c(NA,0,5, 15, 25, 30, 32)

# Alustame graafiku joonistamist - paneme paika teljed ja kanname joonisele
#	1. firma andmed:
plot(aasta, Firma1, type="l", lty=1, ylim=c(0,35), main="Meie konkurendid")

# Lisame joonisele 2. ja 3. firma andmed:
lines(aasta, Firma2, lty=2)
lines(aasta, Firma3, lty=3)

# lisame joonetüüpide selgituse
legend(2003, 15, c("Firma 1","Firma 2", "Firma 3"), lty=1:3)

dev.off()
system(teenimi2)





# ---------------------------------------------------
#
# 	Teine pilt - Usalduspiirid
#
# ---------------------------------------------------


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

par(mar=c(4,4,4,2))

#Algandmed:

x=c(1, 3, 4, 10, 11, 13, 13)
y=c(3, 5, 7, 13, 10, 14, 17)

# Mudeli hindamine

m1=lm(y~x)
summary(m1)

# Leiame prognoosid x-i väärtuste x=0, x=0.1, x=0,2,... x=20 jaoks:
abi=predict(m1, data.frame(x=seq(0,20,0.1)), type="response", interval="confidence", level=0.95)

# Teeme esimese joonise - meie vaatlused
plot(x,y, main="Hajuvusgraafik ja regressioonisirge")

# Lisame regressioonisirge ja 95%-usaldusintervalli regressioonisirgele:
lines(seq(0,20,0.1), abi[,1], lty=1)
lines(seq(0,20,0.1), abi[,2], lty=2)
lines(seq(0,20,0.1), abi[,3], lty=2)

dev.off()
system(teenimi2)


# ---------------------------------------------------
#
# 	Kolmas pilt - joonetüübid
#
# ---------------------------------------------------


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


par(mar=c(1,1,1,1))
plot(NA,NA, xlab="",ylab="", main="", xaxt="n", yaxt="n", xlim=c(0,10), ylim=c(0,7))
for (i in 1:6){
lines(c(2,10), c(7-i,7-i), lty=i)
text(1, 7-i, paste("lty=",i, sep=""), adj=c(0.5,0.5), cex=1.2)
}

dev.off()
system(teenimi2)