integratePower {astrochron} R Documentation

## Determine the total power within a given bandwidth

### Description

Determine the total power within a given bandwidth, and also the ratio of this power to the total power in the spectrum (or up to a specified frequency). If bandwidth is not specified, generate interactive plots for bandwidth selection. For use with the function eha, integratePower can process spectrograms (time-frequency) or single spectra.

### Usage

integratePower(spec,flow=NULL,fhigh=NULL,fmax=NULL,unity=F,f0=T,xmin=NULL,
palette=6,ncolors=100,h=6,w=9,ln=F,genplot=T,verbose=T)


### Arguments

 spec Spectral results to evaluate. If the data frame contains time-frequency results, it must have the following format: column 1=frequency; remaining columns (2 to n)=power; titles for columns 2 to n must be the location (depth or height). Note that this format is ouput by function eha. If the data frame contains one spectrum, it must have the following format: column 1=frequency, column 2=power. flow Low frequency cutoff for integration. If flow or fhigh are not specified, interactive plotting is activated. fhigh High frequency cutoff for integration. If flow or fhigh are not specified, interactive plotting is activated. fmax Integrate total power up to this frequency. unity Normalize spectra such that total variance (up to fmax) is unity. (T of F) f0 Is f(0) included in the spectra? (T or F) xmin Minimum frequency for PLOTTING. xmax Maximum frequency for PLOTTING. ymin Minimum depth/height/time for PLOTTING. Only used if processing time-frequency results. ymax Maximum depth/height/time for PLOTTING. Only used if processing time-frequency results. npts The number of points in the processed time series window. This is needed for proper spectrum normalization. pad The total padded length of the processed time series window. This is needed for proper spectrum normalization. ydir Direction for y-axis in plots (depth or height). -1 = values increase downwards (slower plotting!), 1 = values increase upwards. Only used if processing time-frequency results. palette What color palette would you like to use? (1) rainbow, (2) grayscale, (3) blue, (4) red, (5) blue-white-red, (6) viridis ncolors Number of colors to use in plot. Only used if processing time-frequency results. h Height of plot in inches. w Width of plot in inches. ln Plot natural log of spectral results? (T or F) genplot Generate summary plots? (T or F) verbose Verbose output? (T or F)

### Details

Depending on the normalization used, you may want to preprocess the power spectra prior to integration.

eha

### Examples

# generate etp signal over past 10 Ma
ex=etp(tmax=10000)

# evolutive power