R/clnorm.R

In musto101/wilcox_R: What the package does (short line)

clnorm <-
function(n,epsilon=.1,k=10){
#
# generate n observations from a contaminated lognormal
# distribution
#
#  Using default values, median is approximately 1.14 and 20% trimmed mean is 1.33
if(epsilon>1)stop('epsilon must be less than or equal to 1')
if(epsilon<0)stop('epsilon must be greater than or equal to 0')
if(k<=0)stop('k must be greater than 0')
val<-rlnorm(n)
uval<-runif(n)
flag<-(uval<=1-epsilon)
val[!flag]<-k*val[!flag]
val
}