default_priors() fails if the input strphylo does not contain a migration history

[IanPRoberts]

library( StructCoalescent ); treedata <- treeio::read.beast.newick(file = 'tree.newick.txt'); try( default_priors( as.strphylo( treedata ), n_deme = 2 ) )
fails with error
Error in DemeDecompC(ED, max(ED[, 5]), NodeIndicesC(ED)) : Not a matrix.

tree.newick.txt

[IanPRoberts]

Temporary Workaround

Error is caused by calling default_priors() without including locations and types for migration events. Sampling a single migration history to generate a full structured phylogeny avoids the issue

strphylo <- as.strphylo( treedata )
fit_mig_mat <- matrix(c(0,1,1,0), 2, 2)
strphylo <- DTA_sampling( strphylo, fit_mig_mat, 1 )
try( default_priors( strphylo, n_deme = 2 ) )

passes default_priors()