test_that("ped -> matrix -> ped", {
x = ped(id=1:3, fid=c(0,0,1), mid=c(0,0,2), sex=c(1,2,1))
y = restorePed(as.matrix(x))
expect_identical(x, y)
x2 = reorderPed(x, 3:1)
y2 = restorePed(as.matrix(x2))
expect_identical(x2, y2)
})
test_that("as.ped.data.frame() keeps genotypes", {
# Without FAMID
df = data.frame(id=1:2, fid=0, mid=0, sex=1, A1=1, A2=2)
p = as.ped(df)
expect_identical(genotype(p[[1]],id = 1, markers = 1), c("1", "2"))
expect_identical(genotype(p[[2]],id = 2, markers = 1), c("1", "2"))
# With FAMID, same ID
df2 = data.frame(famid=1:2, id=1, fid=0, mid=0, sex=1, A1=1, A2=2)
pp = as.ped(df2)
expect_identical(genotype(pp[[1]],id = 1, markers = 1), c("1", "2"))
expect_identical(genotype(pp[[2]],id = 1, markers = 1), c("1", "2"))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.