I have a six-item psychometric scale each item having 6 answer categories. In the pilot testing I used EFA and CFA on the Pearson matrix on a sample n=511 and they supported a bifactorial structure with F1(3 items) F2(3 items) good loadings, no cross-loadings, good fit indices, factors correlated r=0,4. Now, in the validation phase (sample smaller n=200) I decided to us polychoric matrix (ordinal level data) and varied factor diagnostic (BIC, PA, Hull method) and extraction methods (ULS, ML, MWLS, Minimum Rank Factor Analysis) and many combination of them in FACTOR software. Whether I configured the test for 1 factor or two, hals of these method combination suggested 1 dimension and about half 2. Typically when I used PA it suggested 1. But if I used ULS and DWLS they suggested 2 when used with any of the three factor retention methods. Explained variance was about 50% in the 1 factor models and about 70% in the 2 factor models. The CFA done in R lavaan package supported two factors. Finally I did Mokken scaling. It really turned everything upside down. It seems to support a clear unidimensional scale with coef H=0,47. I checked this for the two possible subscales too an they also showed monotonicity with coefH 0,5 and 0,7 respectively. Now I have no idea whether my scale is on dimensional or two. Can you help please?