@@ -177,41 +177,28 @@ def create_time_lineage_fate_correlation_plot(
177
177
178
178
if dataset_name == "larry_multilineage" :
179
179
logger .info (
180
- "Creating multilineage clone trajectory from mono and neu subsets"
180
+ "Creating multilineage clone trajectory from mono and neu datasets"
181
+ )
182
+
183
+ from pyrovelocity .io .datasets import larry_mono , larry_neu
184
+
185
+ mono_adata = larry_mono ()
186
+ neu_adata = larry_neu ()
187
+
188
+ logger .info (
189
+ f" - Generating mono trajectory with { mono_adata .n_obs } cells"
190
+ )
191
+ mono_clone = get_clone_trajectory (mono_adata )
192
+
193
+ logger .info (
194
+ f" - Generating neu trajectory with { neu_adata .n_obs } cells"
195
+ )
196
+ neu_clone = get_clone_trajectory (neu_adata )
197
+
198
+ logger .info (" - Concatenating mono and neu trajectories" )
199
+ clone_trajectories [dataset_name ] = mono_clone .concatenate (
200
+ neu_clone
181
201
)
182
-
183
- if "state_info" in adata_pyrovelocity .obs :
184
- mono_mask = adata_pyrovelocity .obs ["state_info" ].str .contains (
185
- "Mono" , case = False , na = False
186
- )
187
- neu_mask = adata_pyrovelocity .obs ["state_info" ].str .contains (
188
- "Neu" , case = False , na = False
189
- )
190
-
191
- mono_adata = adata_pyrovelocity [mono_mask ].copy ()
192
- neu_adata = adata_pyrovelocity [neu_mask ].copy ()
193
-
194
- logger .info (
195
- f" - Generating mono trajectory with { mono_adata .n_obs } cells"
196
- )
197
- mono_clone = get_clone_trajectory (mono_adata )
198
-
199
- logger .info (
200
- f" - Generating neu trajectory with { neu_adata .n_obs } cells"
201
- )
202
- neu_clone = get_clone_trajectory (neu_adata )
203
-
204
- logger .info (" - Concatenating mono and neu trajectories" )
205
- clone_trajectories [dataset_name ] = mono_clone .concatenate (
206
- neu_clone
207
- )
208
- else :
209
- logger .warning (
210
- "Could not identify mono/neu cells in multilineage dataset. Generating unified trajectory."
211
- )
212
- clone_trajectories [dataset_name ] = get_clone_trajectory (
213
- adata_pyrovelocity
214
- )
215
202
else :
216
203
logger .info (
217
204
f"Generating clone trajectory for { dataset_name } with { adata_pyrovelocity .n_obs } cells"
0 commit comments