IDoubleArray hmmChi = pmm.getOutputParameters();
io.writeString(cmd.outdir+"/hmmChi.dat", doubles.toString(hmmChi,"\t","\n"));
}
else if (cmd.hmmTimescalesEstimate)
{
IIntList lagtimes = intsNew.list(0);
lagtimes.append(cmd.taumin);
for (double tau = cmd.taumin; tau <= cmd.taumax; tau *= cmd.taumult)
{
int lag = (int)tau;
if (lag != lagtimes.get(lagtimes.size()-1))
lagtimes.append(lag);
}
PrintStream itsout = new PrintStream(cmd.outdir+"/hmm-its.dat");
IDoubleArray lastTC = cmd.init_T;
IDoubleArray lastChi = cmd.init_Chi;
for (int i=0; i<lagtimes.size(); i++)
{
//set lag and timeshift
int lag = lagtimes.get(i);
System.out.println("\ntau = "+lag+"\n");
// in direct mode, erase the results from last lag. Otherwise use as initialization.
if (cmd.direct)
{