Displaying 1 result from an estimated 1 matches for "ticksatlog".
2007 May 27
0
weibplot (Weibull plot) for R
...quot;x",axes=FALSE,frame.plot=TRUE,xlab="data",ylab="probability")
# Let the confidence interval code below plot the regression curve
# curve(log(-log(exp(-(x/eta)^beta))),log="x",add=TRUE)
# annotate the graph
ticklabels=c(.1,.3,.63,.9)
ticksat=log(-log(1-ticksatlog))
axis(2,at=ticksat,labels=ticklabels)
axis(1)
grid()
# plot 95% confidence intervals
x<-data
out=lm(log(-log(1-ppoints(n,a=.5)))~log(x))
curve(predict(out,newdata=data.frame(x=x)),add=TRUE)
curve(predict(out,newdata=data.frame(x=x),level=confidence_level,interval="confidence&quo...