| 135 | } |
| 136 | |
| 137 | float get_medianNotNA(const std::vector<float> & myvector) { |
| 138 | vector<float> data (myvector); |
| 139 | |
| 140 | data.erase(std::remove(data.begin(), data.end(), NA), data.end()); |
| 141 | |
| 142 | int ndatapoints = data.size(); |
| 143 | if (ndatapoints==0) { |
| 144 | cerr << "Error: zero values to calculate medians..\n"; |
| 145 | exit(-1); |
| 146 | } |
| 147 | sort(data.begin(),data.end()); |
| 148 | // Get median |
| 149 | float median_value = ndatapoints % 2 == 1 ? data[(ndatapoints-1)/2] : (data[ndatapoints/2 - 1] + data[ndatapoints/2])/float(2.0); |
| 150 | data.clear(); |
| 151 | return median_value; |
| 152 | } |
| 153 | |
| 154 | // Calculate sd across individuals around given mean |
| 155 | // --------------------------------------------------------------------------- |