From 8dc8da95179d0f0654b5cb9e189117197aec7c55 Mon Sep 17 00:00:00 2001
From: Benjamin Auder <benjamin@auder>
Date: Mon, 28 Oct 2019 22:45:24 +0100
Subject: [PATCH] Save multistart results in different files

---
 reports/multistart.R    |  2 +-
 reports/run_accu_cl.sh  | 18 ++++++++----------
 reports/run_multi_cl.sh | 24 ++++++++++++++++++++++++
 3 files changed, 33 insertions(+), 11 deletions(-)
 create mode 100644 reports/run_multi_cl.sh

diff --git a/reports/multistart.R b/reports/multistart.R
index 5fa80ad..c7db8f2 100644
--- a/reports/multistart.R
+++ b/reports/multistart.R
@@ -107,4 +107,4 @@ ms <- testMultistart(N, n, K, p, beta, b, link, nstart, ncores)
 ms_params <- list("N"=N, "nc"=ncores, "n"=n, "K"=K, "d"=d, "link"=link,
 	"p"=c(p,1-sum(p)), "beta"=beta, "b"=b, "nstart"=nstart)
 
-save("ms", "ms_params", file="multistart.RData")
+save("ms", "ms_params", file=paste("res_",n,"_",d,"_",link,"_",nstart,".RData",sep=""))
diff --git a/reports/run_accu_cl.sh b/reports/run_accu_cl.sh
index 0cd877b..04f1519 100644
--- a/reports/run_accu_cl.sh
+++ b/reports/run_accu_cl.sh
@@ -1,26 +1,24 @@
 #!/bin/bash
 
+# Lancement: qsub -o .output -j y run_accu_cl.sh
+
 #$ -N morpheus
+#$ -wd /workdir2/auder/morpheus/reports
 #$ -m abes
 #$ -M benjamin@auder.net
-#$ -pe make 50
+#$ -pe make 10
 #$ -l h_vmem=1G
-#$ -j y
-#$ -o .output
 rm -f .output
 
-WORKDIR=/workdir2/auder/morpheus/reports
-cd $WORKDIR
-
 module load R/3.6.1
 
-N=10
+N=100
 n=1e5
-nc=50
+nc=10
 
-for d in 2 5 10 20; do
+for d in 2 5 10; do
 	for link in "logit" "probit"; do
-		for weights in "1,1,1" "6,3,1"; do
+		for weights in "1,1,0"; do
 			R --slave --args N=$N n=$n nc=$nc d=$d link=$link weights=$weights <accuracy.R >out_${n}_${link}_${d}_${weights} 2>&1
 		done
 	done
diff --git a/reports/run_multi_cl.sh b/reports/run_multi_cl.sh
new file mode 100644
index 0000000..cd5e493
--- /dev/null
+++ b/reports/run_multi_cl.sh
@@ -0,0 +1,24 @@
+#!/bin/bash
+
+# Lancement: qsub -o .output -j y run_multi_cl.sh
+
+#$ -N morpheus
+#$ -wd /workdir2/auder/morpheus/reports
+#$ -m abes
+#$ -M benjamin@auder.net
+#$ -pe make 50
+#$ -l h_vmem=1G
+rm -f .output
+
+module load R/3.6.1
+
+N=100
+n=1e5
+nc=50
+nstart=5
+
+for d in 2 5 10; do
+	for link in "logit" "probit"; do
+		R --slave --args N=$N n=$n nc=$nc d=$d link=$link nstart=$nstart <multistart.R >out_${n}_${link}_${d}_${nstart} 2>&1
+	done
+done
-- 
2.44.0