powered by
Builds a list containing CUPM, DUPM, DLPM, CLPM and the overall covariance matrix.
PM.matrix(LPM_degree, UPM_degree, target, variable, pop_adj, norm = FALSE)
A list: $cupm, $dupm, $dlpm, $clpm, $cov.matrix.
numeric; lower partial moment degree (0 = freq, 1 = area).
numeric; upper partial moment degree (0 = freq, 1 = area).
numeric vector; thresholds for each column (defaults to colMeans).
numeric matrix or data.frame.
logical; TRUE adjusts population vs. sample moments.
logical; default FALSE. If TRUE, each quadrant matrix is cell-wise normalized so their sum is 1 at each (i,j).
Partial Moment Matrix
set.seed(123) A <- cbind(rnorm(100), rnorm(100), rnorm(100)) PM.matrix(1, 1, NULL, A, TRUE) # uses norm = FALSE by default PM.matrix(1, 1, NULL, A, TRUE, TRUE) # enable normalization
Run the code above in your browser using DataLab