Ticket #1400: ClusterUtils.tol

NameBlock Rhclust =
[[
  Text _.cwd = GetAbsolutePath(".");
  Text _.ClusterUtils.r = _.cwd+"/"+ "ClusterUtils.r";
  Text _.tmpdir = Replace(TmpDir, "\\", "/");
  Text _.dist.bbm = _.tmpdir+"/"+"dist%S.bbm";
  Text _.result.bbm = _.tmpdir+"/"+"result%S.bbm";
  
  Matrix Run( Matrix dist, Text method, Real size )
  {
    Real verbose = False;
    
    Set args.info =
      [[
        [[ method, "method" ]],
        [[ size, "size", "int" ]],
        [[ _.result.bbm, "result.bbm"]]
      ]];
    Set args.bbm =
    [[
      [[ _.dist.bbm, dist ]]
    ]];

    Text argList = Rapi::CreateArgList("Hclust_args", args.info);
    Set result.names =
    [[
      _.result.bbm
    ]];
    Text script =
      "source('"+_.ClusterUtils.r+"'); "
      "Hclust.TOL('"+_.dist.bbm+"','%A')";
    //WriteLn( "argList = " << argList );
    //WriteLn( "script = " << script );
    Set rrs = Rapi::RunRScript( script,
                                args.bbm, argList, result.names, verbose );
    Matrix If(Card(rrs), 
    {
       Matrix groups = rrs[1]
    },
    {
      WriteLn("[Hclust::Run] Cannot run R script");
      UnknownMatrix
    })
  }
]];