Skip to contents

Test if a degree of relatedness is consistent with an observed PMR

Usage

test_degree(in_tibble, row, degree, verbose = TRUE)

Arguments

in_tibble

a tibble that is the output of the callRelatedness() function.

row

either the row number or pair name for which the posterior distribution is to be plotted.

degree

the degree of relatedness to be tested.

verbose

a logical (boolean) for whether all test output should be printed to screen.

Value

the associated p-value for the test

Examples

test_degree(relatedness_example, 1, 1)
#> Testing H0       : "Ind1 - Ind2" are 1st-degree relatives.
#> Expected PMR     : 0.1635
#> Observed PMR     : 0.2042
#> Estimated degree : 2.9826
#> p-value          : 2.666909e-05
#> Decision         : Reject H0
#> [1] 2.666909e-05