i=0
maxit=${1-10}
while [ $i -lt $maxit ]
do
  let i=$i+1
  echo "Iteration $i"
  twinsim < make300.in > /dev/null
  sib-pair < fam300.in 
  merlin2R
  R --slave < varcom.R
#  R --slave < shap.R 
done
