attach(kyphosis)

freq=table(Start, Kyphosis) # FALSE （手術成功）と TRUE（手術失敗）の度数からなる各行
total=apply(freq, 1, sum)

plot(as.numeric(dimnames(freq)[[1]]), freq[,1]/total, xlab="", ylab="") #手術成功率のプロット

glmStart=glm(Kyphosis ~ Start, family = binomial) #症状が残ったかどうかの論理データベクトルに対して glm を当てはめ

logistic=function(x, coef){
	y=cbind(1, x) %*% coef
	exp(y)/(1+exp(y))
	} #逆リンク関数

x=1:18
lines(x, 1-logistic(x, glmStart$coef)) #当てはめ曲線の重ね描き