Yi<-c(2.1,2.5,4.9,5.5,7,8.4,9.6,10.2,11.4,12.5,13.1,14.6,17,16.8,18.6,19.7,21.3,21.6)
Xi<-c(1,1.5,2,3,4,5,6,7.5,8.5,10,12.5,15,17.5,20,25,30,35,40)
Y<-1/Yi
X<-1/Xi
lm(Y~X)
nls(Yi~(a*Xi)/(b+Xi),start=list(a=29.62085,b=13.44816),trace=TRUE)
plot(Xi,Yi)
curve((28.14*x)/(12.57+x),add=TRUE,0,40)
summary( nls(Yi~(a*Xi)/(b+Xi),start=list(a=29.62085,b=13.44816),trace=TRUE))
res<-resid(nls(Yi~(a*Xi)/(b+Xi),start=list(a=29.62085,b=13.44816),trace=TRUE))
fit<-fitted(nls(Yi~(a*Xi)/(b+Xi),start=list(a=29.62085,b=13.44816),trace=TRUE))
plot(fit, res)
plot(Xi, res)
qqnorm(res)
confint(nls(Yi~(a*Xi)/(b+Xi),start=list(a=29.62085,b=13.44816),trace=TRUE))

Yi<-c(0,0,1,0,0,1,1,1,0,1,1,0,0,1,0,0,1,0,1,0,0,1,1,0,0,1,0,0,1,0,0,1,0)
Xi1<-c(32,45,60,53,25,68,82,38,67,92,72,21,26,40,33,45,61,16,18,22,27,35,40,10,24,15,23,19,22,61,21,32,17)
Xi2<-c(3,2,2,1,4,1,2,5,2,2,3,5,3,4,3,1,2,3,4,6,3,3,3,4,3,4,3,5,2,2,3,5,1)
mod1<-lm(Yi~Xi1+Xi2)
mod2<-glm(Yi~Xi1+Xi2,family=binomial)
mod2
sum((Yi-33*fitted(mod2))^2 / 33*fitted(mod2)*(1-fitted(mod2)))
qchisq(.99,30)
res<-resid(mod2)
plot(c(1:33),res)
anova(mod2, mod1,test="Chisq")

Xi<-c(1, 0, 2, 0, 3, 1, 0, 1, 2, 0)
Yi<-c(16, 9, 17, 12, 22, 13, 8, 15, 19, 11)
plot(Xi,Yi)
glm(Yi~Xi,poisson)
anova(glm(Yi~Xi,poisson))
summary(glm(Yi~Xi,poisson))
deviance(glm(Yi~Xi,poisson))
sum(((Yi-fitted(glm(Yi~Xi,poisson))) / fitted(glm(Yi~Xi,poisson))^.5 )^2)
res<-resid(glm(Yi~Xi,poisson))
res
plot(c(1:10),res)
plot(glm(Yi~Xi,poisson))
predict(glm(Yi~Xi,poisson), data.frame(Xi=0), type="response")
mod2<-lm(Yi~Xi)
predict(lm(Yi~Xi), data.frame(Xi=0), type="response")

x<-c(1, 0, 2, 0, 3, 1, 0, 1, 2, 0)
y<-c(16, 9, 17, 12, 22, 13, 8, 15, 19, 11)
plot(x,y)
abline(lm(y~x))
curve(exp(0.2638*x+2.3529),add=TRUE,col="red",0,3)
predict(glm(Yi~Xi,poisson), data.frame(Xi=0), type="response")

x<-c(0,0,0,0)
y<-c(9,12,8,11)
glm(y~x, poisson)
confint(glm(Yi~Xi,poisson))