R/CalculateApEn.R

Defines functions CalculateApEn

Documented in CalculateApEn

CalculateApEn <-
function(HRVData, indexNonLinearAnalysis = length(HRVData$NonLinearAnalysis), m=2, tau=1, r=0.2, N=1000, verbose=NULL) {
# -------------------------------------
# Calculates Approximate Entropy
# -------------------------------------
	.Deprecated("CalculateSampleEntropy")
  HRVData = HandleVerboseArgument(HRVData, verbose)
  
	npoints = length(HRVData$Beat$niHR)

	if (npoints > N) {
	  DataInt = HRVData$Beat$niHR[(npoints / 2 - N / 2):(npoints / 2 + N / 2)] 
	}
	else{
	  DataInt = HRVData$Beat$niHR
	}
	r = r*sd(DataInt)
	
	VerboseMessage(HRVData$Verbose, "Calculating Approximate Entropy")
	Phi1 = AvgIntegralCorrelation(HRVData,DataInt,m = m,tau = tau,r = r)
	Phi2 = AvgIntegralCorrelation(HRVData,DataInt,m = (m + 1),tau = tau,r = r)
	ApEn = Phi1 - Phi2
	VerboseMessage(HRVData$Verbose, paste("Approximate Entropy: ", ApEn))
	
	HRVData$NonLinearAnalysis[[indexNonLinearAnalysis]]$ApEn=ApEn
	return(HRVData)
}

Try the RHRV package in your browser

Any scripts or data that you put into this service are public.

RHRV documentation built on Nov. 1, 2022, 1:05 a.m.