# NOT RUN {
data(biofam, package="MEDseq")
biofam <- list(covariates = biofam[2L:9L], sequences = biofam[10L:25L] + 1L)
biofam.cov <- biofam$covariates[,colSums(is.na(biofam$covariates)) == 0]
biofam.seq <- seqdef(biofam$sequences,
states = c("P", "L", "M", "L+M", "C", "L+C", "L+M+C", "D"),
labels = c("Parent", "Left", "Married", "Left+Marr", "Child",
"Left+Child", "Left+Marr+Child", "Divorced"))
# }
Run the code above in your browser using DataLab