I am running HOF models in R eHOF package.
library(eHOF) set.seed(100) #create some data occ<-rbinom(100,1,0.5) grad<-seq(1,50,len=100) #run the model H