## Loading required package: boot
## Loading required package: DEoptim
## 
## DEoptim package
## Differential Evolution algorithm in R
## Authors: D. Ardia, K. Mullen, B. Peterson and J. Ulrich
## 
## Loading required package: dplyr
## 
## Attaching package: 'dplyr'
## 
## The following object is masked from 'package:stats':
## 
##     filter
## 
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
## 
## Loading required package: ggplot2
## Loading required package: tidyr

A few examples of the threshlines extending down and working in Alex’s proposed repository change.

library(MPDiR) # contains the Vernier data; use ?Venier for the reference

fit <- quickpsy(Vernier, Phaseshift, NumUpward, N, bootstrap='none',
                grouping = .(Direction, WaveForm, TempFreq))
plotcurves(fit)

gompertz_fun <- function(x, p) exp(-p[1] * exp(-p[2] * x))
fit <- quickpsy(qpdat, phase, resp, fun = gompertz_fun, parini = c(1,.01), bootstrap='none',
                guess=.3, grouping = .(interval))
plot1<-plotcurves(fit)
show(plot1)

A case where the bottom axis limit isn’t 0

compressed<- qpdat
compressed$resp<- qpdat$resp/2 + 0.5
fit <- quickpsy(compressed, phase, resp, bootstrap='none',
                guess=.3, grouping = .(interval))
plot1<-plotcurves(fit)
show(plot1)