rscharpf/MinimumDistance — Issue #14: trioSet Subsetting Problem
Author: syounkin
Created: 2012-02-01 21:18:11+00:00
State: open

When we subset the trioSet object we sometimes get NAs for parental IDs that weren't there before subsetting. See code below.

> sum(is.na(MinimumDistance:::trioIndex(pedigree(pitt.trioSet))$individualId))
[1] 0
> sum(is.na(MinimumDistance:::trioIndex(pedigree(pitt.trioSet[,303,]))$individualId))
[1] 2

See enigma2:~syounkin/md-devel/ for generating script pitt-md-2.R and binary pitt-md.RData.

I think it has to do with duplicate parental IDs.

I'll keep looking into it and post if I make any progress.