# 
# Olemasolevale graafikule elementide lisamine
#
# 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")
}

# ---------------------------------------------------
#
# 	Joonte lisamine -- põhitõed 
#
# ---------------------------------------------------

algus("pilt11.png", 300, 225)

par(mar=c(2.5,2.5,3.5,1))

# plot(NA,NA, xlab="", ylab="", xaxt="n", yaxt="n", main="Põhitõed", xlim=c(0,5),ylim=c(0,4))

plot(NA,NA, xlab="", ylab="", main="Põhitõed", xlim=c(0,5),ylim=c(0,4))
grid()
lines(c(1,3), c(1,2), lwd=2)
# text(2.5,3, "Põhitõed", adj=c(0.5,0.5), cex=1.5)
dev.off()
system(teenimi2)


# ---------------------------------------------------
#
# 	Joonte lisamine -- erinevad joonetüübid 
#
# ---------------------------------------------------


algus("pilt12.png", 300, 225)

x=runif(100)
y=1+1*x+2*x*x+rnorm(100)*0.4
par(mar=c(1,1,3.5,1))
plot(x,y, xlab="", ylab="", main="Joonetüübid", xaxt="n",yaxt="n", col="gray70")

m1=lm(y~x+I(x^2))
xuus=seq(-0.1,1.1,0.01)
abi=predict(m1, data.frame(x=xuus), type="response", interval="prediction", level=0.90)

# Mudeli prognoos
lines(xuus, abi[,1]) 

# 90%-prognoosiintervall punktiirjoonega - lty=2;
lines(xuus, abi[,2], lty=2) 
lines(xuus, abi[,3], lty=2) 

# 99%-prognoosiintervall punktiirjoonega - lty=2;
abi=predict(m1, data.frame(x=xuus), type="response", interval="prediction", level=0.99)
lines(xuus, abi[,2], lty=3) 
lines(xuus, abi[,3], lty=3) 

legend(0.65, 2.1, c("Prognoos", "90%-PI", "95%-PI"), lty=c(1,2,3), cex=0.75)

dev.off()
system(teenimi2)



# ---------------------------------------------------
#
# 	Joonte lisamine -- joone värvist
#
# ---------------------------------------------------


algus("pilt13.png", 300, 225)

YTalvel=c(1,2,4,4.5,5.2,5.5,8)
YSuvel=c(7,6.9,6.5,4,3.5,3,4)
YSygis=c(7.5,7.3,7.1,7,7,7.2,8.5)
TunnusX=0:6

par(mar=c(1,1,3.5,1))
plot(TunnusX, YTalvel, type="l", col="lightblue1", lwd=3, xlab="", ylab="", main="Joone värv", xlim=c(0,6), ylim=c(0,9), xaxt="n", yaxt="n")
grid()

lines(TunnusX, YSuvel, col="green3", lwd=3)
lines(TunnusX, YSygis, col="gold", lwd=3)

legend(4.3,3, c("Suvi","Sügis","Talv"), col=c("green3","gold","lightblue1"), lty=1, lwd=2, cex=0.8)


dev.off()
system(teenimi2)


# ---------------------------------------------------
#
# 	Joonte lisamine -- Muud trikid
#
# ---------------------------------------------------

algus("pilt14.png", 300, 225)

par(mar=c(2,2,3.5,1))
Konkurent=c(115,100, 90, 100, 75,60)
Meie=c(25, 40, 80, 90, 100, 150)
Aasta=2000:2005
plot(Aasta, Konkurent, type="l", bty="u", ylim=c(0,120), yaxs="i", xaxs="i", main="Muud trikid", xlim=c(2000,2005))

par(xpd=NA)
lines(Aasta,Meie, lwd=4)
lines(c(2000,2004.25, NA, 2004.55,2005),c(120,120,NA,120,120))

text(2001,20, "Meie", cex=1.2)
text(2004.25,60, "Nemad")

dev.off()
system(teenimi2)


# ---------------------------------------------------
#
# 	Joonte lisamine -- Tüüpnäited 
#
# ---------------------------------------------------

algus("pilt15.png", 300, 225)

par(mar=c(2.5,2.5,3.5,1))

x=rnorm(40)

hist(x, probability=TRUE, xlim=c(-4,4), main="Tüüpolukordi", col="gold", border="darkgoldenrod2")
xkeskmine=mean(x)
xshalve=sd(x)

xabi=seq(-4,4,0.1)

lines(xabi, dnorm(xabi, mean=xkeskmine, sd=xshalve), lwd=3, col="darkgoldenrod3")

dev.off()
system(teenimi2)