dS2_nlin2 <- function(theta, y)
{
y1 <- y[-length(y)]
y2 <- y[-1]
res <- (y2 - y1 - theta[3] - theta[1] * y1^theta[2])
r1 <- res[seq(1, floor((length(res) - 1) / 2), 2)]
r2 <- res[seq(2, floor((length(res) - 1) / 2), 2)]
r3 <- res[seq((length(res) - floor((length(res) - 1) / 2) + 1), length(res), 2)]
r4 <- res[seq((length(res) - floor((length(res) - 1) / 2)), length(res), 2)]
m <- min(c(length(r1), length(r2), length(r3)))
r1 <- r1[1:m]
r2 <- r2[1:m]
r3 <- r3[m:1]
r4 <- r4[m:1]
InD <- (r1 > 0) * (r2 < 0) * (r3 > 0) * (r4 < 0) + (r1 < 0) * (r2 > 0) * (r3 < 0) * (r4 > 0) #+(1-(r1!=0)*(r2!=0)*(r3!=0)*(r4!=0))
depth<-1 / (m) * sum(InD)
return(depth)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.