seewave (version 1.5.5)

specprop: Spectral properties

Description

This function returns a list of statistical properties of a frequency spectrum.

Usage

specprop(spec, f, str = FALSE, flim=NULL, plot = FALSE, type = "l", ...)

Arguments

spec

a data set resulting of a spectral analysis obtained with spec or meanspec (not in dB).

f

sampling frequency of spec (in Hz).

str

logical, if TRUE returns the results in a structured table.

flim

a vector of length 2 to specifgy the frequency limits of the analysis (in kHz)

plot

if 1 returns the spectrum , if 2 returns the cumulative spectrum, both of them with the first quartile, the third quartile, the median and the mode plotted (by default FALSE).

type

if plot is TRUE, type of plot that should be drawn. See plot for details (by default "l" for lines).

...

other plot graphical parameters.

Value

• A list of 15 values is returned
• mean
  mean frequency (see mean)
• sd
  standard deviation of the mean (see sd)
• sem
  standard error of the mean
• median
  median frequency (see median)
• mad
  absolute deviation of median (see mad)
• mode
  mode frequency, i.e. the dominant frequency
• Q25
  first quartile (see quantile)
• Q75
  third quartile (see quantile)
• IQR
  interquartile range (see IQR)
• cent
  centro�d, see note
• skewness
  skewness, a measure of asymmetry, see note
• kurtosis
  kurtosis, a measure of peakedness, see note
• sfm
  spectral flatness measure (see sfm)
• sh
  spectral entropy (see sh)
• prec
  frequency precision of the spectrum

Details

The spectrum is converted in a probability mass function (PMF). If a selected value has to be selected with $, the argument str has to be set to FALSE.

Examples

data(orni)
a<-meanspec(orni,f=22050,plot=FALSE)
specprop(a,f=22050)
# to get a single measure of the list
specprop(a,f=22050)$mode
# to get the results structured
specprop(a,f=22050,str=TRUE)
# to limit the analysis between 4 and 6 kHz
specprop(a,f=22050,flim=c(4,6),str=TRUE)
# plots
specprop(a,f=22050,plot=1)
specprop(a,f=22050,plot=2)