Nothing
matsquiggle<-function(XMAT, dt1, dist=NULL, thick=1 , FLIP=FALSE, filcol='blue',
tracecol="black", add=FALSE, PLOT=TRUE, xpd=TRUE, plotdir=1)
{
##### plot a seismic section as in exploration seismology
#### with varian - squiggle display
## source("varsquig.R")
d = dim(XMAT)
x = seq(from=0, length=d[1], by=dt1)
if(missing(dist))
{
y1 = seq(from=0, by=1, length=d[2])
}
else
{
y1 = dist
}
### message(y1)
Ry = range(y1)
### fatness = (y1[2]-y1[1])/2
if(missing(thick))
{
fatness = diff(Ry)/length(y1)
}
else
{
fatness = thick
}
if(is.null(fatness))
{
fatness = diff(Ry)/length(y1)
}
if(missing(add)) { add= FALSE }
if(missing(PLOT)) { PLOT=TRUE }
### if add = FALSE do not start a new plot
if(!add){
if(FLIP)
{
plot(range(y1), range(x), type='n', axes=FALSE, xlab="", ylab="" )
}
else
{
plot(range(x),range(y1) , type='n', axes=FALSE, xlab="", ylab="" )
}
}
if(missing(filcol)) { filcol=rep('blue', length(y1)) }
if(missing(tracecol)) { tracecol=rep('blue', length(y1)) }
if(length(filcol)==1) { filcol=rep(filcol , length(y1)) }
if(length(tracecol)==1) { tracecol=rep(tracecol , length(y1)) }
if(PLOT)
{
if(plotdir==1) { pord = 1:length(y1) }
if(plotdir==0) { pord = length(y1):1 }
for(i in pord)
{
if(FLIP)
{
L = list(y=range(x) , x =c(y1[i]-fatness, y1[i]+fatness ) )
}
else
{
L = list(x=range(x) , y =c(y1[i]-fatness, y1[i]+fatness ) )
}
varsquig(x, XMAT[,i], L=L , FLIP=FLIP, filcol=filcol[i], tracecol=tracecol[i], var=TRUE, xpd=xpd)
}
}
invisible(list(thick=fatness, rx = range(x), ry=range(y1) ) )
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.