gam_models/fit_orchestration/
fit_config.rs1use super::*;
2
3fn normalize_optional_column(value: Option<String>, field: &str) -> Result<Option<String>, String> {
4 value
5 .map(|value| {
6 let value = value.trim();
7 if value.is_empty() {
8 Err(format!("{field} must be a non-empty column name"))
9 } else {
10 Ok(value.to_string())
11 }
12 })
13 .transpose()
14}
15
16pub fn validate_survival_baseline_config(
18 likelihood_mode: SurvivalLikelihoodMode,
19 baseline_target: &str,
20 baseline_scale: Option<f64>,
21 baseline_shape: Option<f64>,
22 baseline_rate: Option<f64>,
23 baseline_makeham: Option<f64>,
24) -> Result<(), String> {
25 if likelihood_mode == SurvivalLikelihoodMode::Weibull {
26 if baseline_rate.is_some() || baseline_makeham.is_some() {
27 return Err(
28 "survival likelihood 'weibull' does not use baseline_rate or baseline_makeham"
29 .to_string(),
30 );
31 }
32 if !matches!(baseline_target, "linear" | "weibull") {
33 return Err(
34 "survival likelihood 'weibull' supports only baseline_target 'linear' or 'weibull'"
35 .to_string(),
36 );
37 }
38 return Ok(());
39 }
40
41 match baseline_target {
42 "linear" => {
43 if baseline_scale.is_some()
44 || baseline_shape.is_some()
45 || baseline_rate.is_some()
46 || baseline_makeham.is_some()
47 {
48 return Err("baseline_target 'linear' does not use baseline parameters".to_string());
49 }
50 }
51 "weibull" => {
52 if baseline_rate.is_some() || baseline_makeham.is_some() {
53 return Err(
54 "baseline_target 'weibull' does not use baseline_rate or baseline_makeham"
55 .to_string(),
56 );
57 }
58 }
59 "gompertz" => {
60 if baseline_scale.is_some() || baseline_makeham.is_some() {
61 return Err(
62 "baseline_target 'gompertz' does not use baseline_scale or baseline_makeham"
63 .to_string(),
64 );
65 }
66 }
67 "gompertz-makeham" => {
68 if baseline_scale.is_some() {
69 return Err(
70 "baseline_target 'gompertz-makeham' does not use baseline_scale".to_string(),
71 );
72 }
73 }
74 other => {
75 return Err(format!(
76 "unsupported baseline_target '{other}'; use linear, weibull, gompertz, or gompertz-makeham"
77 ));
78 }
79 }
80 Ok(())
81}
82
83impl FitConfig {
84 pub fn resolve(mut self) -> Result<Self, String> {
90 self.family = self.family.and_then(|value| {
91 let value = value.trim();
92 (!value.eq_ignore_ascii_case("auto")).then(|| value.to_string())
93 });
94 self.survival_likelihood = self.survival_likelihood.trim().to_ascii_lowercase();
95 self.baseline_target = self.baseline_target.trim().to_ascii_lowercase();
96 self.link = self.link.and_then(|value| {
97 let value = value.trim();
98 (!value.is_empty()).then(|| value.to_string())
99 });
100 self.offset_column = normalize_optional_column(self.offset_column, "offset_column")?;
101 self.noise_offset_column =
102 normalize_optional_column(self.noise_offset_column, "noise_offset_column")?;
103 self.weight_column = normalize_optional_column(self.weight_column, "weight_column")?;
104 self.z_column = normalize_optional_column(self.z_column, "z_column")?;
105
106 if !self.ridge_lambda.is_finite() || self.ridge_lambda < 0.0 {
107 return Err("ridge_lambda must be finite and >= 0".to_string());
108 }
109 if self.outer_max_iter == Some(0) {
110 return Err("outer_max_iter must be >= 1".to_string());
111 }
112 self.frailty.validate().map_err(|error| error.to_string())?;
113 self.spatial_optimization.validate()?;
114 let likelihood_mode = parse_survival_likelihood_mode(&self.survival_likelihood)?;
115 validate_survival_baseline_config(
116 likelihood_mode,
117 &self.baseline_target,
118 self.baseline_scale,
119 self.baseline_shape,
120 self.baseline_rate,
121 self.baseline_makeham,
122 )?;
123 Ok(self)
124 }
125}
126
127#[cfg(test)]
128mod tests {
129 use super::*;
130
131 #[test]
132 fn resolve_normalizes_front_end_spellings() {
133 let resolved = FitConfig {
134 family: Some(" AUTO ".to_string()),
135 survival_likelihood: " Transformation ".to_string(),
136 baseline_target: " Linear ".to_string(),
137 ..FitConfig::default()
138 }
139 .resolve()
140 .unwrap();
141 assert_eq!(resolved.family, None);
142 assert_eq!(resolved.survival_likelihood, "transformation");
143 assert_eq!(resolved.baseline_target, "linear");
144 }
145
146 #[test]
147 fn resolve_rejects_invalid_shared_fields() {
148 assert!(
149 FitConfig {
150 ridge_lambda: f64::NAN,
151 ..FitConfig::default()
152 }
153 .resolve()
154 .is_err()
155 );
156 assert!(
157 FitConfig {
158 outer_max_iter: Some(0),
159 ..FitConfig::default()
160 }
161 .resolve()
162 .is_err()
163 );
164 assert!(
165 FitConfig {
166 weight_column: Some(" ".to_string()),
167 ..FitConfig::default()
168 }
169 .resolve()
170 .is_err()
171 );
172 }
173}