I\'m just trying to run this simple model
#simulate some data set.seed(2000) x<-seq(0,50,1) y<-((runif(1,10,20)*x)/(runif(1,0,10)+x))+rnorm(51,0,1) #for sim