disease=c(rep(1,105),rep(0,0),rep(1,1),rep(0,0),rep(1,4),rep(0,0),rep(1,4),rep(0,3),rep(1,10),rep(0,1),rep(1,3),rep(0,3),rep(1,62),rep(0,16),rep(1,444),rep(0,724)) tired=c(rep(1,105),rep(0,1),rep(1,4),rep(0,7),rep(1,11),rep(0,6),rep(1,78),rep(0,1168)) dishonoured=c(rep(1,106),rep(0,11),rep(1,17),rep(0,1246)) bankrupt=c(rep(1,117),rep(0,1263)) suicide=cbind(disease,tired,dishonoured,bankrupt) lca(suicide,2) lca(suicide,3) lca(suicide,4)