DoubleGauss2 function used in fitting nlme curve for observations
doubleGauss2(dat, y, time, params = NULL, concave = TRUE, ...)subject data to be used
outcome variable, character vector
time variable, character vector
NULL unless user wants to specify starting parameters for gnls
Boolean
just in case
User should only have to worry about setting concavity of this function. Presently only work for time series scaled out to 2000ms
y ~ (time < mu) * (exp(-1 * (time - mu) ^ 2
/ (2 * sig1 ^ 2)) * (ht - base1) + base1)
+ (mu <= time) * (exp(-1 * (time - mu) ^ 2
/ (2 * sig2 ^ 2)) * (ht - base2) + base2)