https://github.com/cran/ftsa
Raw File
Tip revision: 592f61af0e8cfd418ea2dfc4fb50aaa23bd972cc authored by Han Lin Shang on 03 January 2019, 07:30:06 UTC
version 5.3
Tip revision: 592f61a
metri.p.R
metri.p <- function(x1, x2, p = 2, w = 1)
{
           if (length(x1) != length(x2)) 
               stop("Dimensions of matrices do not match\n")
           if (length(w) == 1) 
               w = rep(w,length(x1))
           if (length(w) != length(x1)) 
               stop("Length of the vector does not match with the matrix\n")
           if (p > 0){
               (sum(w * (abs(x1 - x2) ^ p))) ^ (1/p)
           } 
           else {
               max(w * abs(x1 - x2))
           }
}
back to top