MODEL:
! A 70 gates 995 Flights Assignment Problem;
SETS:
   GATES: CAPACITY;
   FLIGHTS: DEMAND;
LINKS( GATES, FLIGHTS): COST, VOLUME;
ENDSETS
! Here is the data;
DATA:
   !set members;
   GATES= GATE1
GATE2
GATE3
GATE4

Yi<-c(2.1,2.5,4.9,5.5,7,8.4,9.6,10.2,11.4,12.5,13.1,14.6,17,16.8,18.6,19.7,21.3,21.6)
Xi<-c(1,1.5,2,3,4,5,6,7.5,8.5,10,12.5,15,17.5,20,25,30,35,40)
Y<-1/Yi
X<-1/Xi
lm(Y~X)
nls(Yi~(a*Xi)/(b+Xi),start=list(a=29.62085,b=13.44816),trace=TRUE)
plot(Xi,Yi)
curve((28.14*x)/(12.57+x),add=TRUE,0,40)
summary( nls(Yi~(a*Xi)/(b+Xi),start=list(a=29.62085,b=13.44816),trace=TRUE))
res<-resid(nls(Yi~(a*Xi)/(b+Xi),start=list(a=29.62085,b=13.44816),trace=TRUE))
fit<-fitted(nls(Yi~(a*Xi)/(b+Xi),start=list(a=29.62085,b=13.44816),trace=TRUE))
plot(fit, res)
plot(Xi, res)
qqnorm(res)
confint(nls(Yi~(a*Xi)/(b+Xi),start=list(a=29.62085,b=13.44816),trace=TRUE))

Yi<-c(0,0,1,0,0,1,1,1,0,1,1,0,0,1,0,0,1,0,1,0,0,1,1,0,0,1,0,0,1,0,0,1,0)
Xi1<-c(32,45,60,53,25,68,82,38,67,92,72,21,26,40,33,45,61,16,18,22,27,35,40,10,24,15,23,19,22,61,21,32,17)
Xi2<-c(3,2,2,1,4,1,2,5,2,2,3,5,3,4,3,1,2,3,4,6,3,3,3,4,3,4,3,5,2,2,3,5,1)
mod1<-lm(Yi~Xi1+Xi2)
mod2<-glm(Yi~Xi1+Xi2,family=binomial)
mod2
sum((Yi-33*fitted(mod2))^2 / 33*fitted(mod2)*(1-fitted(mod2)))
qchisq(.99,30)
res<-resid(mod2)
plot(c(1:33),res)
anova(mod2, mod1,test="Chisq")

Xi<-c(1, 0, 2, 0, 3, 1, 0, 1, 2, 0)
Yi<-c(16, 9, 17, 12, 22, 13, 8, 15, 19, 11)
plot(Xi,Yi)
glm(Yi~Xi,poisson)
anova(glm(Yi~Xi,poisson))
summary(glm(Yi~Xi,poisson))
deviance(glm(Yi~Xi,poisson))
sum(((Yi-fitted(glm(Yi~Xi,poisson))) / fitted(glm(Yi~Xi,poisson))^.5 )^2)
res<-resid(glm(Yi~Xi,poisson))
res
plot(c(1:10),res)
plot(glm(Yi~Xi,poisson))
predict(glm(Yi~Xi,poisson), data.frame(Xi=0), type="response")
mod2<-lm(Yi~Xi)
predict(lm(Yi~Xi), data.frame(Xi=0), type="response")

x<-c(1, 0, 2, 0, 3, 1, 0, 1, 2, 0)
y<-c(16, 9, 17, 12, 22, 13, 8, 15, 19, 11)
plot(x,y)
abline(lm(y~x))
curve(exp(0.2638*x+2.3529),add=TRUE,col="red",0,3)
predict(glm(Yi~Xi,poisson), data.frame(Xi=0), type="response")

x<-c(0,0,0,0)
y<-c(9,12,8,11)
glm(y~x, poisson)
confint(glm(Yi~Xi,poisson))

int nbFlt = …;
int nbGate = …;

range Gate 1 .. nbGate;
range Flt 1 .. nbFlt;

var
   int assign[Flt,Gate] in 0..1,
   int y[Flt,Flt] in 0..1;

float+ arrtm[Flt] = …;
float+ dptm[Flt] = …;
//var float+ interval[Flt, Flt];
                     
minimize
     sum (i in Flt, j in Flt : arrtm[j] - dptm[i] >0 & (j-i<=70 \/ i-j<=70) ) y[i,j] / (arrtm[j] - dptm[i])
      //sum (i in Flt, j in Flt : arrtm[j] - dptm[i] >0 /\ (j-i<=70 \/ i-j<=70 ) y[i,j] *(-1)* (arrtm[j] - dptm[i])
 
 
subject to {
   
           //forall(i in Flt, j in Flt : ij )
          //– forall(i in Flt, j in Flt : j-i<=70 \/ i-j<=70 )
            forall(i in Flt, j in Flt : arrtm[j] - dptm[i] >0 => j-i<=70 \/ i-j<=70  )
             {
              //i j;
              sum(k in Gate)assign[i,k] * assign[j,k] = y[i,j];        
               };
               
            forall(i in Flt)
            sum(k in  Gate) assign[i,k] = 1;
                    
           
           //forall(i in Flt, j in Flt)
          //   i j;
                     
           };