/*
written by Patrick E. Meyer under the supervision of prof. Manolis Kellis
UniversitŽ Libre de Bruxelles - Massachusetts Institute of Technology
2010-2011
This file is part of Modencodenet.
Modencodenet is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
Modencodenet is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with Modencodenet. If not, see .
*/
#include "fly.h"
int main(int argc, char ** argv) {
vector temp;
char a[100];
const int n = 7;
int m = 10211202;
double** values = new double*[n];
for(int i = 0; i < n; ++i)
values[i] = new double[m];
for(int i = 0; i < n; ++i) {
read_vector(argv[i+1],temp);
values[i] = vectorString_to_vectorDouble(temp);
cout << i << endl;
}
double w[n] = {2.2950,1.0125,2.049,1.109,0.9751,0.3797,0.4526}; //parameters learned with R using CV on REDFly data
ofstream ofs(argv[n+1]);
double mod;
for(int i = 0; i < (m-1); ++i) {
mod= -1.5790;//-4.7142; //parameters learned with R
for(int j = 0; j