M<-1794
m<-138
T<-13
case<-matrix(0,nrow=138,ncol=13)
pop<-matrix(0,nrow=138,ncol=13)
for( i in 1: m){

case[i,1]<-CASE1[i]
case[i,2]<-CASE2[i]
case[i,3]<-CASE3[i]
case[i,4]<-CASE4[i]
case[i,5]<-CASE5[i]
case[i,6]<-CASE6[i]
case[i,7]<-CASE7[i]
case[i,8]<-CASE8[i]
case[i,9]<-CASE9[i]
case[i,10]<-CASE10[i]
case[i,11]<-CASE11[i]
case[i,12]<-CASE12[i]
case[i,13]<-CASE13[i]
pop[i,1]<-POP1[i]
pop[i,2]<-POP2[i];
pop[i,3]<-POP3[i]
pop[i,4]<-POP4[i]
pop[i,5]<-POP5[i]
pop[i,6]<-POP6[i]
pop[i,7]<-POP7[i]
pop[i,8]<-POP8[i]
pop[i,9]<-POP9[i]
pop[i,10]<-POP10[i]
pop[i,11]<-POP11[i]
pop[i,12]<-POP12[i]
pop[i,13]<-POP13[i]}

yL<-rep(0,M)
pL<-rep(0,M)
T=13
for (i in 1:138){
for (j in 1:13){
k<-j+T*(i-1)
yL[k]<-case[i,j]
pL[k]<-pop[i,j]
}}
timeP<-rep(1:13,length=M)
region<-rep(1:138,each=13)
region2<-region
ind2<-rep(1:M)
data<-data.frame(yL,pL,T,m,M,timeP,region,region2,ind2)


data<-data.frame(yL,pL,T,m,M)