minruns <- 1
maxruns <- 100

output_index <- c('A','B','C','D','E','F')

write("!#/usr/local/bin/bash",file="BeagleRuns.sh",append=FALSE,ncolumns=100)
write("!#/usr/local/bin/bash",file="PlinkRuns.sh",append=FALSE,ncolumns=100)

for (j in 1:length(output_index)) {

for (i in minruns:maxruns) {

 out1 <- paste("./Brun",output_index[j],i,".sh > OutBrun",output_index[j],i,".txt 2>&1 &",sep="")
write(out1,file="BeagleRuns.sh",append=TRUE,ncolumns=100)

 out2 <- paste("./Prun",output_index[j],i,".sh > OutPrun",output_index[j],i,".txt 2>&1 &",sep="")
write(out2,file="PlinkRuns.sh",append=TRUE,ncolumns=100)

 
 com <- paste(
"echo 'minruns <- ",i,"' > ddu
echo 'maxruns <- ",i,"' >> ddu
echo 'output_index <- \"",output_index[j],i,"\"' >> ddu",sep="")

system(com)

if (j<5) {
com1 <- paste("cat be1 ddu be2 > compareMethodsBeagle",output_index[j],i,".R",sep="")
} else {
com1 <- paste("cat be1 ddu be2A > compareMethodsBeagle",output_index[j],i,".R",sep="")
}

system(com1)

if (j<5) {
com1a <- paste("cat pe1 ddu pe2 > compareMethodsPlink",output_index[j],i,".R",sep="")
} else {
com1a <- paste("cat pe1 ddu pe2A > compareMethodsPlink",output_index[j],i,".R",sep="")
}

system(com1a)


com2 <- paste(
"echo '!#/usr/local/bin/bash' > Brun",output_index[j],i,".sh

echo 'cd /seppdata/sepp/linkage/SFSsimulation/sim",output_index[j],"long' >> Brun",output_index[j],i,".sh

echo 'cp /seppdata/sepp/linkage/SFSsimulation/Lscripts/compareMethodsBeagle",output_index[j],i,".R .'  >> Brun",output_index[j],i,".sh

echo '/system/apps/biosoft/R-2.14.0/bin/R -f compareMethodsBeagle",output_index[j],i,".R' >> Brun",output_index[j],i,".sh

echo 'cd /seppdata/sepp/linkage/SFSsimulation' >> Brun",output_index[j],i,".sh",sep="")

system(com2)

com2a <- paste(
"echo '!#/usr/local/bin/bash' > Prun",output_index[j],i,".sh

echo 'cd /seppdata/sepp/linkage/SFSsimulation/sim",output_index[j],"long' >> Prun",output_index[j],i,".sh

echo 'cp /seppdata/sepp/linkage/SFSsimulation/Lscripts/compareMethodsPlink",output_index[j],i,".R .'  >> Prun",output_index[j],i,".sh

echo '/system/apps/biosoft/R-2.14.0/bin/R -f compareMethodsPlink",output_index[j],i,".R' >> Prun",output_index[j],i,".sh

echo 'cd /seppdata/sepp/linkage/SFSsimulation' >> Prun",output_index[j],i,".sh",sep="")

system(com2a)

}

}



minruns <- 1
maxruns <- 100

output_index <- c('A','B','C','D','E','F')

write("!#/usr/local/bin/bash",file="PlinkRunsNew.sh",append=FALSE,ncolumns=100)

blocks <- 8
classR <- (minruns:maxruns)%/%8
noClassR <- classR[maxruns]+1

elemR <- rep(0,noClassR)
for (c in 1:noClassR) {
a <- which(classR==(c-1))
elemR[c] <- length(a)
}

act <- 1


for (c in 1:noClassR) {

minrunsT <- act
act <- act+elemR[c]
maxrunsT <- act-1

for (j in 1:length(output_index)) {

for (i in minrunsT:maxrunsT) {


 out2 <- paste("./Prun",output_index[j],i,".sh > OutPrun",output_index[j],i,".txt 2>&1 &",sep="")
write(out2,file="PlinkRunsNew.sh",append=TRUE,ncolumns=100)

}
}
}




#########COBRA#####################
#########COBRA#####################
#########COBRA#####################
#########COBRA#####################





minruns <- 1
maxruns <- 100

output_index <- c('A','B','C','D','E','F')

write("!#/usr/local/bin/bash",file="BeagleRunsCobra.sh",append=FALSE,ncolumns=100)
write("!#/usr/local/bin/bash",file="PlinkRunsCobra.sh",append=FALSE,ncolumns=100)

for (j in 1:length(output_index)) {

for (i in minruns:maxruns) {

 out1 <- paste("./Brun",output_index[j],i,".sh > OutBrun",output_index[j],i,".txt 2>&1 &",sep="")
write(out1,file="BeagleRunsCobra.sh",append=TRUE,ncolumns=100)

 out2 <- paste("./Prun",output_index[j],i,".sh > OutPrun",output_index[j],i,".txt 2>&1 &",sep="")
write(out2,file="PlinkRunsCobra.sh",append=TRUE,ncolumns=100)

 
 com <- paste(
"echo 'minruns <- ",i,"' > ddu
echo 'maxruns <- ",i,"' >> ddu
echo 'output_index <- \"",output_index[j],i,"\"' >> ddu",sep="")

system(com)

if (j<5) {
com1 <- paste("cat be1C ddu be2C > compareMethodsBeagleCobra",output_index[j],i,".R",sep="")
} else {
com1 <- paste("cat be1C ddu be2AC > compareMethodsBeagleCobra",output_index[j],i,".R",sep="")
}

system(com1)

if (j<5) {
com1a <- paste("cat pe1C ddu pe2C > compareMethodsPlinkCobra",output_index[j],i,".R",sep="")
} else {
com1a <- paste("cat pe1C ddu pe2AC > compareMethodsPlinkCobra",output_index[j],i,".R",sep="")
}

system(com1a)


com2 <- paste(
"echo '!#/usr/local/bin/bash' > Brun",output_index[j],i,".sh

echo 'cd /system/user/seppdata/sepp/linkage/SFSsimulation/sim",output_index[j],"long' >> Brun",output_index[j],i,".sh

echo 'cp /system/user/seppdata/sepp/linkage/SFSsimulation/Lscripts/compareMethodsBeagleCobra",output_index[j],i,".R .'  >> Brun",output_index[j],i,".sh

echo '/system/apps/biosoft/R-2.14.0/bin/R -f compareMethodsBeagleCobra",output_index[j],i,".R' >> Brun",output_index[j],i,".sh

echo 'cd /system/user/seppdata/sepp/linkage/SFSsimulation' >> Brun",output_index[j],i,".sh",sep="")

system(com2)

com2a <- paste(
"echo '!#/usr/local/bin/bash' > Prun",output_index[j],i,".sh

echo 'cd /system/user/seppdata/sepp/linkage/SFSsimulation/sim",output_index[j],"long' >> Prun",output_index[j],i,".sh

echo 'cp /system/user/seppdata/sepp/linkage/SFSsimulation/Lscripts/compareMethodsPlinkCobra",output_index[j],i,".R .'  >> Prun",output_index[j],i,".sh

echo '/system/apps/biosoft/R-2.14.0/bin/R -f compareMethodsPlinkCobra",output_index[j],i,".R' >> Prun",output_index[j],i,".sh

echo 'cd /system/user/seppdata/sepp/linkage/SFSsimulation' >> Prun",output_index[j],i,".sh",sep="")

system(com2a)

}

}



minruns <- 1
maxruns <- 100

output_index <- c('A','B','C','D','E','F')

write("!#/usr/local/bin/bash",file="PlinkRunsNewCobra.sh",append=FALSE,ncolumns=100)

blocks <- 8
classR <- (minruns:maxruns)%/%8
noClassR <- classR[maxruns]+1

elemR <- rep(0,noClassR)
for (c in 1:noClassR) {
a <- which(classR==(c-1))
elemR[c] <- length(a)
}

act <- 1


for (c in 1:noClassR) {

minrunsT <- act
act <- act+elemR[c]
maxrunsT <- act-1

for (j in 1:length(output_index)) {

for (i in minrunsT:maxrunsT) {


 out2 <- paste("./Prun",output_index[j],i,".sh > OutPrun",output_index[j],i,".txt 2>&1 &",sep="")
write(out2,file="PlinkRunsNewCobra.sh",append=TRUE,ncolumns=100)

}
}
}











/system/user

echo $TMPDIR

export TMPDIR=/BIGtmp/scratch/k3307/k337870/tmp