R Under development (unstable) (2025-09-02 r88773 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> #######################################
> ## Name: test.pedigree.unrelated.r
> ## Purpose: Test Suite for pedigree.unrelated
> ## Created: 3/29/2011
> ## Last Updated: 7/13/2011
> ## Author: Jason Sinnwell, MS
> ########################################
> 
> ## examples from help file, available with
> ## > example(pedigree.unrelated)
> 
> require(kinship2)
Loading required package: kinship2
Loading required package: Matrix
Loading required package: quadprog
Warning: kinship2 package is deprecated for R <= 4.5; switch functionality to Pedixplorer from BioConductor> #library(kinship2, lib.loc="~/Rdir/library")
> #citation("kinship2")
> 
> data(sample.ped)
> 
> 
> pedAll <- pedigree(sample.ped$id, sample.ped$father, sample.ped$mother,
+                    sample.ped$sex, famid=sample.ped$ped,
+        affected=cbind(sample.ped$affected, sample.ped$avail))
>                    
> ped1 <- pedAll['1']
> ped2 <- pedAll['2']
> 
> ## to see plot:
> ## plot.pedigree(ped1, align=FALSE)
> set.seed(10)
> id1 <- pedigree.unrelated(ped1, avail=ped1$affected[,2])
> 
> ## some possible vectors
> id1
[1] "109" "113" "133" "141"
> # "109" "113" "133"
> # "109" "110" "130"
> # "109" "118" "141"
> 
> set.seed(10)
> id2 <- pedigree.unrelated(ped2, avail=ped2$affected[,2])
> 
> ## some possible vectors
> id2
[1] "203" "206"
> ##[1] "203" "206"
> ##[1] "203" "213"
> ##[1] "203" "204"
> 
> 
> proc.time()
   user  system elapsed 
   0.93    0.14    1.06