MOEAD_Net_getFinalCommunity_multiError.m

function  MOEAD_Net_getFinalCommunity_multiError(path,name,real_path,overlapping,isRealWorld,c)
%%%%%
global edgeMatrix 
global AdjMatrix 
global edgeNum
global degree

%%%%%
name1=name;
AdjMatrix = [];
tic
%% 导入邻接矩阵
%% 自动根据读取的网络数据格式确定邻接矩阵和有无真是划分%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 hasReal =  exist(real_path,'file');%% 有无真实划分 大于0，有
 networkData = sprintf('adjMatrix_coreNodes/%s/%s.mat',name1,name1);
 hasNetWorkAdj =  exist(networkData,'file');
if hasNetWorkAdj
    load(networkData);
else
    AdjMatrix = single(load(path));
end
networkData2 = sprintf('adjMatrix_coreNodes/%s/coreNodes.mat',name1);
load(networkData2);
if isempty(AdjMatrix)||size(AdjMatrix,2) == 2
%      AdjMatrix = single(load(path));
         AdjMatrix_size = size(AdjMatrix);
         if AdjMatrix_size(2)>2 %%邻接矩阵表示
        %     edgeNum = sum(sum(AdjMatrix));
            [edgeMatrix1,edgeMatrix2] = (find(AdjMatrix==1));
            edgeMatrix = [edgeMatrix2,edgeMatrix1];
         else                   %%边表表示
             edgeMatrix = AdjMatrix;
             needAddOne = 0;    %%是否需要加1
             numVar=(max(max(AdjMatrix(:,1)),max(AdjMatrix(:,2))));
             if find(AdjMatrix==0)>0 %% 从0开始编号
                  needAddOne = 1;
                  numVar=numVar+1;
             end
             edgeNum = AdjMatrix_size(1);
              AdjMatrix = Adjreverse(AdjMatrix,numVar,needAddOne);
         end
end

 numVar=single(size(AdjMatrix,1));
 if hasReal >0
     if overlapping == 0 %%非重叠
         Datalabel=(load(real_path));
         if size(Datalabel,2)==2  %%社团划分为“点--》社团”的2列形式
             Datalabel=(Datalabel(:,2)');
         end
     else                %%重叠
         if isRealWorld ==1 %%真实网络
              Datalabel=(load(real_path));
            realCommunity = label2community(load(real_path));
         else
            [realCommunity,~,~] = LFR_community2community(real_path);
            for k = 1:length(realCommunity)
                Datalabel(1,realCommunity{k}) = k;
            end
         end
     end
 else
     if overlapping == 0 %%非重叠
         Datalabel= false(1,numVar);
     else                %%重叠
         realCommunity = {};
     end
     
     
 end
%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 clear EDGE EDGE_Data edgeMatrix1 edgeMatrix2;
 
%  degree=sum(AdjMatrix,1);
%  index=find(degree==0);
%  AdjMatrix(index,:)=[];
%  AdjMatrix(:,index)=[];
 
% degree = single(sum(AdjMatrix));
% avgDegree = floor(sum(degree)/length(AdjMatrix));
degree=single(sum(AdjMatrix,2));
edgeNum=sum(degree)/2;

community = {};
networkDataRoot = sprintf('results/%s/coreCommunityFirstFront/coreCommunity_%s/',name1,num2str(c));%%获取第c次扩充的核心点
    path = networkDataRoot;  
    fileExt = '*.mat';  
    files = dir(fullfile(path,fileExt));  
    for i=1:size(files,1)  
        fileName = strcat(path,files(i,1).name);  
        coreCommunity = load(fileName);
%         community = [community {coreCommunity.coreCommunity}];
community = coreCommunity.cellCommunity;
    end; 
    
     pathTime = sprintf('results/%s/time/postProcessTime/',name1);
    if ~isdir(pathTime) %判断路径是否存在
        mkdir(pathTime);
    end
    pathResult = sprintf('results/%s/Q_NMI/',name1);
    if ~isdir(pathResult) %判断路径是否存在
        mkdir(pathResult);
    end
     pathResultError = sprintf('results/%s/errorNodes/',name1);
    if ~isdir(pathResultError) %判断路径是否存在
        mkdir(pathResultError);
    end
    
    clock1 = clock;
    
    
    
    
    
    [communityIndex,bestQ] = run_EA(name,community);

    findNodes = [];
    for ii = 1:length(communityIndex)
        myCommunity{ii} = community{ii}{communityIndex(ii)}(1,1:end-2);
        findNodes = unique([findNodes myCommunity{ii}]);
    end
    save myCommunity;
    errorNodes = zeros(1,length(myCommunity),'single');
    for oo = 1:length(myCommunity)
        errorNodes(oo) = length(unique(Datalabel(1,myCommunity{oo})))-1;%%与真实结果的误差数量
    end
    
     avgErrorNodes = sum(errorNodes)/length(myCommunity);
     pathResultError=sprintf('results/%s/errorNodes/errorNodes(%s_%s=%s)_%s.txt',name1,num2str(sum(errorNodes)),num2str(length(myCommunity)),num2str(avgErrorNodes),num2str(c));
     savedata1(pathResultError,errorNodes);
    
    
    
    myLabel = zeros(1,length(Datalabel),'single');
    for ii = 1:length(myCommunity)
        myLabel(myCommunity{ii}) = ii;
    end
    modnmi0 = [];
    modnmi0(1,1) = bestQ ;%%不考虑离散点直接计算
    if sum(Datalabel) ~= 0
      modnmi0(1,2) = nmi(Datalabel,myLabel(1,:));
    else
      modnmi0(1,2) = 0;
    end



remainNodes = setdiff(1:length(AdjMatrix),findNodes);


%% 对剩余的节点进行分派，选取Q最大的
% for i = 1:length(remainNodes)
%     i_Q = [];
%     for j = 1:length(myCommunity)
%         myLabel(1,remainNodes(1,i)) = j;
%         i_Q = [i_Q modularity(myLabel,AdjMatrix)];
%     end
%     [maxQ,index] = max(i_Q);
%      myLabel(1,remainNodes(1,i)) = index;
%      myCommunity{index} = unique([ myCommunity{index} remainNodes(1,i)]);
% end
%% 查找邻居所在团
% for i = 1:length(remainNodes)
%     i_Q = [];
%     for j = 1:length(myCommunity)
%         if ~isempty(find(AdjMatrix(remainNodes(1,i),myCommunity{j})>0)) %%有邻居所在团
%             myLabel(1,remainNodes(1,i)) = j;%%加入邻居所在团
%             i_Q = [i_Q;j,modularity(myLabel,AdjMatrix)];
%         end
%     end
%     if ~isempty(i_Q)
%          [maxQ,index] = max(i_Q(:,2));
%         myLabel(1,remainNodes(1,i)) = i_Q(index,1);
%         myCommunity{i_Q(index,1)} = unique([ myCommunity{i_Q(index,1)} remainNodes(1,i)]);
%     end
%    
% end
%%%%%%%%%加入最多邻居所在社团
% % for i = 1:length(remainNodes)
% %     for j = 1:length(myCommunity)
% %         n = find(AdjMatrix(remainNodes(i),:)>0);
% %         n_c_num = zeros(1,length(myCommunity));
% %         if ~isempty(n)
% %             for k = 1:length(n)
% %                 if ismember(n(k),myCommunity{j})
% %                     n_c_num(j) =  n_c_num(j)+1;
% %                 end
% %             end
% %         end
% %     end
% %     [~,remain_i_label] = max(n_c_num);
% %     myLabel(1,remainNodes(1,i)) = remain_i_label;%%加入邻居所在团
% % end


% tic
%     for i = 1:length(remainNodes)
%         i_nebor = find(AdjMatrix(remainNodes(i),:)>0);
%         i_nebor_label = myLabel(i_nebor);
%         i_nebor_label(i_nebor_label==0) =[];
%         if ~isempty(i_nebor_label)
%             X = unique(i_nebor_label);
%             [M,N]=hist(i_nebor_label,X);
%             [~,index] = max(M);
%             myLabel(remainNodes(i)) = N(index);
%         else
%            myLabel(remainNodes(i)) = randi(length(myCommunity)); 
%         end
%     end
% toc











topTree = {};
 finalCommunity = {};
     treeIndex = 1;
     topTree{treeIndex} = myCommunity;

         COR = zeros(length(myCommunity),length(myCommunity),'single');
         treeIndex = 1;
         for top_j = 1:length(myCommunity)-1
             for top_k = top_j+1:length(myCommunity)
                 COR(top_j,top_k) = length(intersect(myCommunity{top_j},myCommunity{top_k}))/min(length(myCommunity{top_j}),length(myCommunity{top_k}));
                 if COR(top_j,top_k)>0.4
                     myCommunity{top_j} = union(myCommunity{top_j},myCommunity{top_k});
                     myCommunity{top_k} = [];
                     
                 end
             end
         end
         myCommunity(cellfun(@isempty,myCommunity))=[];
         topTree{treeIndex} = myCommunity;

     maxQ = 0;
     Q = [];
     for tt = 1:length(topTree)
         myLabel = zeros(1,length(AdjMatrix),'single');
         for yy = 1:length(topTree{tt})
            myLabel(topTree{tt}{yy}) = yy;
         end
         Q = [Q,modularity(myLabel,AdjMatrix,degree,edgeNum)];
     end
     [~,sortIndex] = sort(Q,'descend');
     
     [maxQ,maxQ_index] = max(Q);
     
     if length(sortIndex)>=10
         len = 10;
     else
         len = length(sortIndex);
     end


     tic
     tt = 0;
remainNodes2 = [];
     remainNodes3 = [];
    for i = 1:length(remainNodes)
        i_nebor = find(AdjMatrix(remainNodes(i),:)>0);
        i_nebor_label = myLabel(i_nebor);
        i_nebor_label(i_nebor_label==0) =[];
        if ~isempty(i_nebor_label)
            X = unique(i_nebor_label);
            [M,N]=hist(i_nebor_label,X);
            [~,index] = max(M);
            myLabel(remainNodes(i)) = N(index);
        else

        remainNodes2 = [remainNodes2 remainNodes(i)];
        end
    end
    
    for i = 1:length(remainNodes2)
        i_nebor = find(AdjMatrix(remainNodes2(i),:)>0);
        i_nebor_label = myLabel(i_nebor);
        i_nebor_label(i_nebor_label==0) =[];
        if ~isempty(i_nebor_label)
            X = unique(i_nebor_label);
            [M,N]=hist(i_nebor_label,X);
            [~,index] = max(M);
            myLabel(remainNodes2(i)) = N(index);
        else

         remainNodes3 = [remainNodes3 remainNodes2(i)];
        end
    end
    
    for i = 1:length(remainNodes3)
        i_nebor = find(AdjMatrix(remainNodes3(i),:)>0);
        i_nebor_label = myLabel(i_nebor);
        i_nebor_label(i_nebor_label==0) =[];
        if ~isempty(i_nebor_label)
            X = unique(i_nebor_label);
            [M,N]=hist(i_nebor_label,X);
            [~,index] = max(M);
            myLabel(remainNodes3(i)) = N(index);
        else
            tt = tt+1;
           myLabel(remainNodes3(i)) = randi(length(myCommunity)); 
        end
    end
toc
     

root1 = sprintf('results/%s/remainNodes',name1);
if ~isdir(root1) %判断路径是否存在
    mkdir(root1);
end
 path1=sprintf('results/%s/remainNodes/remainNodes(remainNodesNum=%s)_%s.txt',name1,num2str(length(remainNodes)),num2str(c));
savedata1(path1,remainNodes);
root2 = sprintf('results/%s/addRemainNodesCommunity',name1);
if ~isdir(root2) %判断路径是否存在
    mkdir(root2);
end
 path2=sprintf('results/%s/addRemainNodesCommunity/addRemainNodesCommunity_%s.txt',name1,num2str(c));
savedata1(path2,myLabel);





     modnmi1 = [];
     for i = 1:len
         finalCommunity{i} = topTree{sortIndex(i)};
         qqq = single(modularity(myLabel(i,:),AdjMatrix,degree,edgeNum));%%分派离散点后直接计算Q
        modnmi1(i,1) = qqq ;
        if sum(Datalabel) ~= 0
          modnmi1(i,2) = nmi(Datalabel,myLabel(i,:));
       else
          modnmi1(i,2) = 0;
       end
         
         allNodes = cell2mat(finalCommunity{i});
         A=tabulate(allNodes);
% % % %          errorNodes=find(A(:,2)>1);
         
          
           [myLabel,Q2]=find_error2(single(myLabel),AdjMatrix,qqq,degree,edgeNum);
          
          modnmi = zeros(size(myLabel,1),2,'single');
            for k = 1:size(myLabel,1)
            modnmi(k,1) = (modularity(myLabel(k,:),AdjMatrix,degree,edgeNum)) ;%%纠错后计算
               if sum(Datalabel) ~= 0
                  modnmi(k,2) = nmi(Datalabel,myLabel(k,:));
              else
                  modnmi(k,2) = 0;
              end
           end
            modnmi

     end
     
     clock2 = clock;
     etime(clock2,clock1);
     pathTime=sprintf('results/%s/time/postProcessTime/postProcessTime_%s.txt',name1,num2str(c));
     savedata1(pathTime, etime(clock2,clock1));
     
    
% % % % %      labelIndex = unique(myLabel);
% % % % %      finalCommunityNum = length(labelIndex);
% % % % %      COMMUNITY = {};
% % % % %      for iii = 1:finalCommunityNum
% % % % %          COMMUNITY{iii} = find(myLabel == labelIndex(iii));
% % % % %      end
% % % % %      
% % % % %       afterErrorCorrectionErrorNodes = zeros(1,finalCommunityNum,'single');
% % % % %     for oo = 1:finalCommunityNum
% % % % %         afterErrorCorrectionErrorNodes(oo) = length(unique(Datalabel(1,COMMUNITY{oo})))-1;%%与真实结果的误差数量
% % % % %     end
% % % % %     
% % % % %      avgErrorNodes = sum(afterErrorCorrectionErrorNodes)/finalCommunityNum;%%每个团的错误节点数
% % % % %      
% % % % %      %%所有错误节点总数除以社团数
% % % % %      pathResultError2=sprintf('results/%s/errorNodes/afterErrorCorrectionErrorNodes(%s_%s=%s)_%s.txt',name1,num2str(sum(afterErrorCorrectionErrorNodes)),num2str(finalCommunityNum),num2str(avgErrorNodes),num2str(c));
% % % % %      savedata1(pathResultError2,afterErrorCorrectionErrorNodes);
% % % % %      

   root3 = sprintf('results/%s/findErrorCommunity',name1);
    if ~isdir(root3) %判断路径是否存在
        mkdir(root3);
    end
 path3=sprintf('results/%s/findErrorCommunity/findErrorCommunity(cNum=%s)_%s.txt',name1,num2str(length(unique(myLabel))),num2str(c));
savedata1(path3,myLabel);
pathResult = sprintf('results/%s/Q_NMI/Q_NMI_%s.txt',name1,num2str(c));
result = [modnmi0;modnmi1;modnmi];
savedata1(pathResult,[result;0 0 ;mean(result(:,1)) mean(result(:,2));max(result(:,1)) max(result(:,2))]);

 end 
 


function [currentLabel,Q2]=find_error2(currentLabel,AdjMatrix,Q,degree,edgeNum)
    Q2 = Q;
%         M=[];
        
%         numVar=size(AdjMatrix,1);
%             change_node=[];
%             lable= currentLabel;
            %% %%后处理：2种方式-最大的Q lable= ParetoFront1(in,1:numVar);
            
            
            T = 5;
           while T>0
            originalLabel = currentLabel;
            communityNum=unique(currentLabel);
            for i=1:length(communityNum)
                currentCommunity = find(currentLabel==communityNum(1,i));
                if(length(currentCommunity)>=3)
                    for j=1:length(currentCommunity)
%                         m=CLique{i}(j);
%                         k=lable(m);
% k = i;
%                         index_last=find(lable==i);
%                         A=setdiff(index_last,m);
                        neighbors=find(AdjMatrix(currentCommunity(1,j),:));
                        com_max=multi_label(neighbors,currentLabel);
                        currentLabel(1,currentCommunity(1,j))=com_max;
%                         if com_max~=i
%                             change_node=[change_node{in} m];
%                         end
                    end
                end
            end
            QQ = modularity(currentLabel,AdjMatrix,degree,edgeNum);
            if QQ>Q                                                      %有疑问
%                 M=[M in];
               
                Q2 = QQ;
%                 finalLabel = currentLabel;
            else
               currentLabel =  originalLabel;
            end
             T = T-1;
           end
%         error_node=[];
%         for i=1:length(M)
%             error_node=[error_node change_node{M(i)}];
%         end
%         remove=unique(single(error_node));
        
    end