# 
# mar, xpd graafikaparameetrite võimaluste demonstreerimine
#
# Viimati muudetud - 19. juulil 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\\komb\\"

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")
}



# Järgnev osa joonistab pildi "jooniste täiendamise" lehekülje tarvis - logistiline regr.

algus("pilt30.png", 600, 450)

# Algandmete genereerimine
tootlus=rep(6:12, 8)
onn=rbinom(7*8, 1, (tootlus-5)/8)

# Hindame logistilise regressiooni mudeli
mudel=glm(onn~tootlus, family=binomial)

# Kasutades mudelit arvutame prognoosid
x1=seq(6,12,0.1)
x2=c(seq(4,16,0.1))
y1=predict(mudel, data.frame(tootlus=x1), type="response")
y2=predict(mudel, data.frame(tootlus=x2), type="response")


# Asume joonist joonistama
windows(width=6, height=4.5)

# Jätame joonise kohale vaba ruumi, hetkel lülitame ka pildist väljapoole jäävate
# joonte ja punktide "äralõikamise" reziimi sisse (xpd=FALSE)

par(mar=c(5,5,7,2)+0.1, xpd=FALSE)

plot(x1,y1*100, type="l", xlim=c(4,16), ylim=c(0,100), bty="n", xlab="Pingo-Pango concentration (mg)",ylab="Probability of success (%)", las=1, yaxs="i",lwd=2)

lines(x2,100*y2, lty=2, lwd=2)

# Nüüd lubame joonistada ka väljapoole joonistuspiirkonda:
par(xpd=NA)

koht=120
lines(c(4,16),c(koht,koht))
text(15.7,koht+2,"Succesful", adj=c(1,0), cex=1)
text(15.7,koht-2,"Unsuccesful", adj=c(1,1), cex=1)

aa=table(tootlus[onn==1])
xx=as.numeric(names(aa))
a=function(x){return(seq(1,x))}
yy=c(lapply(aa,a), recursive=TRUE)
points(rep(xx,aa),yy*1.5+koht, pch=20)

aa=table(tootlus[onn==0])
xx=as.numeric(names(aa))
yy=c(lapply(aa,a), recursive=TRUE)
points(rep(xx,aa),-yy*1.5+koht, pch=20) 

dev.off()
system(teenimi2)