kin_multi_stage_time_variant_2sex.R needs updating to reflect Focal's fathers multi-state distribution.

[ButterickJoe]

Hi Ivan,

I have recently been working on UK parity specific kinship with two sexes, and noticed that the DemoKin code I provided you does not correctly implement the age*stage distribution of Focal's father. As such I have adjusted the file kin_multi_stage_time_variant_2sex.R. Note that the only kin affected by the previous error is Focal’s father: no other kin within Focal’s network are affected.

The changes in the file can be found in functions all_kin_dy and all_kin_dy_TV, where in the matrix projections, I have changed X_parents[,1] (i.e., the distribution of Focal's parents in her first age-class) to reflect both mother and father (lines 400 and 618), moving from the input: parents_joint_age_stage <- population_age_stage_of_parenting[[2]] to parents_joint_age_stage <- population_age_stage_of_parenting[[1]] (lines 342, 354, 562, and 575).

In doing so I noticed that trying to join two sparse vectors using c(.) does not work so have moved to rbind(). These changes occur in the functions pi_mix, pi_mix_TV, pi_mix_parity and pi_mix_TV_parity (lines 887, 915, 958, and 997 resp.).

Sorry for this. I guess sometimes we find mistakes in code which we do not expect! I am sending this request with an updated version of kin_multi_stage_time_variant_2sex.R. The only changes involve what I mention above.

Hope that you are well.

Joe

[IvanWilli]

Hi Joe
Good detection. Yes, because fertility is linked always to woman, then father´s initial distribution dont have an impact.
Dont worry! we did a lot of correction ealier in other functions ;). And also i did not dectect that when revised the function.
Thanks!