R/MRRln3p.r

Defines functions MRRln3p

Documented in MRRln3p

## MRRln3p.r file
 ##
## A quick fit method using defaults, with options to provide confidence interval bounds
## and/or a graphical display.



MRRln3p<-function(x, s=NULL, bounds=FALSE, show=FALSE)  {
# permit one to forget that second argument was reserved for a suspensions vector
if(is.logical(s)) {
	show<-bounds
	bounds<-s
	s<-NULL
}

obj<-wblr.fit(wblr(x,s), dist="lognormal3p", modify.by.t0=T, col="red2")
fit<-obj$fit[[1]]$fit_vec
attributes(fit)$data_types<-NULL
## bounds are not prepared for 3p fits, so the bounds argument is simply ignored.
#if(bounds==TRUE) {
#obj<-wblr.conf(obj, dq="minitab", col="pink2")
#bnds<-obj$fit[[1]]$conf[[1]]$bounds
#ret<-list(fit,bnds)
#stitle<-"MRR fit with 90% double-sided pivotal bounds"
#}else{
ret<-fit
stitle<-"Lognormal MRR 3-parameter fit"
#}

if(show==TRUE) {
plot(obj, canvas="lognormal", in.legend.blives=FALSE, sub=stitle)
}

ret
}

Try the WeibullR package in your browser

Any scripts or data that you put into this service are public.

WeibullR documentation built on July 27, 2018, 3 p.m.