library(R2admb)  ## requires 0.7.2
sessionInfo()

K<-5
n<-7
y<-matrix(data = NA, nrow = K, ncol = n)
x <- c(118.00, 484.00, 664.00, 1004.00, 1231.00, 1372.00, 1582.00)
phi2<-728.8
phi3<-353.5
sigma<-22.3
mu1<-192.1
sigma1=31.6
phi1<-rnorm(K,mean=mu1,sd=sigma1)
for(i in 1:K) {
y[i,]<-phi1[i] / (1 + exp(-(x-phi2)/phi3))
}
y<-y+matrix(rnorm(n*K, mean = 0, sd = sigma),nrow = K, ncol = n)


out<-list(K=K,n=n,x=x,y=y)

filen<-'OrangeTree.dat'

cat('# init_int K\n',K,'\n', file=filen)
cat('# init_int n\n',n,'\n', file=filen,append=TRUE)
cat('# init_vector x(1,n)\n',x,'\n', file=filen,append=TRUE)
cat('# init_matrix Y(1,K,1,n)\n', file=filen,append=TRUE)
write.table(y,row.names=FALSE,col.names=FALSE, quote=FALSE, append=TRUE, file=filen)
cat('\n', file=filen,append=TRUE)
#write_dat("OrangeTree",out)