integratePower {astrochron} | R Documentation |

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.

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)

`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) |

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

# 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)

[Package *astrochron* version 1.1 Index]