KP_2PVR_obtain_extAssemblies_and_extPairs.m

function [extraAssemblies,extraPairs,extAssemblies,matExtraAssemblies,matExtAssemblies] = KP_2PVR_obtain_extAssemblies_and_extPairs(assemblies,pairsInfo)
extraAssemblies = cell(size(assemblies));
extraPairs = cell(size(assemblies));
extAssemblies = assemblies;
for a = 1:length(assemblies)
    extraAssemblies{a} = [];
    extraPairs{a} = [];
    assembly = assemblies{a};
    sizeAssembly = length(assembly);
    for n = 1:sizeAssembly
        neuron = assembly(n);
        otherNeurons = assembly(~ismember(assembly,neuron));
        for i = 1:size(pairsInfo,1)
            neuronPair = pairsInfo(i,[7,8]);
            resAnchor = find(neuronPair == neuron);
            if ~isempty(resAnchor)
                if resAnchor == 1
                    queryNeuron = neuronPair(2);
                elseif resAnchor == 2
                    queryNeuron = neuronPair(1);
                end
                if ~ismember(queryNeuron,otherNeurons)
                    sectorA = pairsInfo(i,5);
                    sectorB = pairsInfo(i,6);
                    if sectorA == sectorB
                        extraAssemblies{a} = [extraAssemblies{a},queryNeuron];
                        extraPairs{a} = [extraPairs{a};pairsInfo(i,:)];
                    end
                end
            end
        end
    end
    if ~isempty(extraAssemblies{a})
        extraAssemblies{a} = unique(extraAssemblies{a});
        extAssemblies{a} = [extAssemblies{a},extraAssemblies{a}];
    end
end

% Final matrix generation
maxSize = max(cellfun(@length, extraAssemblies)); 
matExtraAssemblies = NaN(length(assemblies), maxSize);
for a = 1:length(extraAssemblies)
    extraAssembly = extraAssemblies{a};
    matExtraAssemblies(a, 1:length(extraAssembly)) = extraAssembly;
end

maxSize = max(cellfun(@length, extAssemblies)); 
matExtAssemblies = NaN(length(assemblies), maxSize);
for a = 1:length(extAssemblies)
    extAssembly = extAssemblies{a};
    matExtAssemblies(a, 1:length(extAssembly)) = extAssembly;
end
end