#' @import fda
#' @import fdakma
#' @importFrom matrixStats rowWeightedMedians rowWeightedMeans
#' @import robustX
#' @import FscanStats
#'
#' @export
###############################################################################
## This function returns a weighted consensus for iterated registration
###############################################################################
getWeightedConsensus <- function( Curves, abscissa ){
medianCurves <- robustX::L1median( t( Curves ) )$estimate
options( warn = -1 )
Sim_All <- FscanStats::getKmaSimilarityWithTemplate(
curves = Curves,
template = medianCurves,
abscissa = abscissa,
similarity.method = c( "d1.pearson" ),
Deriv = FALSE
)
options( warn = 0 )
Weights <- Sim_All
Weights[ Weights < 0 ] <- 0
medianWeighted <- matrixStats::rowWeightedMedians( x = Curves, w = Weights )
meanWeighted <- matrixStats::rowWeightedMeans(x = Curves, w = Weights)
return( list( Mean_Weighted = meanWeighted, Median_Weighted = medianWeighted ) )
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.