# integratePower: Determine the total power within a given bandwidth In astrochron: A Computational Tool for Astrochronology

 integratePower 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,
xmax=NULL,ymin=NULL,ymax=NULL,npts=NULL,pad=NULL,ydir=1,
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
pwr=eha(ex,win=500,fmax=.1,pad=2000,output=2,pl=2)

# integrate power from main obliquity term
integratePower(pwr,flow=0.02,fhigh=0.029,npts=501,pad=2000)
``````

astrochron documentation built on Aug. 26, 2023, 5:07 p.m.