R/kin2gds.R

In UW-GAC/GENESIS: GENetic EStimation and Inference in Structured samples (GENESIS): Statistical methods for analyzing genetic data from samples with population structure and/or relatedness

mat2gds <- function(mat, gdsfile) {
    ## if we have a Matrix object, coerce
    if (!is.matrix(mat)) mat <- as.matrix(mat)
    gds <- createfn.gds(gdsfile)
    add.gdsn(gds, "sample.id", colnames(mat))
    add.gdsn(gds, "kinship", mat)
    closefn.gds(gds)
}


kin2gds <- function(ibdobj, gdsfile) {
    gds <- createfn.gds(gdsfile)
    add.gdsn(gds, "sample.id", ibdobj$sample.id)
    add.gdsn(gds, "kinship", ibdobj$kinship)
    closefn.gds(gds)
}