as.raw(Y$SEQ), as.double(Y$ANC), as.double(Y$w),
as.integer(phy$edge[, 1]), as.integer(phy$edge[, 2]),
as.double(phy$edge.length), as.integer(npart),
as.raw(Y$SEQ), as.double(Y$ANC), as.double(Y$w),
as.integer(phy$edge[, 1]), as.integer(phy$edge[, 2]),
as.double(phy$edge.length), as.integer(npart),