## Defining a state sequence object
data(mvad)
mvad.seq <- seqdef(mvad[, 17:86])
## Building dissimilarities
mvad.lcs <- seqdist(mvad.seq, method="LCS")
## Compute distance to center according to group gcse5eq
dc <- disscenter(mvad.lcs,group=mvad$gcse5eq)
## Ploting distribution of dissimilarity to center
boxplot(dc~mvad$gcse5eq, col="cyan")
Run the code above in your browser using DataLab