harm_pow: Calculates harmonic power given a list from e/t/rfourier

In vbonhomme/Momocs: Morphometrics using R

Calculates harmonic power given a list from e/t/rfourier

Description

Given a list with an, bn (and eventually cn and dn), returns the harmonic power.

Usage

harm_pow(xf)

Arguments

xf	A list with an, bn (and cn, dn) components, typically from a e/r/tfourier passed on coo_

Value

Returns a vector of harmonic power

Examples

ef <- efourier(bot[1], 24)
rf <- efourier(bot[1], 24)
harm_pow(ef)
harm_pow(rf)

plot(cumsum(harm_pow(ef)[-1]), type='o',
  main='Cumulated harmonic power without the first harmonic',
  ylab='Cumulated harmonic power', xlab='Harmonic rank')

vbonhomme/Momocs documentation built on Nov. 13, 2023, 8:54 p.m.