Perform the following steps in R:

>install.packages("TraMineR")
> library(TraMineR)
> data("mvad")
> head(mvad)
>alphabet = seqstatl(mvad[,17:86])
Output
    [1] "employment" "FE" "HE" "joblessness" "school" 
    "training" 
    > fulllabel<- c("employment", "further education",
    "higher education","joblessness", "school", "training")
    > shortlabel<- c("EM", "FE", "HE", "JL", "SC", "TR")
    > seq_mvad<- seqdef(mvad[, 17:86], alphabet = alphabet,
    states = shortlabel, labels = fulllabel, weights = mvad$weight,
    xtstep = 6)
    Output
    [>] state coding:
     [alphabet] [label] [long label] 
    1 employment EM employment
    2 FEFE further education
    3 HEHE higher education
    4 joblessness JL joblessness
    5 school SC school
    6 training TR training
     [>] sum of weights: 711.57 - min/max: 0.13/4.46
     [>] 712 sequences in the data set
     [>] min/max sequence length: 70/70
    > seq_mvad[1:2,]
    Output
    Sequence 
    1 EM-EM-EM-EM-TR-TR-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM
    -EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-
    EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM
    -EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM-EM
    2 FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-
    FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-FE-
    FE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-
    HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE-HE
    > print(seq_mvad[1:2,],format="SPS")
    Output
     Sequence 
    [1] (EM,4)-(TR,2)-(EM,64)
    [2] (FE,36)-(HE,34)