M5two: Pooled central moment estimates - two-sample

In innager/edgee: Edgeworth Expansions and Higher-Order Inference

Pooled central moment estimates - two-sample

Description

Calculate pooled unbiased estimates of central moments and their powers and products.

Usage

M5two(m2, m3, m5, n_x, n_y)

Arguments

m2	naive biased variance estimate m_2 = 1/(n_x + n_y) \sum_{i = 1}^{n_x} ((X_i - \bar{X})^2 + \sum_{i = 1}^{n_y} ((Y_i - \bar{Y})^2 for vectors X and Y.
m3	naive biased third central moment estimate m_3 = 1/(n_x + n_y) \sum_{i = 1}^{n_x} ((X_i - \bar{X})^3 + \sum_{i = 1}^{n_y} ((Y_i - \bar{Y})^3 for vectors X and Y.
m5	naive biased fifth central moment estimate m_5 = 1/(n_x + n_y) \sum_{i = 1}^{n_x} ((X_i - \bar{X})^5 + \sum_{i = 1}^{n_y} ((Y_i - \bar{Y})^5 for vectors X and Y.
n_x	number of observations in the first group.
n_y	number of observations in the second group.

Value

Pooled estimate of a fifth central moment.

See Also

Other unbiased estimates: M2M3one(), M2M3two(), M2M4one(), M2M4two(), M2one(), M2pow2one(), M2pow2two(), M2pow3one(), M2pow3two(), M2two(), M3one(), M3pow2one(), M3pow2two(), M3two(), M4one(), M4two(), M5one(), M6one(), M6two()

Examples

n1 <- 10
n2 <- 8
shp <- 3
smp1 <- rgamma(n1, shape = shp) - shp
smp2 <- rgamma(n2, shape = shp)
for (j in 2:6) {
  assign(paste("m", j, sep = ""), 
         mean(c((smp1 - mean(smp1))^j, (smp2 - mean(smp2))^j)))
}          
M5two(m2, m3, m5, n1, n2)

innager/edgee documentation built on April 24, 2024, 8:14 p.m.