Nothing
`GETARAIC` <-
function(z4, DT=0.008, Mar=8, O1=2, O2=0.2, WW= 2, T1=1 , PLOT=FALSE )
{
Nz4 = length(z4)
### the 10 here is 10 samples before and after the window
if(missing(DT)) { DT=0.008 }
if(missing(Mar)) { Mar = 8 }
if(missing(O1)) { O1=10*DT }
if(missing(O2)) { O2=10*DT }
if(missing(WW)) { WW=10*DT }
if(missing(PLOT)) { PLOT=FALSE }
oldpar <- par(no.readonly = TRUE)
on.exit(par(oldpar))
### WW = window width (s)
### O1 = time before (s)
### O2 = time after (s)
### T1 = sample of guess
if(missing(T1)) { T1 = O1/DT }
aout = rep(0, Nz4)
### k1=T1+floor(O2/DT)+1;
### k2 = k1+IW;
### message(c(k1, k2, IW, Nz4));
ary = .C("CALL_ARAIC", PACKAGE = "RSEIS",
as.double(z4), as.integer(Nz4),as.double(DT), as.integer(Mar),
as.integer(T1), as.double(O1), as.double(O2), as.double(WW), as.double(aout))
kaic = ary[[9]]
kaic[kaic==0]=NA
Taic =TFIN= which.min(kaic)
if(PLOT==TRUE)
{
par(mfrow=c(2,1))
plot.ts(z4)
if(any(!is.na(kaic)))
{
xkaic = 1:length(kaic)
plot(xkaic,kaic, type='l')
lm1 = lm(kaic ~ cbind(xkaic, xkaic^2, xkaic^3, xkaic^4))
lines(xkaic[!is.na(kaic)], lm1$fitted.values, col=2)
vline(Taic, COL=rgb(.4,.8,1) )
}
else
{
warning("ERROR (GETARAIC) - had some NA in kaic.")
}
}
return(Taic)
}
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.