The only part of this function differring from my previous functional one is adding a variable to be calculated in each generation. NDeme=0.5(Nmax+Nmin)+0.5(Nmax+Nmin)*sin(2pi*t/period) #period=12 genetation in setting #oscillating population size As NDeme can be 0 at 9th (9+12n) generation, I set…

During the process of making minimum working example, I found it’s due to the lowest point from this sin() function. After migration, the NDeme (sometimes serve as denominator) could be 0, which leads to NaNs in output.

Output matrix in last (t-1) generation: normal; Indexing the same matrix in this generation: NaNs in the (t-1) column

General Usage

Tamas_Papp February 24, 2020, 6:35am 2

I think there is a lot of context missing here, please post a self-contained MWE.

1 Like