censored Weibull MLE for Type II right-censored data

Description

Calculate the censored Weibull MLE for the Type II censored data with the algorithm described in ASTM 5457(2004). Return the estimates of the shape and scale parameters. A comprehesive description of this algorithm can be found in Liu (2012).

Usage

1
cenWbMLE.T2(dat, n, useC = FALSE, conCr = 1e-09, nIter = 1000)

Arguments

dat

A vector of the observations after censoring. None NA or negative values allowed.

n

The original sample size, including both the censored and uncensored observations.

useC

See cenWbMLE.T1

conCr

See cenWbMLE.T1

nIter

See cenWbMLE.T1

Value

See cenWbMLE.T1

Note

Please report the numerical problems and inconvenience when using this function to the author.

Please notice that the ways of inputing data in cenWbMLE.T1 and cenWbMLE.T2 are differnt. For cenWbMLE.T1, the algorithm require a full "orginal" data set (with the uncensored observations as NA or a arbitary value larger than the threshold) and the original sample size is decided as the length of the input data, while cenWbMLE.T2 requires the observed data points and the original sample size.

Author(s)

Yang (Seagle) Liu <yang.liu@stat.ubc.ca>

References

See cenWbMLE.T1

See Also

rweibull, cenWbMLE.T1, emCenWbMix.T2

Examples

1
2
3
4
set.seed(1)
y <- sort(rweibull(100, 7, 7)) ##Generate the data
cenWbMLE.T2(y, 100) #The MLE for the complete data
cenWbMLE.T2(y[1:10], 100) #Censor the largerst 90% of the data.