R Graphical Manual

Browse All

Last data update: 2014.03.03

R: plot a matrix of several seismograms

PLOT.MATN

R Documentation

plot a matrix of several seismograms

Description

Matrix of several seismograms

Usage

PLOT.MATN(ascd, tim=1, dt=1,  WIN=c(0,1), labs="",
notes=notes, sfact=1,ampboost=0,  shift=NULL, LOG="",
COL='red', add=1, AXES=1, units=NULL, VS=FALSE)

Arguments

`ascd`	N by K matrix of seismograms where
`tim`	time values fo x-axis
`dt`	sample interval, seconds
`WIN`	vector, time window for zoom
`labs`	vector of labels for each panel
`notes`	vector of notes for each panel
`sfact`	scaling factor, 1=window, 2=trace
`ampboost`	increase each amplitude by this multiplier
`shift`	vector, shift each trace by these time
`LOG`	log x-axis
`COL`	vector of colors or indexes to colors
`add`	numeric, to existing plot. add = 1,2,3 if add=1 plot and add traces, add =2 plot, but no traces, add = 3 no plot, but add traces. DEFAULT=1
`AXES`	numeric, 0,1,2,3,4; default=1
`units`	label for units of Y-axis
`VS`	var-squiggle display

Details

Plots a matrix of seismograms that each have the same starting time. For the AXES argument, 0 = no axes, AXES=1 plot scale for largest amplitude band and a multiplier for all others, AXES=2 left side, AXES=3 right side, AXES=4 alternate sides

Value

Graphical side effects and,

`n`	n
`windiv`	matrix of n rows, with columns=(window Y min, window Y max, user Y min, user Y max)

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

Examples


dt  <-  0.001

 t  <-  seq(0, 6, by=0.001)

thefreqs  <-  seq(from=10, to=100, by=10)
theamps  <-   runif(length(thefreqs))

#   sample rate is 1000 Hz, 0.001 seconds 601 samples
 x  <-  NULL

for(i in 1:length(thefreqs))
{
x  <-  cbind(x, theamps[i]*sin(2*pi*thefreqs[i]*t))


}

PLOT.MATN(x,  dt = dt)

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(RSEIS)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/RSEIS/PLOT.MATN.Rd_%03d_medium.png", width=480, height=480)
> ### Name: PLOT.MATN
> ### Title: plot a matrix of several seismograms
> ### Aliases: PLOT.MATN
> ### Keywords: hplot
> 
> ### ** Examples
> 
> 
> dt  <-  0.001
> 
>  t  <-  seq(0, 6, by=0.001)
> 
> thefreqs  <-  seq(from=10, to=100, by=10)
> theamps  <-   runif(length(thefreqs))
> 
> #   sample rate is 1000 Hz, 0.001 seconds 601 samples
>  x  <-  NULL
> 
> for(i in 1:length(thefreqs))
+ {
+ x  <-  cbind(x, theamps[i]*sin(2*pi*thefreqs[i]*t))
+ 
+ 
+ }
> 
> PLOT.MATN(x,  dt = dt)
> 
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>