X=read.table("Orings.txt",header=TRUE)
y=ifelse(X[,2]>0,1,0)
m1=glm(y~X[,1],family="binomial")
logit=function(x) 1/(1+exp(-m1$coef[1]-m1$coef[2]*x))
curve(logit(x),40,90,ylab="P-stwo awarii pierscienia",xlab="Temp")
points(X[,1],y,col=2)
print(summary(m1)) 
print(logit(37.5))