Changeset 83


Ignore:
Timestamp:
May 26, 2004, 4:42:25 PM (19 years ago)
Author:
cecilia
Message:

added test output

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/Merge.cc

    r82 r83  
    4141    weights = calculate_weights();
    4242
    43     // för varje annot: kolla vilka rader den finns på
     43    // för varje annot: kolla vilka rader den finns på
    4444    map<string, vector<u_int> > name;
    4545    for (u_int i=0; i<annots_.size(); i++)
    4646      name[annots_[i]].push_back(i);
    4747   
    48     // för varje namn:
    49    
     48    // för varje annot
    5049    for (map<string,vector<u_int> >::const_iterator k=name.begin(); k!=name.end(); k++)
    5150      merging((*k).second, weights);
     
    5352
    5453  void Merge::merging(vector<u_int> index, gslapi::matrix& weights){
    55     // nya rader i error_ och data_
    56    
    57     vector<double> new_data;
    58     vector<double> new_error;
     54
    5955    // merge genes for each assay
    6056    for (u_int j=0; j<nof_assays_; j++){
     
    6864      }
    6965      double m = sum_wx/sum_w;
    70       new_data.push_back(m);
    7166     
    7267      // calc. U
     
    8075      }
    8176      double U = 1/sum_u+sum_var/(sum_w*sum_w);
    82       new_error.push_back(U);
    8377     
    8478      // insert new values
Note: See TracChangeset for help on using the changeset viewer.