ped.param is a list of 1 element:
alpha a list of d=4 (the number of measurements) entries, each is a matrix of K=3 (the number of classes) times S[j]
(the number of possible values of measurement j), a row alpha[[j]][k,] contains the logistic coefficients of the measurement j for class k.