Sure! Here is my code:
##Group A Data
M_GrpA = read.csv("M_GRPA.csv", header = T)
M_GrpA_Items = M_GrpA[ , c(6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26)]
M_GrpA_Items = na.omit(M_GrpA_Items)
##Group B Data
M_GrpB = read.csv("M_GRPB.csv", header = T)
M_GrpB_Items = M_GrpB[ , c(6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26)]
M_GrpB_Items = na.omit(M_GrpB_Items)
##Binding retrieved data into single file
M_GrpAvB_items = rbind(M_GrpA_Items, M_GrpB_Items)
groups = c(rep('A',1820), rep('B', 481))
model ='TMI = 1-21'
#Configural invariance model
AvB_config = multipleGroup(M_GrpAvB_items, model, group = groups,
technical = list(NCYCLES=5000))
itemfit(AvB_config)
Thank you for your help!