attach(kyphosis)

freq=table(Number, Kyphosis)
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) #当てはめ直線を重ね描き