Estimate parameters of the normal distribution accounting for censoring

norm_M3.Rd

Mean and standard deviation of a normal distribution are estimated based on maximum likelihood estimation. For samples above the LLOQ, the likelihood is determined from the normal density dnorm(x, mean, sd). For values below the LLOQ, the likelihood is the integrated density pnorm(x, mean, sd). Mean and standard deviation are estimated by optimization of -2 LL.

norm_M3(x, LLOQ = -Inf, return_LLOQ = -Inf)

Arguments

x

numeric vector, the normally distributed sample. Censored values should be coded by numbers <= LLOQ.

LLOQ

value or numeric valus of length of x indicating the lower limit of quantification (left censoring).

return_LLOQ

numeric, the return value for the case that all elements of x are below LLOQ.

Value

List with mean (mean) and standard deviation (sd).

Details

The return values for mean and standard deviation are the estimated mean and N/(N-1) times the estimated standard deviation, where N is the length of the sample. I.e., the Bessel corrected standard deviation is returned.

See also

Other Auxiliary: IQRloadCSVdata(), IQRsaveCSVdata(), and(), aux_explode(), aux_explodePC(), aux_fileparts(), aux_fileread(), aux_filewrite(), aux_getRelPath(), aux_mkdir(), aux_na_locf(), aux_postFillChar(), aux_preFillChar(), aux_quantilenumber(), aux_rmdir(), aux_simplifypath(), aux_splitVectorEqualPieces(), aux_strFindAll(), aux_strrep(), aux_strtrim(), aux_unlevel(), aux_version(), calcAICBIC(), clusterX(), compare_IQRmodel_IQRsysModel_simulation(), fit_EmaxModel(), format_GUM(), ge(), gen_aux_version(), geocv(), geomean(), geosd(), ginv(), gt(), interp0(), interp1(), interpcs(), inv_logit(), le(), logit(), lt(), mod(), mvrnorm(), or(), piecewise(), progressBar(), remove_duplicates(), run_silent_IQR(), stopIQR(), tempdirIQR(), tempfileIQR(), warningIQR()

Examples

if (FALSE) { # \dontrun{
norm_M3(rnorm(10000), LLOQ = -Inf, return_LLOQ = -Inf)
norm_M3(rnorm(10000), LLOQ = 0, return_LLOQ = -Inf)
x <- rnorm(10000)
x[x<0] <- 0
norm_M3(x, LLOQ = 0, return_LLOQ = -Inf)
} # }