#-------------------------------------------------------------
#
# Licensed to the Apache Software Foundation (ASF) under one
# or more contributor license agreements.  See the NOTICE file
# distributed with this work for additional information
# regarding copyright ownership.  The ASF licenses this file
# to you under the Apache License, Version 2.0 (the
# "License"); you may not use this file except in compliance
# with the License.  You may obtain a copy of the License at
#
#   http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing,
# software distributed under the License is distributed on an
# "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
# KIND, either express or implied.  See the License for the
# specific language governing permissions and limitations
# under the License.
#
#-------------------------------------------------------------

# data preparation
FXY = read("./src/test/resources/datasets/Salaries.csv",
  data_type="frame", format="csv", header=TRUE);
F = FXY[,1:ncol(FXY)-1];
y = as.matrix(FXY[,ncol(FXY)]);
jspec= "{ ids:true, recode:[1,2,3,6],bin:[" 
      +"{id:4, method:equi-width, numbins:14},"
      +"{id:5, method:equi-width, numbins:12}]}"
[X,M] = transformencode(target=F, spec=jspec);
X = X[,2:ncol(X)]

m = nrow(X)
n = ncol(X)
fdom = colMaxs(X);
foffb = t(cumsum(t(fdom))) - fdom;
foffe = t(cumsum(t(fdom)))
rix = matrix(seq(1,m)%*%matrix(1,1,n), m*n, 1)
cix = matrix(X + foffb, m*n, 1);
X2 = table(rix, cix); #one-hot encoded

# learn model
B = lm(X=X2, y=y, verbose=FALSE);
yhat = X2 %*% B;
e = (y-yhat)^2;

# call slice finding
[TS,TR,d] = slicefinder(X=X, e=e, k=10,
  alpha=0.95, minSup=4, tpEval=TRUE, verbose=TRUE);

print("TS:\n" + toString(TS))
print("TR:\n" + toString(TR))
print("Debug matrix:\n" + toString(d))
