Quick Example of Simulated Kinships with Partial Parentage

Introduction

This vignette provides a quick example of how to estimate kinship coefficients using simulation. One simulation uses 100 iterations and the other uses the same setup with 500 simulations to illustrate the type of kinship coefficient estimate variance you can expect with a simple pedigree exhibiting partial parentage for some of the members.

Anticipation of further uses of this kinship estimation method

Retention of founder alleles is a primary driving force behind this package of utilities. This example is intentionally overly simplistic to clearly illustrate this one aspect of the methodology used with realistic pedigrees.

In this example, some of the potential parents are not members of the example pedigree and thus, by definition, have a kinship coeficient of 0.0 with all other pedigree members.

In practice most or all of the potential parents will be members of the same pedigree and will potentially have nonzero kinship coeficients with other pedigree members.

Identification of potential parents

This tutorial assumes knowledge of potential parents and does not present methods for identifying potential parents.

Planned capabilities¹ will provide users the ability to fully automate the identification of potential parents. However, prior to that capability being available, users can use other methods to identify potential parents.

Creation of example pedigree²

The example is based on the following simple pedigree setup. In the pedigree given below, all of the original animals have identifiers that are single letters of the alphabet.:

Only those IDs that have unknown parents are included in the lists. For those IDs with one known parent, the known parent is included in the . See for ID . Note also that potential parents can either come from the pedigree being augmented by simulation or from outside the pedigree. See , , , and for examples of this.

Currently, if you want to provide differential weights for the different potential parents, you will need to do this directly by modifying the number of times each parent is included in the list.

Animals , , and

knitr::opts_chunk$set(echo = TRUE)
library(kableExtra) # nolint: undesirable_function_linter
library(nprcgenekeepr) # nolint: undesirable_function_linter
library(stringi) # nolint: undesirable_function_linter

ped <- nprcgenekeepr::smallPed
simParent_1 <- list( # nolint: object_name_linter
  id = "A", # nolint: object_name_linter
  sires = "Q",
  dams = c("d1_1", "d1_2", "d1_3", "d1_4")
)
simParent_2 <- list( # nolint: object_name_linter
  id = "B", # nolint: object_name_linter
  sires = c("s1_1", "s1_2", "s1_3"),
  dams = c("d1_1", "d1_2", "d1_3", "d1_4")
)
simParent_3 <- list( # nolint: object_name_linter
  id = "E", # nolint: object_name_linter
  sires = c("A", "C", "s1_1"),
  dams = c("d3_1", "B")
)
simParent_4 <- list( # nolint: object_name_linter
  id = "J", # nolint: object_name_linter
  sires = c("A", "C", "s1_1"),
  dams = c("d3_1", "B")
)
simParent_5 <- list( # nolint: object_name_linter
  id = "K", # nolint: object_name_linter
  sires = c("A", "C", "s1_1"),
  dams = c("d3_1", "d1_2")
)
simParent_6 <- list( # nolint: object_name_linter
  id = "N", # nolint: object_name_linter
  sires = c("A", "C", "s1_2"),
  dams = c("d3_1", "B")
)
allSimParents <- list(
  simParent_1, simParent_2, simParent_3,
  simParent_4, simParent_5, simParent_6
)

extractKinship <- function(simKinships, id1, id2, simulation) {
  ids <- dimnames(simKinships[[simulation]])[[1L]]
  simKinships[[simulation]][
    seq_along(ids)[ids == id1],
    seq_along(ids)[ids == id2]
  ]
}

extractKValue <- function(kValue, id1, id2, simulation) {
  kValue[kValue$id_1 == id1 & kValue$id_2 == id2, paste0("sim_", simulation)]
}

Small Example

This is the simulation. I am only printing out rows with kinship values that vary.

Before running these simulations, take time to look at the included function descriptions to see what they are expecting as arguments and what they return.

?createSimKinships
?kinshipMatricesToKValues
?extractKValue

# Only set this seed if you want to get the same simulation results each time.
set.seed(1L)
n <- 10L
simKinships <- createSimKinships(ped, allSimParents, pop = ped$id, n = n)
kValues <- kinshipMatricesToKValues(simKinships)
extractKValue(kValues, id1 = "A", id2 = "F", simulation = 1L:n)

##  [1] "sim_1"  "sim_2"  "sim_3"  "sim_4"  "sim_5"  "sim_6"  "sim_7"  "sim_8" 
##  [9] "sim_9"  "sim_10"

counts <- countKinshipValues(kValues)
counts$kinshipIds[1L:3L]

## NULL

counts$kinshipValues[1L:3L]

## NULL

counts$kinshipCounts[1L:3L]

## NULL

stats_10 <- summarizeKinshipValues(counts)
nrow(stats_10[stats_10$sd > 0.0, ])

## [1] 78

kable(stats_10[stats_10$sd > 0.0, ], longtable = TRUE) |>
  kable_styling(
    latex_options = c("striped", "repeat_header"),
    repeat_header_method = "replace",
    repeat_header_text = "\\textit{(continued)}"
  )

	id_1	id_2	min	secondQuartile	mean	median	thirdQuartile	max	sd
5	A	E	0.00000	0.00000	0.075000	0.000000	0.25000	0.25000	0.1207615
6	A	F	0.12500	0.12500	0.162500	0.125000	0.25000	0.25000	0.0603807
7	A	G	0.12500	0.12500	0.162500	0.125000	0.25000	0.25000	0.0603807
9	A	I	0.25000	0.25000	0.312500	0.312500	0.37500	0.37500	0.0658808
10	A	J	0.00000	0.00000	0.125000	0.125000	0.25000	0.25000	0.1317616
14	A	N	0.00000	0.00000	0.050000	0.000000	0.00000	0.25000	0.1054093
21	B	E	0.00000	0.00000	0.150000	0.250000	0.25000	0.25000	0.1290994
22	B	F	0.12500	0.12500	0.200000	0.250000	0.25000	0.25000	0.0645497
23	B	G	0.12500	0.12500	0.200000	0.250000	0.25000	0.25000	0.0645497
25	B	I	0.00000	0.00000	0.075000	0.125000	0.12500	0.12500	0.0645497
26	B	J	0.00000	0.00000	0.150000	0.250000	0.25000	0.25000	0.1290994
30	B	N	0.00000	0.00000	0.050000	0.000000	0.00000	0.25000	0.1054093
36	C	E	0.00000	0.00000	0.112500	0.125000	0.12500	0.25000	0.0922331
37	C	F	0.12500	0.12500	0.181250	0.187500	0.18750	0.25000	0.0461165
38	C	G	0.12500	0.12500	0.181250	0.187500	0.18750	0.25000	0.0461165
40	C	I	0.12500	0.18750	0.193750	0.187500	0.18750	0.25000	0.0354779
41	C	J	0.00000	0.12500	0.137500	0.125000	0.12500	0.25000	0.0709558
45	C	N	0.00000	0.00000	0.050000	0.000000	0.12500	0.25000	0.0874007
50	D	E	0.00000	0.00000	0.112500	0.125000	0.12500	0.25000	0.0922331
51	D	F	0.25000	0.25000	0.306250	0.312500	0.31250	0.37500	0.0461165
52	D	G	0.25000	0.25000	0.306250	0.312500	0.31250	0.37500	0.0461165
54	D	I	0.12500	0.18750	0.193750	0.187500	0.18750	0.25000	0.0354779
55	D	J	0.00000	0.12500	0.137500	0.125000	0.12500	0.25000	0.0709558
59	D	N	0.00000	0.00000	0.050000	0.000000	0.12500	0.25000	0.0874007
64	E	F	0.25000	0.25000	0.306250	0.312500	0.31250	0.37500	0.0461165
65	E	G	0.25000	0.25000	0.306250	0.312500	0.31250	0.37500	0.0461165
66	E	H	0.00000	0.00000	0.112500	0.125000	0.12500	0.25000	0.0922331
67	E	I	0.00000	0.00000	0.087500	0.062500	0.18750	0.18750	0.0790569
68	E	J	0.00000	0.00000	0.100000	0.125000	0.12500	0.25000	0.0790569
69	E	K	0.00000	0.00000	0.062500	0.062500	0.12500	0.12500	0.0658808
70	E	L	0.06250	0.06250	0.087500	0.062500	0.12500	0.12500	0.0322749
71	E	M	0.00000	0.00000	0.037500	0.000000	0.12500	0.12500	0.0603807
72	E	N	0.00000	0.00000	0.037500	0.000000	0.12500	0.12500	0.0603807
74	E	P	0.00000	0.00000	0.018750	0.000000	0.06250	0.06250	0.0301904
75	E	Q	0.00000	0.00000	0.037500	0.000000	0.12500	0.12500	0.0603807
76	F	F	0.50000	0.50000	0.556250	0.562500	0.56250	0.62500	0.0461165
77	F	G	0.25000	0.25000	0.306250	0.312500	0.31250	0.37500	0.0461165
78	F	H	0.12500	0.12500	0.181250	0.187500	0.18750	0.25000	0.0461165
79	F	I	0.06250	0.12500	0.140625	0.140625	0.18750	0.18750	0.0423127
80	F	J	0.00000	0.12500	0.118750	0.125000	0.12500	0.18750	0.0547247
81	F	K	0.00000	0.00000	0.031250	0.031250	0.06250	0.06250	0.0329404
82	F	L	0.09375	0.09375	0.106250	0.093750	0.12500	0.12500	0.0161374
83	F	M	0.06250	0.06250	0.081250	0.062500	0.12500	0.12500	0.0301904
84	F	N	0.00000	0.00000	0.043750	0.000000	0.12500	0.12500	0.0592927
86	F	P	0.03125	0.03125	0.040625	0.031250	0.06250	0.06250	0.0150952
87	F	Q	0.06250	0.06250	0.081250	0.062500	0.12500	0.12500	0.0301904
88	G	G	0.50000	0.50000	0.556250	0.562500	0.56250	0.62500	0.0461165
89	G	H	0.12500	0.12500	0.181250	0.187500	0.18750	0.25000	0.0461165
90	G	I	0.06250	0.12500	0.140625	0.140625	0.18750	0.18750	0.0423127
91	G	J	0.00000	0.12500	0.118750	0.125000	0.12500	0.18750	0.0547247
92	G	K	0.00000	0.00000	0.031250	0.031250	0.06250	0.06250	0.0329404
93	G	L	0.09375	0.09375	0.106250	0.093750	0.12500	0.12500	0.0161374
94	G	M	0.06250	0.06250	0.081250	0.062500	0.12500	0.12500	0.0301904
95	G	N	0.00000	0.00000	0.043750	0.000000	0.12500	0.12500	0.0592927
97	G	P	0.03125	0.03125	0.040625	0.031250	0.06250	0.06250	0.0150952
98	G	Q	0.06250	0.06250	0.081250	0.062500	0.12500	0.12500	0.0301904
100	H	I	0.12500	0.18750	0.193750	0.187500	0.18750	0.25000	0.0354779
101	H	J	0.00000	0.12500	0.137500	0.125000	0.12500	0.25000	0.0709558
105	H	N	0.00000	0.00000	0.050000	0.000000	0.12500	0.25000	0.0874007
109	I	I	0.50000	0.50000	0.562500	0.562500	0.62500	0.62500	0.0658808
110	I	J	0.25000	0.25000	0.312500	0.312500	0.37500	0.37500	0.0658808
111	I	K	0.00000	0.00000	0.006250	0.000000	0.00000	0.06250	0.0197642
112	I	L	0.06250	0.09375	0.100000	0.093750	0.12500	0.12500	0.0197642
113	I	M	0.12500	0.12500	0.156250	0.156250	0.18750	0.18750	0.0329404
114	I	N	0.00000	0.00000	0.056250	0.000000	0.12500	0.25000	0.0952427
116	I	P	0.06250	0.06250	0.078125	0.078125	0.09375	0.09375	0.0164702
117	I	Q	0.12500	0.12500	0.156250	0.156250	0.18750	0.18750	0.0329404
119	J	K	0.00000	0.00000	0.012500	0.000000	0.00000	0.12500	0.0395285
120	J	L	0.00000	0.06250	0.075000	0.062500	0.12500	0.12500	0.0395285
121	J	M	0.00000	0.00000	0.062500	0.062500	0.12500	0.12500	0.0658808
122	J	N	0.00000	0.00000	0.062500	0.000000	0.12500	0.25000	0.1062296
124	J	P	0.00000	0.00000	0.031250	0.031250	0.06250	0.06250	0.0329404
125	J	Q	0.00000	0.00000	0.062500	0.062500	0.12500	0.12500	0.0658808
129	K	N	0.00000	0.00000	0.025000	0.000000	0.00000	0.12500	0.0527046
135	L	N	0.00000	0.00000	0.037500	0.000000	0.06250	0.12500	0.0527046
140	M	N	0.00000	0.00000	0.025000	0.000000	0.00000	0.12500	0.0527046
146	N	P	0.00000	0.00000	0.012500	0.000000	0.00000	0.06250	0.0263523
147	N	Q	0.00000	0.00000	0.025000	0.000000	0.00000	0.12500	0.0527046

A larger simulation

set.seed(1L)
n <- 100L
simKinships <- createSimKinships(ped, allSimParents, pop = ped$id, n = n)
kValues <- kinshipMatricesToKValues(simKinships)
extractKValue(kValues, id1 = "A", id2 = "F", simulation = 1L:10L)

##  [1] "sim_1"  "sim_2"  "sim_3"  "sim_4"  "sim_5"  "sim_6"  "sim_7"  "sim_8" 
##  [9] "sim_9"  "sim_10"

counts <- countKinshipValues(kValues)
counts$kinshipIds[1L:3L]

## NULL

counts$kinshipValues[1L:3L]

## NULL

counts$kinshipCounts[1L:3L]

## NULL

stats_100 <- summarizeKinshipValues(counts)
nrow(stats_100[stats_100$sd > 0.0, ])

## [1] 94

kable(stats_100[stats_100$sd > 0.0, ], longtable = TRUE) |>
  kable_styling(
    latex_options = c("striped", "repeat_header"),
    repeat_header_method = "replace",
    repeat_header_text = "\\textit{(continued)}"
  )

	id_1	id_2	min	secondQuartile	mean	median	thirdQuartile	max	sd
5	A	E	0.00000	0.00000	0.0825000	0.00000	0.25000	0.25000	0.1181454
6	A	F	0.12500	0.12500	0.1662500	0.12500	0.25000	0.25000	0.0590727
7	A	G	0.12500	0.12500	0.1662500	0.12500	0.25000	0.25000	0.0590727
9	A	I	0.25000	0.25000	0.2987500	0.25000	0.37500	0.37500	0.0612759
10	A	J	0.00000	0.00000	0.0975000	0.00000	0.25000	0.25000	0.1225518
11	A	K	0.00000	0.00000	0.0700000	0.00000	0.25000	0.25000	0.1128152
12	A	L	0.12500	0.12500	0.1600000	0.12500	0.25000	0.25000	0.0564076
14	A	N	0.00000	0.00000	0.0800000	0.00000	0.25000	0.25000	0.1172065
21	B	E	0.00000	0.00000	0.1325000	0.25000	0.25000	0.25000	0.1254034
22	B	F	0.12500	0.12500	0.1912500	0.25000	0.25000	0.25000	0.0627017
23	B	G	0.12500	0.12500	0.1912500	0.25000	0.25000	0.25000	0.0627017
25	B	I	0.00000	0.00000	0.0625000	0.06250	0.12500	0.12500	0.0628149
26	B	J	0.00000	0.00000	0.1250000	0.12500	0.25000	0.25000	0.1256297
30	B	N	0.00000	0.00000	0.1100000	0.00000	0.25000	0.25000	0.1247219
36	C	E	0.00000	0.00000	0.1075000	0.12500	0.12500	0.25000	0.0870751
37	C	F	0.12500	0.12500	0.1787500	0.18750	0.18750	0.25000	0.0435375
38	C	G	0.12500	0.12500	0.1787500	0.18750	0.18750	0.25000	0.0435375
40	C	I	0.12500	0.15625	0.1806250	0.18750	0.18750	0.25000	0.0386145
41	C	J	0.00000	0.06250	0.1112500	0.12500	0.12500	0.25000	0.0772291
42	C	K	0.00000	0.00000	0.0350000	0.00000	0.12500	0.12500	0.0564076
43	C	L	0.25000	0.25000	0.2675000	0.25000	0.31250	0.31250	0.0282038
45	C	N	0.00000	0.00000	0.0950000	0.12500	0.12500	0.25000	0.0925235
50	D	E	0.00000	0.00000	0.1075000	0.12500	0.12500	0.25000	0.0870751
51	D	F	0.25000	0.25000	0.3037500	0.31250	0.31250	0.37500	0.0435375
52	D	G	0.25000	0.25000	0.3037500	0.31250	0.31250	0.37500	0.0435375
54	D	I	0.12500	0.15625	0.1806250	0.18750	0.18750	0.25000	0.0386145
55	D	J	0.00000	0.06250	0.1112500	0.12500	0.12500	0.25000	0.0772291
56	D	K	0.00000	0.00000	0.0350000	0.00000	0.12500	0.12500	0.0564076
57	D	L	0.12500	0.12500	0.1425000	0.12500	0.18750	0.18750	0.0282038
59	D	N	0.00000	0.00000	0.0950000	0.12500	0.12500	0.25000	0.0925235
64	E	F	0.25000	0.25000	0.3037500	0.31250	0.31250	0.37500	0.0435375
65	E	G	0.25000	0.25000	0.3037500	0.31250	0.31250	0.37500	0.0435375
66	E	H	0.00000	0.00000	0.1075000	0.12500	0.12500	0.25000	0.0870751
67	E	I	0.00000	0.00000	0.0762500	0.06250	0.12500	0.25000	0.0730275
68	E	J	0.00000	0.00000	0.0700000	0.00000	0.12500	0.25000	0.0800883
69	E	K	0.00000	0.00000	0.0387500	0.00000	0.12500	0.12500	0.0581029
70	E	L	0.00000	0.06250	0.0731250	0.06250	0.12500	0.18750	0.0533409
71	E	M	0.00000	0.00000	0.0412500	0.00000	0.12500	0.12500	0.0590727
72	E	N	0.00000	0.00000	0.0625000	0.00000	0.12500	0.25000	0.0743235
74	E	P	0.00000	0.00000	0.0206250	0.00000	0.06250	0.06250	0.0295363
75	E	Q	0.00000	0.00000	0.0412500	0.00000	0.12500	0.12500	0.0590727
76	F	F	0.50000	0.50000	0.5537500	0.56250	0.56250	0.62500	0.0435375
77	F	G	0.25000	0.25000	0.3037500	0.31250	0.31250	0.37500	0.0435375
78	F	H	0.12500	0.12500	0.1787500	0.18750	0.18750	0.25000	0.0435375
79	F	I	0.06250	0.09375	0.1284375	0.12500	0.15625	0.25000	0.0425788
80	F	J	0.00000	0.06250	0.0906250	0.06250	0.12500	0.25000	0.0636727
81	F	K	0.00000	0.00000	0.0368750	0.00000	0.06250	0.12500	0.0436055
82	F	L	0.06250	0.09375	0.1078125	0.09375	0.12500	0.18750	0.0321447
83	F	M	0.06250	0.06250	0.0831250	0.06250	0.12500	0.12500	0.0295363
84	F	N	0.00000	0.00000	0.0787500	0.06250	0.12500	0.25000	0.0719607
86	F	P	0.03125	0.03125	0.0415625	0.03125	0.06250	0.06250	0.0147682
87	F	Q	0.06250	0.06250	0.0831250	0.06250	0.12500	0.12500	0.0295363
88	G	G	0.50000	0.50000	0.5537500	0.56250	0.56250	0.62500	0.0435375
89	G	H	0.12500	0.12500	0.1787500	0.18750	0.18750	0.25000	0.0435375
90	G	I	0.06250	0.09375	0.1284375	0.12500	0.15625	0.25000	0.0425788
91	G	J	0.00000	0.06250	0.0906250	0.06250	0.12500	0.25000	0.0636727
92	G	K	0.00000	0.00000	0.0368750	0.00000	0.06250	0.12500	0.0436055
93	G	L	0.06250	0.09375	0.1078125	0.09375	0.12500	0.18750	0.0321447
94	G	M	0.06250	0.06250	0.0831250	0.06250	0.12500	0.12500	0.0295363
95	G	N	0.00000	0.00000	0.0787500	0.06250	0.12500	0.25000	0.0719607
97	G	P	0.03125	0.03125	0.0415625	0.03125	0.06250	0.06250	0.0147682
98	G	Q	0.06250	0.06250	0.0831250	0.06250	0.12500	0.12500	0.0295363
100	H	I	0.12500	0.15625	0.1806250	0.18750	0.18750	0.25000	0.0386145
101	H	J	0.00000	0.06250	0.1112500	0.12500	0.12500	0.25000	0.0772291
102	H	K	0.00000	0.00000	0.0350000	0.00000	0.12500	0.12500	0.0564076
103	H	L	0.12500	0.12500	0.1425000	0.12500	0.18750	0.18750	0.0282038
105	H	N	0.00000	0.00000	0.0950000	0.12500	0.12500	0.25000	0.0925235
109	I	I	0.50000	0.50000	0.5487500	0.50000	0.62500	0.62500	0.0612759
110	I	J	0.25000	0.25000	0.2987500	0.25000	0.37500	0.37500	0.0612759
111	I	K	0.00000	0.00000	0.0468750	0.00000	0.12500	0.18750	0.0649002
112	I	L	0.06250	0.09375	0.1137500	0.09375	0.12500	0.21875	0.0359529
113	I	M	0.12500	0.12500	0.1493750	0.12500	0.18750	0.18750	0.0306379
114	I	N	0.00000	0.00000	0.0637500	0.00000	0.12500	0.25000	0.0794453
116	I	P	0.06250	0.06250	0.0746875	0.06250	0.09375	0.09375	0.0153190
117	I	Q	0.12500	0.12500	0.1493750	0.12500	0.18750	0.18750	0.0306379
119	J	K	0.00000	0.00000	0.0237500	0.00000	0.00000	0.12500	0.0492847
120	J	L	0.00000	0.06250	0.0675000	0.06250	0.12500	0.18750	0.0467370
121	J	M	0.00000	0.00000	0.0487500	0.00000	0.12500	0.12500	0.0612759
122	J	N	0.00000	0.00000	0.0475000	0.00000	0.12500	0.25000	0.0727785
124	J	P	0.00000	0.00000	0.0243750	0.00000	0.06250	0.06250	0.0306379
125	J	Q	0.00000	0.00000	0.0487500	0.00000	0.12500	0.12500	0.0612759
127	K	L	0.25000	0.25000	0.2675000	0.25000	0.31250	0.31250	0.0282038
128	K	M	0.00000	0.00000	0.0350000	0.00000	0.12500	0.12500	0.0564076
129	K	N	0.00000	0.00000	0.0287500	0.00000	0.00000	0.12500	0.0528691
131	K	P	0.00000	0.00000	0.0175000	0.00000	0.06250	0.06250	0.0282038
132	K	Q	0.00000	0.00000	0.0350000	0.00000	0.12500	0.12500	0.0564076
133	L	L	0.50000	0.50000	0.5175000	0.50000	0.56250	0.56250	0.0282038
134	L	M	0.06250	0.06250	0.0800000	0.06250	0.12500	0.12500	0.0282038
135	L	N	0.00000	0.00000	0.0618750	0.06250	0.12500	0.18750	0.0536654
137	L	P	0.03125	0.03125	0.0400000	0.03125	0.06250	0.06250	0.0141019
138	L	Q	0.06250	0.06250	0.0800000	0.06250	0.12500	0.12500	0.0282038
140	M	N	0.00000	0.00000	0.0400000	0.00000	0.12500	0.12500	0.0586033
146	N	P	0.00000	0.00000	0.0200000	0.00000	0.06250	0.06250	0.0293016
147	N	Q	0.00000	0.00000	0.0400000	0.00000	0.12500	0.12500	0.0586033

A much larger simulation

set.seed(1L)
n <- 500L
simKinships <- createSimKinships(ped, allSimParents, pop = ped$id, n = n)
kValues <- kinshipMatricesToKValues(simKinships)
extractKValue(kValues, id1 = "A", id2 = "F", simulation = 1L:10L)

##  [1] "sim_1"  "sim_2"  "sim_3"  "sim_4"  "sim_5"  "sim_6"  "sim_7"  "sim_8" 
##  [9] "sim_9"  "sim_10"

counts <- countKinshipValues(kValues)
counts$kinshipIds[1L:3L]

## NULL

counts$kinshipValues[1L:3L]

## NULL

counts$kinshipCounts[1L:3L]

## NULL

stats_500 <- summarizeKinshipValues(counts)
nrow(stats_500[stats_500$sd > 0.0, ])

## [1] 94

kable(stats_500[stats_500$sd > 0.0, ], longtable = TRUE) |>
  kable_styling(
    latex_options = c("striped", "repeat_header"),
    repeat_header_method = "replace",
    repeat_header_text = "\\textit{(continued)}"
  )

	id_1	id_2	min	secondQuartile	mean	median	thirdQuartile	max	sd
5	A	E	0.00000	0.00000	0.0865000	0.00000	0.25000	0.25000	0.1190424
6	A	F	0.12500	0.12500	0.1682500	0.12500	0.25000	0.25000	0.0595212
7	A	G	0.12500	0.12500	0.1682500	0.12500	0.25000	0.25000	0.0595212
9	A	I	0.25000	0.25000	0.2920000	0.25000	0.37500	0.37500	0.0591015
10	A	J	0.00000	0.00000	0.0840000	0.00000	0.25000	0.25000	0.1182030
11	A	K	0.00000	0.00000	0.0800000	0.00000	0.25000	0.25000	0.1167358
12	A	L	0.12500	0.12500	0.1650000	0.12500	0.25000	0.25000	0.0583679
14	A	N	0.00000	0.00000	0.0820000	0.00000	0.25000	0.25000	0.1174888
21	B	E	0.00000	0.00000	0.1255000	0.25000	0.25000	0.25000	0.1251242
22	B	F	0.12500	0.12500	0.1877500	0.25000	0.25000	0.25000	0.0625621
23	B	G	0.12500	0.12500	0.1877500	0.25000	0.25000	0.25000	0.0625621
25	B	I	0.00000	0.00000	0.0582500	0.00000	0.12500	0.12500	0.0624178
26	B	J	0.00000	0.00000	0.1165000	0.00000	0.25000	0.25000	0.1248356
30	B	N	0.00000	0.00000	0.1205000	0.00000	0.25000	0.25000	0.1250441
36	C	E	0.00000	0.00000	0.1060000	0.12500	0.12500	0.25000	0.0845772
37	C	F	0.12500	0.12500	0.1780000	0.18750	0.18750	0.25000	0.0422886
38	C	G	0.12500	0.12500	0.1780000	0.18750	0.18750	0.25000	0.0422886
40	C	I	0.12500	0.12500	0.1751250	0.18750	0.18750	0.25000	0.0406799
41	C	J	0.00000	0.00000	0.1002500	0.12500	0.12500	0.25000	0.0813599
42	C	K	0.00000	0.00000	0.0400000	0.00000	0.12500	0.12500	0.0583679
43	C	L	0.25000	0.25000	0.2700000	0.25000	0.31250	0.31250	0.0291840
45	C	N	0.00000	0.00000	0.1012500	0.12500	0.12500	0.25000	0.0875787
50	D	E	0.00000	0.00000	0.1060000	0.12500	0.12500	0.25000	0.0845772
51	D	F	0.25000	0.25000	0.3030000	0.31250	0.31250	0.37500	0.0422886
52	D	G	0.25000	0.25000	0.3030000	0.31250	0.31250	0.37500	0.0422886
54	D	I	0.12500	0.12500	0.1751250	0.18750	0.18750	0.25000	0.0406799
55	D	J	0.00000	0.00000	0.1002500	0.12500	0.12500	0.25000	0.0813599
56	D	K	0.00000	0.00000	0.0400000	0.00000	0.12500	0.12500	0.0583679
57	D	L	0.12500	0.12500	0.1450000	0.12500	0.18750	0.18750	0.0291840
59	D	N	0.00000	0.00000	0.1012500	0.12500	0.12500	0.25000	0.0875787
64	E	F	0.25000	0.25000	0.3030000	0.31250	0.31250	0.37500	0.0422886
65	E	G	0.25000	0.25000	0.3030000	0.31250	0.31250	0.37500	0.0422886
66	E	H	0.00000	0.00000	0.1060000	0.12500	0.12500	0.25000	0.0845772
67	E	I	0.00000	0.00000	0.0740000	0.06250	0.12500	0.25000	0.0751310
68	E	J	0.00000	0.00000	0.0615000	0.00000	0.12500	0.25000	0.0761041
69	E	K	0.00000	0.00000	0.0290000	0.00000	0.00000	0.12500	0.0528165
70	E	L	0.00000	0.00000	0.0675000	0.06250	0.12500	0.18750	0.0529976
71	E	M	0.00000	0.00000	0.0432500	0.00000	0.12500	0.12500	0.0595212
72	E	N	0.00000	0.00000	0.0562500	0.00000	0.12500	0.25000	0.0716062
74	E	P	0.00000	0.00000	0.0216250	0.00000	0.06250	0.06250	0.0297606
75	E	Q	0.00000	0.00000	0.0432500	0.00000	0.12500	0.12500	0.0595212
76	F	F	0.50000	0.50000	0.5530000	0.56250	0.56250	0.62500	0.0422886
77	F	G	0.25000	0.25000	0.3030000	0.31250	0.31250	0.37500	0.0422886
78	F	H	0.12500	0.12500	0.1780000	0.18750	0.18750	0.25000	0.0422886
79	F	I	0.06250	0.09375	0.1245625	0.12500	0.15625	0.25000	0.0473352
80	F	J	0.00000	0.00000	0.0808750	0.06250	0.12500	0.25000	0.0654838
81	F	K	0.00000	0.00000	0.0345000	0.00000	0.06250	0.12500	0.0443840
82	F	L	0.06250	0.09375	0.1062500	0.09375	0.12500	0.18750	0.0322061
83	F	M	0.06250	0.06250	0.0841250	0.06250	0.12500	0.12500	0.0297606
84	F	N	0.00000	0.00000	0.0787500	0.06250	0.12500	0.25000	0.0675096
86	F	P	0.03125	0.03125	0.0420625	0.03125	0.06250	0.06250	0.0148803
87	F	Q	0.06250	0.06250	0.0841250	0.06250	0.12500	0.12500	0.0297606
88	G	G	0.50000	0.50000	0.5530000	0.56250	0.56250	0.62500	0.0422886
89	G	H	0.12500	0.12500	0.1780000	0.18750	0.18750	0.25000	0.0422886
90	G	I	0.06250	0.09375	0.1245625	0.12500	0.15625	0.25000	0.0473352
91	G	J	0.00000	0.00000	0.0808750	0.06250	0.12500	0.25000	0.0654838
92	G	K	0.00000	0.00000	0.0345000	0.00000	0.06250	0.12500	0.0443840
93	G	L	0.06250	0.09375	0.1062500	0.09375	0.12500	0.18750	0.0322061
94	G	M	0.06250	0.06250	0.0841250	0.06250	0.12500	0.12500	0.0297606
95	G	N	0.00000	0.00000	0.0787500	0.06250	0.12500	0.25000	0.0675096
97	G	P	0.03125	0.03125	0.0420625	0.03125	0.06250	0.06250	0.0148803
98	G	Q	0.06250	0.06250	0.0841250	0.06250	0.12500	0.12500	0.0297606
100	H	I	0.12500	0.12500	0.1751250	0.18750	0.18750	0.25000	0.0406799
101	H	J	0.00000	0.00000	0.1002500	0.12500	0.12500	0.25000	0.0813599
102	H	K	0.00000	0.00000	0.0400000	0.00000	0.12500	0.12500	0.0583679
103	H	L	0.12500	0.12500	0.1450000	0.12500	0.18750	0.18750	0.0291840
105	H	N	0.00000	0.00000	0.1012500	0.12500	0.12500	0.25000	0.0875787
109	I	I	0.50000	0.50000	0.5420000	0.50000	0.62500	0.62500	0.0591015
110	I	J	0.25000	0.25000	0.2920000	0.25000	0.37500	0.37500	0.0591015
111	I	K	0.00000	0.00000	0.0523750	0.00000	0.12500	0.18750	0.0671421
112	I	L	0.06250	0.09375	0.1137500	0.09375	0.12500	0.21875	0.0393995
113	I	M	0.12500	0.12500	0.1460000	0.12500	0.18750	0.18750	0.0295507
114	I	N	0.00000	0.00000	0.0670000	0.06250	0.12500	0.25000	0.0734629
116	I	P	0.06250	0.06250	0.0730000	0.06250	0.09375	0.09375	0.0147754
117	I	Q	0.12500	0.12500	0.1460000	0.12500	0.18750	0.18750	0.0295507
119	J	K	0.00000	0.00000	0.0247500	0.00000	0.00000	0.12500	0.0498614
120	J	L	0.00000	0.00000	0.0625000	0.06250	0.09375	0.18750	0.0498934
121	J	M	0.00000	0.00000	0.0420000	0.00000	0.12500	0.12500	0.0591015
122	J	N	0.00000	0.00000	0.0520000	0.00000	0.12500	0.25000	0.0728465
124	J	P	0.00000	0.00000	0.0210000	0.00000	0.06250	0.06250	0.0295507
125	J	Q	0.00000	0.00000	0.0420000	0.00000	0.12500	0.12500	0.0591015
127	K	L	0.25000	0.25000	0.2700000	0.25000	0.31250	0.31250	0.0291840
128	K	M	0.00000	0.00000	0.0400000	0.00000	0.12500	0.12500	0.0583679
129	K	N	0.00000	0.00000	0.0275000	0.00000	0.00000	0.12500	0.0518326
131	K	P	0.00000	0.00000	0.0200000	0.00000	0.06250	0.06250	0.0291840
132	K	Q	0.00000	0.00000	0.0400000	0.00000	0.12500	0.12500	0.0583679
133	L	L	0.50000	0.50000	0.5200000	0.50000	0.56250	0.56250	0.0291840
134	L	M	0.06250	0.06250	0.0825000	0.06250	0.12500	0.12500	0.0291840
135	L	N	0.00000	0.00000	0.0643750	0.06250	0.12500	0.18750	0.0540035
137	L	P	0.03125	0.03125	0.0412500	0.03125	0.06250	0.06250	0.0145920
138	L	Q	0.06250	0.06250	0.0825000	0.06250	0.12500	0.12500	0.0291840
140	M	N	0.00000	0.00000	0.0410000	0.00000	0.12500	0.12500	0.0587444
146	N	P	0.00000	0.00000	0.0205000	0.00000	0.06250	0.06250	0.0293722
147	N	Q	0.00000	0.00000	0.0410000	0.00000	0.12500	0.12500	0.0587444

Comparing the values and variation found for the various kinship values:

stats_short <- stats_10[stats_10$sd > 0.0, ]
stats_long <- stats_500[stats_500$sd > 0.0, ]

## Align the short and long summaries by animal pair: the sd > 0 sets can
## differ in size between simulation lengths, so join on the ID pair rather
## than by row position.
merged <- merge(stats_short, stats_long, by = c("id_1", "id_2"),
                suffixes = c("_short", "_long"))

comprison <- data.frame(
  id_1 = merged$id_1,
  id_2 = merged$id_2,
  meanKin_short = merged$mean_short,
  meanKin_long = merged$mean_long,
  meanKinDelta = abs(merged$mean_short - merged$mean_long),
  sdKin_short = merged$sd_short,
  sdKin_long = merged$sd_long,
  sdKinDelta = abs(merged$sd_short - merged$sd_long)
)

kable(comprison,
  longtable = TRUE,
  digits = c(0L, 0L, 4L, 4L, 4L, 4L, 4L, 4L),
  caption = stri_c(
    "Comparision of estimated kinships between simulations ",
    "of 10 (short) and 500 (long)"
  )
) |>
  kable_styling(
    latex_options = c("striped", "repeat_header"),
    repeat_header_method = "replace",
    repeat_header_text = "\\textit{(continued)}",
    font_size = 10L
  )

Comparision of estimated kinships between simulations of 10 (short) and 500 (long)
id_1	id_2	meanKin_short	meanKin_long	meanKinDelta	sdKin_short	sdKin_long	sdKinDelta
A	E	0.0750	0.0865	0.0115	0.1208	0.1190	0.0017
A	F	0.1625	0.1683	0.0058	0.0604	0.0595	0.0009
A	G	0.1625	0.1683	0.0058	0.0604	0.0595	0.0009
A	I	0.3125	0.2920	0.0205	0.0659	0.0591	0.0068
A	J	0.1250	0.0840	0.0410	0.1318	0.1182	0.0136
A	N	0.0500	0.0820	0.0320	0.1054	0.1175	0.0121
B	E	0.1500	0.1255	0.0245	0.1291	0.1251	0.0040
B	F	0.2000	0.1878	0.0123	0.0645	0.0626	0.0020
B	G	0.2000	0.1878	0.0123	0.0645	0.0626	0.0020
B	I	0.0750	0.0583	0.0167	0.0645	0.0624	0.0021
B	J	0.1500	0.1165	0.0335	0.1291	0.1248	0.0043
B	N	0.0500	0.1205	0.0705	0.1054	0.1250	0.0196
C	E	0.1125	0.1060	0.0065	0.0922	0.0846	0.0077
C	F	0.1812	0.1780	0.0033	0.0461	0.0423	0.0038
C	G	0.1812	0.1780	0.0033	0.0461	0.0423	0.0038
C	I	0.1938	0.1751	0.0186	0.0355	0.0407	0.0052
C	J	0.1375	0.1003	0.0373	0.0710	0.0814	0.0104
C	N	0.0500	0.1013	0.0513	0.0874	0.0876	0.0002
D	E	0.1125	0.1060	0.0065	0.0922	0.0846	0.0077
D	F	0.3062	0.3030	0.0033	0.0461	0.0423	0.0038
D	G	0.3062	0.3030	0.0033	0.0461	0.0423	0.0038
D	I	0.1938	0.1751	0.0186	0.0355	0.0407	0.0052
D	J	0.1375	0.1003	0.0373	0.0710	0.0814	0.0104
D	N	0.0500	0.1013	0.0513	0.0874	0.0876	0.0002
E	F	0.3062	0.3030	0.0033	0.0461	0.0423	0.0038
E	G	0.3062	0.3030	0.0033	0.0461	0.0423	0.0038
E	H	0.1125	0.1060	0.0065	0.0922	0.0846	0.0077
E	I	0.0875	0.0740	0.0135	0.0791	0.0751	0.0039
E	J	0.1000	0.0615	0.0385	0.0791	0.0761	0.0030
E	K	0.0625	0.0290	0.0335	0.0659	0.0528	0.0131
E	L	0.0875	0.0675	0.0200	0.0323	0.0530	0.0207
E	M	0.0375	0.0432	0.0057	0.0604	0.0595	0.0009
E	N	0.0375	0.0562	0.0188	0.0604	0.0716	0.0112
E	P	0.0187	0.0216	0.0029	0.0302	0.0298	0.0004
E	Q	0.0375	0.0432	0.0057	0.0604	0.0595	0.0009
F	F	0.5562	0.5530	0.0032	0.0461	0.0423	0.0038
F	G	0.3062	0.3030	0.0033	0.0461	0.0423	0.0038
F	H	0.1812	0.1780	0.0033	0.0461	0.0423	0.0038
F	I	0.1406	0.1246	0.0161	0.0423	0.0473	0.0050
F	J	0.1187	0.0809	0.0379	0.0547	0.0655	0.0108
F	K	0.0312	0.0345	0.0033	0.0329	0.0444	0.0114
F	L	0.1062	0.1062	0.0000	0.0161	0.0322	0.0161
F	M	0.0813	0.0841	0.0029	0.0302	0.0298	0.0004
F	N	0.0437	0.0788	0.0350	0.0593	0.0675	0.0082
F	P	0.0406	0.0421	0.0014	0.0151	0.0149	0.0002
F	Q	0.0813	0.0841	0.0029	0.0302	0.0298	0.0004
G	G	0.5562	0.5530	0.0032	0.0461	0.0423	0.0038
G	H	0.1812	0.1780	0.0033	0.0461	0.0423	0.0038
G	I	0.1406	0.1246	0.0161	0.0423	0.0473	0.0050
G	J	0.1187	0.0809	0.0379	0.0547	0.0655	0.0108
G	K	0.0312	0.0345	0.0033	0.0329	0.0444	0.0114
G	L	0.1062	0.1062	0.0000	0.0161	0.0322	0.0161
G	M	0.0813	0.0841	0.0029	0.0302	0.0298	0.0004
G	N	0.0437	0.0788	0.0350	0.0593	0.0675	0.0082
G	P	0.0406	0.0421	0.0014	0.0151	0.0149	0.0002
G	Q	0.0813	0.0841	0.0029	0.0302	0.0298	0.0004
H	I	0.1938	0.1751	0.0186	0.0355	0.0407	0.0052
H	J	0.1375	0.1003	0.0373	0.0710	0.0814	0.0104
H	N	0.0500	0.1013	0.0513	0.0874	0.0876	0.0002
I	I	0.5625	0.5420	0.0205	0.0659	0.0591	0.0068
I	J	0.3125	0.2920	0.0205	0.0659	0.0591	0.0068
I	K	0.0063	0.0524	0.0461	0.0198	0.0671	0.0474
I	L	0.1000	0.1138	0.0137	0.0198	0.0394	0.0196
I	M	0.1562	0.1460	0.0103	0.0329	0.0296	0.0034
I	N	0.0562	0.0670	0.0108	0.0952	0.0735	0.0218
I	P	0.0781	0.0730	0.0051	0.0165	0.0148	0.0017
I	Q	0.1562	0.1460	0.0103	0.0329	0.0296	0.0034
J	K	0.0125	0.0248	0.0122	0.0395	0.0499	0.0103
J	L	0.0750	0.0625	0.0125	0.0395	0.0499	0.0104
J	M	0.0625	0.0420	0.0205	0.0659	0.0591	0.0068
J	N	0.0625	0.0520	0.0105	0.1062	0.0728	0.0334
J	P	0.0312	0.0210	0.0102	0.0329	0.0296	0.0034
J	Q	0.0625	0.0420	0.0205	0.0659	0.0591	0.0068
K	N	0.0250	0.0275	0.0025	0.0527	0.0518	0.0009
L	N	0.0375	0.0644	0.0269	0.0527	0.0540	0.0013
M	N	0.0250	0.0410	0.0160	0.0527	0.0587	0.0060
N	P	0.0125	0.0205	0.0080	0.0264	0.0294	0.0030
N	Q	0.0250	0.0410	0.0160	0.0527	0.0587	0.0060

stats_short <- stats_100[stats_100$sd > 0.0, ]
stats_long <- stats_500[stats_500$sd > 0.0, ]

## Align the short and long summaries by animal pair: the sd > 0 sets can
## differ in size between simulation lengths, so join on the ID pair rather
## than by row position.
merged <- merge(stats_short, stats_long, by = c("id_1", "id_2"),
                suffixes = c("_short", "_long"))

comprison <- data.frame(
  id_1 = merged$id_1,
  id_2 = merged$id_2,
  meanKin_short = merged$mean_short,
  meanKin_long = merged$mean_long,
  meanKinDelta = abs(merged$mean_short - merged$mean_long),
  sdKin_short = merged$sd_short,
  sdKin_long = merged$sd_long,
  sdKinDelta = abs(merged$sd_short - merged$sd_long)
)

kable(comprison,
  longtable = TRUE,
  digits = c(0L, 0L, 4L, 4L, 4L, 4L, 4L, 4L),
  caption = stri_c(
    "Comparision of estimated kinships between simulations ",
    "of 100 (short) and 500 (long)"
  )
) |>
  kable_styling(
    latex_options = c("striped", "repeat_header"),
    repeat_header_method = "replace",
    repeat_header_text = "\\textit{(continued)}",
    font_size = 10L
  )

Comparision of estimated kinships between simulations of 100 (short) and 500 (long)
id_1	id_2	meanKin_short	meanKin_long	meanKinDelta	sdKin_short	sdKin_long	sdKinDelta
A	E	0.0825	0.0865	0.0040	0.1181	0.1190	0.0009
A	F	0.1663	0.1683	0.0020	0.0591	0.0595	0.0004
A	G	0.1663	0.1683	0.0020	0.0591	0.0595	0.0004
A	I	0.2988	0.2920	0.0068	0.0613	0.0591	0.0022
A	J	0.0975	0.0840	0.0135	0.1226	0.1182	0.0043
A	K	0.0700	0.0800	0.0100	0.1128	0.1167	0.0039
A	L	0.1600	0.1650	0.0050	0.0564	0.0584	0.0020
A	N	0.0800	0.0820	0.0020	0.1172	0.1175	0.0003
B	E	0.1325	0.1255	0.0070	0.1254	0.1251	0.0003
B	F	0.1912	0.1878	0.0035	0.0627	0.0626	0.0001
B	G	0.1912	0.1878	0.0035	0.0627	0.0626	0.0001
B	I	0.0625	0.0583	0.0042	0.0628	0.0624	0.0004
B	J	0.1250	0.1165	0.0085	0.1256	0.1248	0.0008
B	N	0.1100	0.1205	0.0105	0.1247	0.1250	0.0003
C	E	0.1075	0.1060	0.0015	0.0871	0.0846	0.0025
C	F	0.1788	0.1780	0.0008	0.0435	0.0423	0.0012
C	G	0.1788	0.1780	0.0008	0.0435	0.0423	0.0012
C	I	0.1806	0.1751	0.0055	0.0386	0.0407	0.0021
C	J	0.1113	0.1003	0.0110	0.0772	0.0814	0.0041
C	K	0.0350	0.0400	0.0050	0.0564	0.0584	0.0020
C	L	0.2675	0.2700	0.0025	0.0282	0.0292	0.0010
C	N	0.0950	0.1013	0.0063	0.0925	0.0876	0.0049
D	E	0.1075	0.1060	0.0015	0.0871	0.0846	0.0025
D	F	0.3038	0.3030	0.0008	0.0435	0.0423	0.0012
D	G	0.3038	0.3030	0.0008	0.0435	0.0423	0.0012
D	I	0.1806	0.1751	0.0055	0.0386	0.0407	0.0021
D	J	0.1113	0.1003	0.0110	0.0772	0.0814	0.0041
D	K	0.0350	0.0400	0.0050	0.0564	0.0584	0.0020
D	L	0.1425	0.1450	0.0025	0.0282	0.0292	0.0010
D	N	0.0950	0.1013	0.0063	0.0925	0.0876	0.0049
E	F	0.3038	0.3030	0.0008	0.0435	0.0423	0.0012
E	G	0.3038	0.3030	0.0008	0.0435	0.0423	0.0012
E	H	0.1075	0.1060	0.0015	0.0871	0.0846	0.0025
E	I	0.0762	0.0740	0.0023	0.0730	0.0751	0.0021
E	J	0.0700	0.0615	0.0085	0.0801	0.0761	0.0040
E	K	0.0388	0.0290	0.0097	0.0581	0.0528	0.0053
E	L	0.0731	0.0675	0.0056	0.0533	0.0530	0.0003
E	M	0.0412	0.0432	0.0020	0.0591	0.0595	0.0004
E	N	0.0625	0.0562	0.0062	0.0743	0.0716	0.0027
E	P	0.0206	0.0216	0.0010	0.0295	0.0298	0.0002
E	Q	0.0412	0.0432	0.0020	0.0591	0.0595	0.0004
F	F	0.5538	0.5530	0.0007	0.0435	0.0423	0.0012
F	G	0.3038	0.3030	0.0008	0.0435	0.0423	0.0012
F	H	0.1788	0.1780	0.0008	0.0435	0.0423	0.0012
F	I	0.1284	0.1246	0.0039	0.0426	0.0473	0.0048
F	J	0.0906	0.0809	0.0097	0.0637	0.0655	0.0018
F	K	0.0369	0.0345	0.0024	0.0436	0.0444	0.0008
F	L	0.1078	0.1062	0.0016	0.0321	0.0322	0.0001
F	M	0.0831	0.0841	0.0010	0.0295	0.0298	0.0002
F	N	0.0788	0.0788	0.0000	0.0720	0.0675	0.0045
F	P	0.0416	0.0421	0.0005	0.0148	0.0149	0.0001
F	Q	0.0831	0.0841	0.0010	0.0295	0.0298	0.0002
G	G	0.5538	0.5530	0.0007	0.0435	0.0423	0.0012
G	H	0.1788	0.1780	0.0008	0.0435	0.0423	0.0012
G	I	0.1284	0.1246	0.0039	0.0426	0.0473	0.0048
G	J	0.0906	0.0809	0.0097	0.0637	0.0655	0.0018
G	K	0.0369	0.0345	0.0024	0.0436	0.0444	0.0008
G	L	0.1078	0.1062	0.0016	0.0321	0.0322	0.0001
G	M	0.0831	0.0841	0.0010	0.0295	0.0298	0.0002
G	N	0.0788	0.0788	0.0000	0.0720	0.0675	0.0045
G	P	0.0416	0.0421	0.0005	0.0148	0.0149	0.0001
G	Q	0.0831	0.0841	0.0010	0.0295	0.0298	0.0002
H	I	0.1806	0.1751	0.0055	0.0386	0.0407	0.0021
H	J	0.1113	0.1003	0.0110	0.0772	0.0814	0.0041
H	K	0.0350	0.0400	0.0050	0.0564	0.0584	0.0020
H	L	0.1425	0.1450	0.0025	0.0282	0.0292	0.0010
H	N	0.0950	0.1013	0.0063	0.0925	0.0876	0.0049
I	I	0.5488	0.5420	0.0067	0.0613	0.0591	0.0022
I	J	0.2988	0.2920	0.0068	0.0613	0.0591	0.0022
I	K	0.0469	0.0524	0.0055	0.0649	0.0671	0.0022
I	L	0.1138	0.1138	0.0000	0.0360	0.0394	0.0034
I	M	0.1494	0.1460	0.0034	0.0306	0.0296	0.0011
I	N	0.0638	0.0670	0.0033	0.0794	0.0735	0.0060
I	P	0.0747	0.0730	0.0017	0.0153	0.0148	0.0005
I	Q	0.1494	0.1460	0.0034	0.0306	0.0296	0.0011
J	K	0.0238	0.0248	0.0010	0.0493	0.0499	0.0006
J	L	0.0675	0.0625	0.0050	0.0467	0.0499	0.0032
J	M	0.0488	0.0420	0.0067	0.0613	0.0591	0.0022
J	N	0.0475	0.0520	0.0045	0.0728	0.0728	0.0001
J	P	0.0244	0.0210	0.0034	0.0306	0.0296	0.0011
J	Q	0.0488	0.0420	0.0067	0.0613	0.0591	0.0022
K	L	0.2675	0.2700	0.0025	0.0282	0.0292	0.0010
K	M	0.0350	0.0400	0.0050	0.0564	0.0584	0.0020
K	N	0.0288	0.0275	0.0013	0.0529	0.0518	0.0010
K	P	0.0175	0.0200	0.0025	0.0282	0.0292	0.0010
K	Q	0.0350	0.0400	0.0050	0.0564	0.0584	0.0020
L	L	0.5175	0.5200	0.0025	0.0282	0.0292	0.0010
L	M	0.0800	0.0825	0.0025	0.0282	0.0292	0.0010
L	N	0.0619	0.0644	0.0025	0.0537	0.0540	0.0003
L	P	0.0400	0.0412	0.0013	0.0141	0.0146	0.0005
L	Q	0.0800	0.0825	0.0025	0.0282	0.0292	0.0010
M	N	0.0400	0.0410	0.0010	0.0586	0.0587	0.0001
N	P	0.0200	0.0205	0.0005	0.0293	0.0294	0.0001
N	Q	0.0400	0.0410	0.0010	0.0586	0.0587	0.0001

nprcgenekeepr: an R Package for the Genetic Management of Colonies

R. Mark Sharp, Ph.D.

8/29/2021

Introduction

Anticipation of further uses of this kinship estimation method

Identification of potential parents

Creation of example pedigree²

Small Example

Quick Example of Simulated Kinships with Partial Parentage

nprcgenekeepr: an R Package for the Genetic Management of Colonies

R. Mark Sharp, Ph.D.

8/29/2021

Introduction

Anticipation of further uses of this kinship estimation method

Identification of potential parents

Creation of example pedigree2

Small Example

Creation of example pedigree²