tests/failure.R

In kinship2: Pedigree Functions

require(kinship2)
#
# Here is a case where the levels fail to line up properly
#

#data(testped2)
data(sample.ped)

# rearrange the founders to get a nicer plot
df1<- sample.ped[sample.ped$ped==1,]

ped1 <- with(df1, pedigree(id, father, mother, sex, affected))

#plot(ped1)

df1reord <- df1[c(35:41,1:34),]
ped1reord <- with(df1reord, pedigree(id, father, mother, 
       sex, affected=affected))


#plot(ped1reord, col=df1reord$avail+1)


# Two brothers married two sisters, which is currently "too much" for
#  the kindepth routine.