You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Simplifying the tree sequence down to a specific subset of individuals and then loading it doesn't work. Boo.
> ts %>% ts_simplify(simplify_to = c("pop_1", "pop_2", "pop_3")) %>% ts_save("/tmp/small2.trees")
> ts_load("/tmp/small2.trees", model)
Error in `dplyr::bind_cols()`:
! Can't recycle `..1` (size 3) to match `..2` (size 100).
Run `rlang::last_trace()` to see where the error occurred.
The error suggests an issue with some join-like operation on some tree sequence table, although I'd think the tree sequence itself should carry every piece of slendr-relevant metadata along side the standard tskit data? Well, clearly not.
Interestingly, loading the small tree sequence without the slendr metadata works without problems. So something is clearly messed up with slendr-specific information not being passed through when a simplified tree sequence is saved.
Yeah, the population had 100 individuals to start with. Simple simplification does not change the number of "sampled individuals" (still 100). But the subsetting of individuals changes the number of sampled individuals (in the code above down to 3), which matches the error from the join operation:
! Can't recycle ..1 (size 3) to match ..2 (size 100).
Let's encode the following trivial model:
Simulate a tree sequence from it, save that tree sequence to disk:
So far so good.
The error suggests an issue with some join-like operation on some tree sequence table, although I'd think the tree sequence itself should carry every piece of slendr-relevant metadata along side the standard tskit data? Well, clearly not.
The text was updated successfully, but these errors were encountered: