update readme

README.md

@@ -159,6 +159,15 @@ p <- diagnostic_plot(list_tmb = l_tmb)
 
 The diagnostic metrics show that the fit is not as good for all genes. 
 
+### Extracts a tidy result table from a list tmb object
+
+The tidy_results function extracts a dataframe containing estimates of ln(fold changes), standard errors, test statistics, p-values, and adjusted p-values for fixed effects. Additionally, it provides access to correlation terms and standard deviations for random effects, offering a detailed view of HTRfit modeling results.
+
+```
+## -- get tidy results
+tidy_res <- tidy_results(l_tmb, coeff_threshold = 0.1, alternative_hypothesis = "greater")
+```
+
 ### Evaluation 
 
 ```