attach(kyphosis)

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

x=as.numeric(dimnames(freq)[[1]])
y=freq[,1]/total  #手術成功率

plot(x, y, xlab="", ylab="") #手術成功率をプロット

lmSuccess=lm(y[1:6]~x[1:6]) # 手術脊椎数が9や10のケースをのぞいて，手術成功率に線形モデルを当てはめる
abline(lmSuccess, lty=3) #当てはめ直線を重ね描き