psdPlot {IRISSeismic} | R Documentation |
Generate plots from a set of PSDs
Description
The psdPlot
function is used to generate plots from the data in a psdList
or psdDF
dataframe.
Usage
psdPlot(PSDs, style='psd', evalresp=NULL, ylo=-200, yhi=-50, showNoiseModel=TRUE,
showMaxMin=TRUE, showMode=TRUE, showMean=FALSE, showMedian=FALSE, ...)
Arguments
PSDs |
either a list as returned by |
style |
character identifier of plot type: |
evalresp |
dataframe of freq, amp, phase information matching output of |
ylo |
numeric setting lower limit of plot y-axis (default= |
yhi |
numeric setting upper limit of plot y-axis (default= |
showNoiseModel |
logical controlling plotting of noise model lines (default= |
showMaxMin |
logical controlling plotting of PSD max and min lines (default= |
showMode |
logical controlling plotting of PDF mode line (default= |
showMean |
logical controlling plotting of PSD mean line (default= |
showMedian |
logical controlling plotting of PSD median line (default= |
... |
arguments to be passed to plotting methods |
Details
The psdPlot
function creates visualizations for sets of PSDs. Plots generated with style='pdf'
mimic the plots presented in the McNamara paper.
Author(s)
Jonathan Callahan jonathan@mazamascience.com
References
Seismic Noise Analysis System Using Power Spectral Density Probability Density Functions (McNamara and Boaz 2005)
See Also
McNamaraPSD
,
psdList
,
psdStatistics
Examples
## Not run:
# Create a new IrisClient
iris <- new("IrisClient", debug=TRUE)
# Get seismic data
starttime <- as.POSIXct("2011-05-05", tz="GMT") # 2011.125
endtime <- starttime + 1*24*3600
st <- getDataselect(iris,"IU","GRFO","--","BHE",starttime,endtime)
# Generate power spectral density for each hour long segment
psdList <- psdList(st)
# 'psd' line plot
psdPlot(psdList,style='psd',type='l',col=adjustcolor('black',0.3))
# McNamara 'pdf' plot
psdPlot(psdList,style='pdf')
## End(Not run)