func (c *L1VM) PredictVector(x *core.Vector) float64 {
s := core.NewSample()
for k, xs := range c.sv {
s.AddFeature(core.Feature{Id: int64(k), Value: RBFKernel(xs, x, c.radius)})
}
return c.ftrl.Predict(s)
}
func (c *L1VM) Train(dataset *core.DataSet) {
c.sv = []*core.Vector{}
kernel_dataset := core.NewDataSet()
positive := []int{}
negative := []int{}
for i, si := range dataset.Samples {
if si.Label > 0.0 {
positive = append(positive, i)
} else {
negative = append(negative, i)
}
}
perm_positive := rand.Perm(len(positive))
for i, k := range perm_positive {
if i > c.count {
break
}
c.sv = append(c.sv, dataset.Samples[positive[k]].GetFeatureVector())
}
perm_negative := rand.Perm(len(negative))
for i, k := range perm_negative {
if i > c.count {
break
}
c.sv = append(c.sv, dataset.Samples[negative[k]].GetFeatureVector())
}
for _, si := range dataset.Samples {
xi := si.GetFeatureVector()
tsample := core.NewSample()
tsample.Label = si.Label
for j, xj := range c.sv {
tsample.AddFeature(core.Feature{Id: int64(j), Value: RBFKernel(xi, xj, c.radius)})
}
kernel_dataset.AddSample(tsample)
}
c.ftrl.Train(kernel_dataset)
}
func main() {
train_path, test_path, pred_path, _, params := hector.PrepareParams()
total := 5
methods := []string{"ftrl", "fm"}
all_methods_predictions := [][]*eval.LabelPrediction{}
all_methods_test_predictions := [][]*eval.LabelPrediction{}
for _, method := range methods {
fmt.Println(method)
average_auc := 0.0
all_predictions := []*eval.LabelPrediction{}
for part := 0; part < total; part++ {
train, test, _ := SplitFile(train_path, total, part)
classifier := hector.GetClassifier(method)
auc, predictions, _ := hector.AlgorithmRun(classifier, train, test, "", params)
fmt.Println("AUC:")
fmt.Println(auc)
average_auc += auc
os.Remove(train)
os.Remove(test)
classifier = nil
for _, pred := range predictions {
all_predictions = append(all_predictions, pred)
}
}
all_methods_predictions = append(all_methods_predictions, all_predictions)
fmt.Println(average_auc / float64(total))
classifier := hector.GetClassifier(method)
fmt.Println(test_path)
_, test_predictions, _ := hector.AlgorithmRun(classifier, train_path, test_path, "", params)
all_methods_test_predictions = append(all_methods_test_predictions, test_predictions)
}
var wait sync.WaitGroup
wait.Add(2)
dataset := core.NewDataSet()
go func() {
for i, _ := range all_methods_predictions[0] {
sample := core.NewSample()
sample.Label = all_methods_predictions[0][i].Label
for j, _ := range all_methods_predictions {
feature := core.Feature{Id: int64(j), Value: all_methods_predictions[j][i].Prediction}
sample.AddFeature(feature)
}
dataset.Samples <- sample
}
close(dataset.Samples)
wait.Done()
}()
ensembler := lr.LinearRegression{}
go func() {
ensembler.Init(params)
ensembler.Train(dataset)
wait.Done()
}()
wait.Wait()
fmt.Println(ensembler.Model)
wait.Add(2)
test_dataset := hector.NewDataSet()
go func() {
for i, _ := range all_methods_test_predictions[0] {
sample := hector.NewSample()
sample.Label = all_methods_test_predictions[0][i].Prediction
for j, _ := range all_methods_test_predictions {
feature := hector.Feature{Id: int64(j), Value: all_methods_test_predictions[j][i].Prediction}
sample.AddFeature(feature)
}
test_dataset.Samples <- sample
}
close(test_dataset.Samples)
wait.Done()
}()
go func() {
pred_file, _ := os.Create(test_path + ".out")
for sample := range test_dataset.Samples {
prediction := sample.Label //ensembler.Predict(sample)
pred_file.WriteString(strconv.FormatFloat(prediction, 'g', 5, 64) + "\n")
}
defer pred_file.Close()
wait.Done()
}()
wait.Wait()
}