The y-bounded r-th raw moment of the distribution equals:
μyr=scalerΓ(shaper+1,(scaley)shape)
where Γ(,) denotes the upper incomplete gamma function.
Examples
## The zeroth truncated moment is equivalent to the probability functionpweibull(2, shape =2, scale =1)mweibull(truncation =2)## The (truncated) first moment is equivalent to the mean of a (truncated) random sample,#for large enough samples.x <- rweibull(1e5, shape =2, scale =1)mean(x)mweibull(r =1, lower.tail =FALSE)sum(x[x > quantile(x,0.1)])/ length(x)mweibull(r =1, truncation = quantile(x,0.1), lower.tail =FALSE)