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) #当てはめ直線を重ね描き