#!/bin/bash
# Georges Kunstler 14/05/2014
# read one variable
export LD_LIBRARY_PATH=/usr/lib64/R/library

mkdir -p Rscript_temp


# test parameter
nbargs=$#
echo "number of arguments="$nbargs
if [ $nbargs -ne 1 ]
then
  echo "need one and only one argument"
  echo " usage :"
  echo " ./launch_all_lmer.sh  sample.size"
  exit 100
fi


samplesize=$1
# "'Seed.mass'" "'Leaf.N'"
 for trait in "'SLA'" "'Wood.density'" "'Max.height'" ; do


# ALL data
  	echo "/usr/local/R/R-3.1.1/bin/Rscript -e \"source('R/analysis/lmer.run.R'); run.multiple.model.for.set.one.trait(model.files.lmer.Tf.3[1], run.lmer,$trait);print('done')\"" > Rscript_temp/allf${trait}.sh
	qsub Rscript_temp/allf${trait}.sh -d ~/trait.competition.workshop -l nodes=1:ppn=1,mem=8gb -N "lmerall2all.f${trait}" -q opt32G -j oe

  	echo "/usr/local/R/R-3.1.1/bin/Rscript -e \"source('R/analysis/lmer.run.R'); run.multiple.model.for.set.one.trait(model.files.lmer.Tf.3[2], run.lmer,$trait);print('done')\"" > Rscript_temp/allf2${trait}.sh
	qsub Rscript_temp/allf2${trait}.sh -d ~/trait.competition.workshop -l nodes=1:ppn=1,mem=8gb -N "lmerall2all.fb${trait}" -q opt32G -j oe

  	echo "/usr/local/R/R-3.1.1/bin/Rscript -e \"source('R/analysis/lmer.run.R'); run.multiple.model.for.set.one.trait(model.files.lmer.Tf.3[3], merge.biomes.TF = TRUE, run.lmer,$trait);print('done')\"" > Rscript_temp/allf3${trait}.sh
	qsub Rscript_temp/allf3${trait}.sh -d ~/trait.competition.workshop -l nodes=1:ppn=1,mem=8gb -N "lmerall2all.fb2${trait}" -q opt32G -j oe

# ALL data All census
  	echo "/usr/local/R/R-3.1.1/bin/Rscript -e \"source('R/analysis/lmer.run.R'); run.multiple.model.for.set.one.trait(model.files.lmer.Tf.4[1], run.lmer,$trait, select.one.census = FALSE);print('done')\"" > Rscript_temp/allcensus${trait}.sh
	qsub Rscript_temp/allcensus${trait}.sh -d ~/trait.competition.workshop -l nodes=1:ppn=1,mem=8gb -N "lmerall2all.census${trait}" -q opt32G -j oe

  	echo "/usr/local/R/R-3.1.1/bin/Rscript -e \"source('R/analysis/lmer.run.R'); run.multiple.model.for.set.one.trait(model.files.lmer.Tf.4[2], run.lmer,$trait, select.one.census = FALSE);print('done')\"" > Rscript_temp/allcensus2${trait}.sh
	qsub Rscript_temp/allcensus2${trait}.sh -d ~/trait.competition.workshop -l nodes=1:ppn=1,mem=8gb -N "lmerall2all.censusb${trait}" -q opt32G -j oe

  	echo "/usr/local/R/R-3.1.1/bin/Rscript -e \"source('R/analysis/lmer.run.R'); run.multiple.model.for.set.one.trait(model.files.lmer.Tf.4[3], merge.biomes.TF = TRUE, run.lmer,$trait,select.one.census = FALSE);print('done')\"" > Rscript_temp/allcensus3${trait}.sh
	qsub Rscript_temp/allcensus3${trait}.sh -d ~/trait.competition.workshop -l nodes=1:ppn=1,mem=8gb -N "lmerall2all.censusb2${trait}" -q opt32G -j oe



done

  	# echo "/usr/local/R/R-3.1.1/bin/Rscript -e \"source('R/analysis/lmer.run.R'); run.multiple.model.for.set.one.trait(c(model.files.lmer.Tf.Multi[1]), run.lmer,'Multi', data.type ='Multi');print('done')\"" > Rscript_temp/Multi1.sh
	# qsub Rscript_temp/Multi1.sh -d ~/trait.competition.workshop -l nodes=1:ppn=1,mem=8gb -N "Multi1" -q opt32G -j oe

#  	echo "/usr/local/R/R-3.1.1/bin/Rscript -e \"source('R/analysis/lmer.run.R'); run.multiple.model.for.set.one.trait(c(model.files.lmer.Tf.Multi[2]), run.lmer,'Multi', data.type ='Multi');print('done')\"" > Rscript_temp/Multi2.sh
#	qsub Rscript_temp/Multi2.sh -d ~/trait.competition.workshop -l nodes=1:ppn=1,mem=8gb -N "Multi2" -q opt32G -j oe

  	# echo "/usr/local/R/R-3.1.1/bin/Rscript -e \"source('R/analysis/lmer.run.R'); run.multiple.model.for.set.one.trait(c(model.files.lmer.Tf.Multi[3]), run.lmer,'Multi', data.type ='Multi');print('done')\"" > Rscript_temp/Multi3.sh
	# qsub Rscript_temp/Multi3.sh -d ~/trait.competition.workshop -l nodes=1:ppn=1,mem=8gb -N "Multi3" -q opt32G -j oe