# Kontingenssitaulukkojen Poissonin regressio (Altham s. 42)
#
# Käytetään Poissonin jakaumaa kaksiulotteisen kontingenssitaulukon
# logaritmisesti lineaariseen analysointiin ja verrataan tulosta vastaavaan
# binomiaaliseen formulointiin.
#   
#  The Daily Telegraph (28.10.1988) julkaisi artikkelin
# "Johtajat alttiita ajamaan humalassa", jossa esitettiin 
# Royal Ascotin ja Henley Regatan yhteydessä järjestettyjen
# puhallustestien tuloksia (Royal Ascot ja Henley Regatta ovat
# tunnetusti "kosteita" englantilaisia urheilutapahtumia kuten
# hevos- ja soutukilpailut). Jääkö RA:ssa todennäköisemmin
# kiinni kuin HR:ssä (jos testataan)? RA:ssa puhallutettiin
# kaikkiaan 24+2210 ja HR:ssä 5+680. Näistä 24 positiivisia
# RA:ssa ja 5 HR:ssä.
# 
r <- scan()
24 2210   
 5 680    

# 1. rivi Royal Ascot
# 2. rivi Henley Regatta

Row <- c(1,1,2,2); Col <- c(1,2,1,2)
# Älä käytä nimeä "row", koska se on funktion nimi R:ssä
ROW <- factor(Row); COL <- factor(Col)

# Col=1 PIDÄTETTY; Col=2 EI pidätetty