Dear members,

I am using lclogit for fitting latent class conditional logit models with the guide provided by Pacifico and Yoo.

To predict the unconditional probability of actual choice and the probability of actual choice conditional on being in class c:

I use the following as instructed.
Code.
Code:
lclogitpr pr, pr
generate byte class = .
forvalues c = 1/`e(nclasses)' {    
    quietly replace class = `c' if cpmax==cp`c'    
}    
forvalues c = 1/`e(nclasses)' {    
    quietly summarize pr if class == `c' & choi ==1
    local n=r(N)
    local a=r(mean)
    quietly summarize pr`c' if class == `c' & choi ==1
    local b=r(mean)
    matrix pr = nullmat(pr) \ `n', `c', `a', `b'    
}
matrix colnames pr = "Obs" "Class" "Uncond_Pr" "Cond_PR"
matlist pr, name(columns)
It works fine when I use model with 5 classes but when I want 3 classes what do I do? because at the moment I get output like this which makes no sense to me.
Obs Class Uncond_Pr Cond_PR
-------------------------------------------
336 1 .400312 .6827738
300 2 .3449114 .6237895
192 3 .4260958 .8268863
336 4 .3319297 .9369147
510 5 .4137588 .6265153
252 1 .3815055 .9111855
216 2 .4584624 .7100105
696 3 .3878322 .4982426
288 4 .4446026 .614101
204 1 .3674326 .5913422
204 2 .4262477 .7003562
798 3 .3721382 .3802904
432 4 .4378642 .6870132
150 1 .3719234 .8538355
330 2 .3474019 .7184852
534 3 .4206686 .4666928
390 4 .4663766 .7157002
234 5 .4264594 .7264037
150 1 .3719234 .8538355
330 2 .3474019 .7184852
534 3 .4206686 .4666928
390 4 .4663766 .7157002
234 5 .4264594 .7264037
150 1 .3719231 .8538312
330 2 .3474018 .7184852
534 3 .4206686 .4666928
390 4 .4663763 .7157001
234 5 .4264591 .7264037
150 1 .3719231 .8538312
330 2 .3474018 .7184852
534 3 .4206686 .4666928
390 4 .4663763 .7157001
234 5 .4264591 .7264037
150 1 .3719231 .8538312
330 2 .3474018 .7184852
534 3 .4206686 .4666928
390 4 .4663763 .7157001
234 5 .4264591 .7264037
252 1 .3815055 .9111855
216 2 .4584624 .7100105
696 3 .3878322 .4982426
288 4 .4446026 .614101
252 1 .3815055 .9111855
216 2 .4584624 .7100105
696 3 .3878322 .4982426
288 4 .4446026 .614101
534 1 .4168518 .705491
396 2 .3826234 .6238397
486 3 .339419 .460339


I would have expected output for 3 classes only.
Thank you