Implement model diagnostics plot for a single feature

# Description

Create a function plot_model_diagnostics(gene_model) to visualize:

- Observed vs. fitted values
- Q-Q plot of residuals
- Histogram with overlaid fitted distribution

This aids in evaluating model quality on a per-feature basis.