Added mean region width

Author: khoroshevskyi

bindings/python/gtars/models/__init__.pyi

@@ -82,6 +82,12 @@ class RegionSet:
         """
         ...
 
+    def mean_region_width(self) -> int:
+        """
+        Mean width of the regions
+        """
+        ...
+
     def __len__(self) -> int:
         """
         Size of the regionset

bindings/python/src/models/region_set.rs

@@ -157,6 +157,10 @@ impl PyRegionSet {
         self.regionset.sort();
         Ok(())
     }
+
+    fn mean_region_width(&self) -> PyResult<u32> {
+        Ok(self.regionset.mean_region_width())
+    }
 }
 
 #[pyclass(name = "TokenizedRegionSet", module = "gtars.models")]

gtars/src/common/models/region_set.rs

@@ -50,7 +50,8 @@ impl TryFrom<&Path> for RegionSet {
 
         let reader = match path.is_file() {
             true => get_dynamic_reader(path).expect("!Can't read file"),
-            false => get_dynamic_reader_from_url(path).expect("!Can't get file neither from path or url!")
+            false => get_dynamic_reader_from_url(path)
+                .expect("!Can't get file neither from path or url!"),
         };
 
         let mut header: String = String::new();
@@ -368,6 +369,20 @@ impl RegionSet {
         false
     }
 
+    pub fn mean_region_width(&self) -> u32 {
+        if self.is_empty() {
+            return 0;
+        }
+        let sum: u32 = self
+            .regions
+            .iter()
+            .map(|region| region.end - region.start)
+            .sum();
+        let count: u32 = self.regions.len() as u32;
+
+        sum / count
+    }
+
     ///
     /// Get number of regions in RegionSet
     ///