allow seeding of V and added documentation

1 files changed, 100 insertions, 0 deletions

diff --git a/src/predMSVMMaj.c b/src/predMSVMMaj.c
new file mode 100644
index 0000000..5d26dc3
--- /dev/null
+++ b/src/predMSVMMaj.c
@@ -0,0 +1,100 @@
+#include "libMSVMMaj.h"
+
+#define MINARGS 3
+
+void print_null(const char *s) {}
+void exit_with_help();
+void parse_command_line(int argc, char **argv, struct Model *model,
+		char *input_filename, char *output_filename, 
+		char *model_filename);
+
+void exit_with_help()
+{
+	printf("This is MSVMMaj, version %1.1f\n\n", VERSION);
+	printf("Usage: predMSVMMaj [options] test_data_file model_file\n");
+	printf("Options:\n");
+	printf("-o output_file : write output to file\n");
+	printf("-q : quiet mode (no output)\n");
+	exit(0);
+}
+
+int main(int argc, char **argv)
+{
+	long *predy;
+	double performance;
+
+	char input_filename[MAX_LINE_LENGTH];
+	char model_filename[MAX_LINE_LENGTH];
+	char output_filename[MAX_LINE_LENGTH];;
+
+	struct Model *model = Malloc(struct Model, 1);
+	struct Data *data = Malloc(struct Data, 1);
+
+	if (argc < MINARGS || check_argv(argc, argv, "-help") || check_argv_eq(argc, argv, "-h") )
+		exit_with_help();
+	parse_command_line(argc, argv, model, input_filename, output_filename,
+			model_filename);
+
+	// TODO: make sure that read_data allows for files without labels
+	read_data(data, input_filename);
+	read_model(model, model_filename);
+
+	// check if the number of attributes in data equals that in model
+	if (data->m != model->m) {
+		fprintf(stderr, "Error: number of attributes in data (%li) "
+				"does not equal the number of attributes in "
+				"model (%li)\n", data->m, model->m);
+		exit(1);
+	}
+
+	predy = Calloc(long, data->n);
+	predict_labels(data, model, predy);
+	if (data->y != NULL) {
+		performance = prediction_perf(data, predy);
+		info("Predictive performance: %3.2f%%\n", performance);
+	}
+
+	if (check_argv_eq(argc, argv, "-o")) {
+		write_predictions(data, predy, output_filename);
+		info("Predictions written to: %s\n", output_filename);
+	}
+
+	free_model(model);
+	free_data(data);
+	free(predy);
+
+	return 0;
+}
+
+void parse_command_line(int argc, char **argv, struct Model *model, 
+		char *input_filename, char *output_filename, char *model_filename)
+{
+	int i;
+	void (*print_func)(const char*) = NULL;	
+
+	for (i=1; i<argc; i++) {
+		if (argv[i][0] != '-') break;
+		if (++i >= argc)
+			exit_with_help();
+		switch (argv[i-1][1]) {
+			case 'o':
+				strcpy(output_filename, argv[i]);
+				break;
+			case 'q':
+				print_func = &print_null;
+				i--;
+			default:
+				fprintf(stderr, "Unknown option: -%c\n", argv[i-1][1]);
+				exit_with_help();
+		}
+	}
+
+	set_print_string_function(print_func);
+
+	if (i >= argc)
+		exit_with_help();
+
+	strcpy(input_filename, argv[i]);
+	i++;
+	strcpy(model_filename, argv[i]);
+}