Trying to write a simple neutral nucleotide simulation and extract the sequences at the end. The simulation runs fine when the output function is omitted. The example code in the R docs throws the same error.
Simulation finished with exit status: 1
ERROR (EidosValue_Object::GetPropertyOfElements): property genome1 is not defined for object element type Individual.
Error on script line 17, character 49:
cat(paste(sim.subpopulations.individuals.genome1.nucleotides()));
^^^^^^^
simple.3.pops.mig <- slim_script(
slim_block(initialize(), {
initializeSLiMOptions(nucleotideBased=T);
initializeAncestralNucleotides(randomNucleotides(1000));
initializeMutationTypeNuc("m1", 0.5, "f", 0.0);
initializeGenomicElementType("g1", m1, 1.0, mmJukesCantor(1e-5));
initializeGenomicElement(g1, 0, 1000 - 1);
initializeRecombinationRate(1e-8);
}),
slim_block(1, early(), {
sim.addSubpop("p1", 1000);
sim.addSubpop("p2", 1000);
sim.addSubpop("p3", 1000);
p1.setMigrationRates(c(p2,p3), c(0.01,0.01));
p2.setMigrationRates(c(p1,p3), c(0.01,0.2));
p3.setMigrationRates(c(p1,p2), c(0.01,0.2));
}),
slim_block(2000, late(), {
r_output_nucleotides(subpops = T);
sim.simulationFinished();
})
)
Trying to write a simple neutral nucleotide simulation and extract the sequences at the end. The simulation runs fine when the output function is omitted. The example code in the R docs throws the same error.
The error:
my code:
Thanks!