Thanks so much for your help ../__\ :::: It is a sensitive data set so I cannot share the data set with you. It is so kind of you to offer my guidance and support. Perhaps I can leave my code here and see if I have done anything from when shortening the scale:
I used the following syntax for testing the bifactor model of the Social Potency (8 facets, 1 factor)
dataset = read.spss(file.choose(), to.data.frame=TRUE)
specific <-c(1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,4,4,4,5,5,5,5,5,5,5,5,5,5,6,6,6,6,6,6,6,6,6,6,7,7,7,7,7,7,7,7,7,7,8,8,8,8,8,8,8,8,8,8)
b_mod1<-bfactor(dataset,specific, SE=TRUE)
coef(b_mod1, printSE=TRUE)
summary(b_mod1)
residuals(b_mod1)
itemfit(b_mod1, QMC=TRUE)
M2(b_mod1, QMC=TRUE)
I selected items with slope values of .9 or higher, paired items with LD values (below 10) . I also checked the model fit.