# TYÖARKKI 7.
# 
# Logistinen regressio ja turvallisuus avaruudessa.
# Swan ja Rigby (GLIM Newsletter 24, 1995) käsittelevät oheista aineistoa
# käyttäen binomiaalista logistista regressiota.
# Lainaus artikkelista: 
# "In 1986 the NASA space shuttle Challenger exploded shortly after it 
# was launched. After an investigation ...
#  ...
# The temperatures for the 23 previous launches ranged from 53 to 81 degrees
# Fahrenheit (ie, -0.6 degrees Centigrade)"
#
# 
data <- read.table(("A:\\Kurssit\\TMallit\\Datat\\Shuttle.txt"),header=T)
six <- rep(6,times=23)	# "O" renkaiden lkm 6 kaikissa laukaisuissa
attach(data)
pfail <- nfail/six		# pettäneiden renkaiden suht. osuus	
pfail			
plot(temp,pfail)	# vastaava lämpötila
plot(temp,pfail,type="l")
summary(data)
sh <- glm(pfail ~ temp,binomial,weights=six)
sh; summary(sh)	# Kommentoi tuloksia
# Piirrä sopivia kuvioita
names(sh)
plot(sh,ask=T)
# Hit <Returnto see next plot: 
# Hit <Returnto see next plot: 
#
# Ovatko jotkut havainnot erityisen vaikuttvia logist. regressiossa?
# Mikä vaikutus tuloksiin on kaikkien sell. hav. Poistaminen,
# joilla nfail=0. 
sh0 <- glm(pfail ~ temp,binomial,weights=six,subset=(nfail>0))
plot(temp,pfail)
sh0
summary(sh0)
plot(sh0,ask=T)
# Hit <Returnto see next plot: 
# Hit <Returnto see next plot:
# Poista joku vaikuttava havainto ja tarkastele tuloksia.
# Tutki joidenkin kahden vaikuttan poistamisen vaikutusta. 
sh1 <- glm(pfail ~ temp,binomial,weights=six,subset=(nfail != 2 | temp !=75))
sh2 <- glm(pfail ~ temp,binomial,weights=six,subset=(nfail != 2))
sh0; sh; sh1; sh2
summary(sh0); summary(sh); summary(sh1); summary(sh2)