## 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)
Run the code above in your browser using DataLab