R/derived.2.r

In hysteresis: Tools for Modeling Rate-Dependent Hysteretic Processes and Ellipses

derived.2 <- function(b.x,b.y,retention,period) {
area<-abs(pi*b.x*retention)
coercion <- b.x/(sqrt(1+(b.y/retention)^2))
lag<-abs(atan2(retention,b.y))*period/(pi*2)
lag <- ifelse(lag > period,NA,lag)
z <- c(area,lag,coercion)
z
}