private void CalHao(SequenceData seqX, SequenceData seqY, Result result, int k, int m = 0)
        {
            if (k < 2)
            {
                return;
            }

            if (result.ContiansKeys("Hao_k" + k))
            {
                return;
            }

            double resultHao = 0;
            double tmpXY     = 0;
            double tmpX      = 0;
            double tmpY      = 0;
            double tmp1      = 0;
            double tmp2      = 0;
            double fXi       = 0;
            double fYi       = 0;

            List <int>    listKtupleX = seqX.GetKtupleData(k);
            List <double> listPosX    = seqX.GetMarkovData(k, k - 2);
            int           totalX      = seqX.GetTotal(k);

            List <int>    listKtupleY = seqY.GetKtupleData(k);
            List <double> listPosY    = seqY.GetMarkovData(k, k - 2);
            int           totalY      = seqY.GetTotal(k);

            for (int i = 0; i < listKtupleY.Count; i++)
            {
                fXi = listKtupleX[i] * 1.0 / totalX;
                fYi = listKtupleY[i] * 1.0 / totalY;
                if (listPosX[i] == 0)
                {
                    tmp1 = -1;
                }
                else
                {
                    tmp1 = (fXi / listPosX[i]) - 1;
                }
                if (listPosY[i] == 0)
                {
                    tmp2 = -1;
                }
                else
                {
                    tmp2 = (fYi / listPosY[i]) - 1;
                }
                tmpXY += tmp1 * tmp2;
                tmpX  += tmp1 * tmp1;
                tmpY  += tmp2 * tmp2;
            }
            tmpX      = Math.Sqrt(tmpX);
            tmpY      = Math.Sqrt(tmpY);
            resultHao = (1 - tmpXY / (tmpX * tmpY)) / 2;
            SaveResult("Hao_k" + k, seqX, seqY, result, resultHao);
            //result.AddToMatrix("Hao_k" + k, seqX.ID, seqY.ID, resultHao);
            //result.AddToMatrix("Hao_k" + k, seqY.ID, seqX.ID, resultHao);
        }
Exemple #2
0
        private void CalCh(SequenceData seqX, SequenceData seqY, Result result, int k, int m = 0)
        {
            //public
            double fXi = 0;
            double fYi = 0;
            double cXi = 0;
            double cYi = 0;

            //Ch
            double resultCh = 0;
            double tmpCh    = 0;

            List <int> listKtupleX = seqX.GetKtupleData(k);
            int        totalX      = seqX.GetTotal(k);

            List <int> listKtupleY = seqY.GetKtupleData(k);
            int        totalY      = seqY.GetTotal(k);

            for (int i = 0; i < listKtupleX.Count; i++)
            {
                //public
                cXi = listKtupleX[i];
                cYi = listKtupleY[i];
                fXi = cXi / totalX;
                fYi = cYi / totalY;

                //Ch
                tmpCh    = Math.Abs(fXi - fYi);
                resultCh = tmpCh > resultCh ? tmpCh : resultCh;
            }
            //Ch
            SaveResult("Ch_k" + k, seqX, seqY, result, resultCh);
            //result.AddToMatrix("Ch_k" + k, seqX.ID, seqY.ID, resultCh);
            //result.AddToMatrix("Ch_k" + k, seqY.ID, seqX.ID, resultCh);
        }
Exemple #3
0
        private void CalD2Star(SequenceData seqX, SequenceData seqY, Result result, int k, int m)
        {
            if (k < 3)
            {
                return;
            }
            //public
            List <int>    listKtupleX = seqX.GetKtupleData(k);
            List <double> listPosX    = seqX.GetMarkovData(k, m);
            double        totalX      = seqX.GetTotal(k);

            List <int>    listKtupleY = seqY.GetKtupleData(k);
            List <double> listPosY    = seqY.GetMarkovData(k, m);
            double        totalY      = seqY.GetTotal(k);

            //D2Star
            double tmpXD2Star   = 0;
            double tmpYD2Star   = 0;
            double resultD2Star = 0;

            for (int i = 0; i < listKtupleX.Count; i++)
            {
                //public
                double cXi     = listKtupleX[i];
                double cYi     = listKtupleY[i];
                double pXi     = listPosX[i];
                double pYi     = listPosY[i];
                double tmpX    = totalX * pXi;
                double tmpY    = totalY * pYi;
                double cXi_bar = cXi - tmpX;
                double cYi_bar = cYi - tmpY;
                double tmp1    = cXi_bar * cXi_bar;
                double tmp2    = cYi_bar * cYi_bar;


                double tmp3 = Math.Sqrt(tmpX * tmpY);
                if (tmpX == 0)
                {
                    tmpX = tmp3 = 1;
                }
                if (tmpY == 0)
                {
                    tmpX = tmp3 = 1;
                }
                resultD2Star += cXi_bar * cYi_bar / tmp3;
                tmpXD2Star   += tmp1 / tmpX;
                tmpYD2Star   += tmp2 / tmpY;
            }

            //D2Star

            tmpXD2Star   = Math.Sqrt(tmpXD2Star);
            tmpYD2Star   = Math.Sqrt(tmpYD2Star);
            resultD2Star = 0.5 * (1 - resultD2Star / (tmpXD2Star * tmpYD2Star));
            SaveResult("D2Star_k" + k + "_M" + m, seqX, seqY, result, resultD2Star);
            //result.AddToMatrix("D2Star_k" + k + "_M" + m, seqX.ID, seqY.ID, resultD2Star);
            //result.AddToMatrix("D2Star_k" + k + "_M" + m, seqY.ID, seqX.ID, resultD2Star);
        }
        private void CalD2S(SequenceData seqX, SequenceData seqY, Result result, int k, int m)
        {
            if (k < 3 || k == m)
            {
                return;
            }
            //public
            List <int>    listKtupleX = seqX.GetKtupleData(k);
            List <double> listPosX    = seqX.GetMarkovData(k, m);
            double        totalX      = seqX.GetTotal(k);

            List <int>    listKtupleY = seqY.GetKtupleData(k);
            List <double> listPosY    = seqY.GetMarkovData(k, m);
            double        totalY      = seqY.GetTotal(k);


            //D2S
            double resultD2S = 0;
            double tmpD2S    = 0;

            //double tmpXD2S = 0;
            //double tmpYD2S = 0;

            for (int i = 0; i < listKtupleX.Count; i++)
            {
                //public
                double cXi     = listKtupleX[i];
                double cYi     = listKtupleY[i];
                double pXi     = listPosX[i];
                double pYi     = listPosY[i];
                double tmpX    = totalX * pXi;
                double tmpY    = totalY * pYi;
                double cXi_bar = cXi - tmpX;
                double cYi_bar = cYi - tmpY;
                double tmp1    = cXi_bar * cXi_bar;
                double tmp2    = cYi_bar * cYi_bar;



                tmpD2S = Math.Sqrt(tmp1 + tmp2);
                if (tmpD2S == 0)
                {
                    tmpD2S = 1;
                }
                resultD2S += cXi_bar * cYi_bar / tmpD2S;
                //tmpXD2S += tmp1 / tmpD2S;
                //tmpYD2S += tmp2 / tmpD2S;
            }

            //tmpXD2S = Math.Sqrt(tmpXD2S);
            //tmpYD2S = Math.Sqrt(tmpYD2S);
            //resultD2S = (1 - resultD2S / (tmpXD2S * tmpYD2S)) * 0.5;

            SaveResult("D2S_k" + k + "_M" + m, seqX, seqY, result, resultD2S);

            //result.AddToMatrix("D2S_k" + k + "_M" + m, seqX.ID, seqY.ID, resultD2S);
            //result.AddToMatrix("D2S_k" + k + "_M" + m, seqY.ID, seqX.ID, resultD2S);
        }
        private void CalEu(SequenceData seqX, SequenceData seqY, Result result, int k, int m = 0)
        {
            if (result.ContiansKeys("Eu_k" + k))
            {
                return;
            }
            //public
            double fXi = 0;
            double fYi = 0;
            double cXi = 0;
            double cYi = 0;

            //Eu
            double resultEu = 0;
            double tmpEu    = 0;


            List <int> listKtupleX = seqX.GetKtupleData(k);
            int        totalX      = seqX.GetTotal(k);

            List <int> listKtupleY = seqY.GetKtupleData(k);
            int        totalY      = seqY.GetTotal(k);

            //Eu
            for (int i = 0; i < listKtupleX.Count; i++)
            {
                //public
                cXi = listKtupleX[i];
                cYi = listKtupleY[i];
                fXi = cXi / totalX;
                fYi = cYi / totalY;

                //Eu
                tmpEu     = (fXi - fYi) * (fXi - fYi);
                resultEu += tmpEu;
            }
            //Eu
            resultEu = Math.Sqrt(resultEu);
            SaveResult("Eu_k" + k, seqX, seqY, result, resultEu);

            //result.AddToMatrix("Eu_k" + k, seqX.ID, seqY.ID, resultEu);
            //result.AddToMatrix("Eu_k" + k, seqY.ID, seqX.ID, resultEu);
        }
        private void CalMa(SequenceData seqX, SequenceData seqY, Result result, int k, int m = 0)
        {
            if (result.ContiansKeys("Ma_k" + k))
            {
                return;
            }
            //public
            double fXi = 0;
            double fYi = 0;
            double cXi = 0;
            double cYi = 0;

            //Ma
            double tmpMa    = 0;
            double resultMa = 0;


            List <int> listKtupleX = seqX.GetKtupleData(k);
            int        totalX      = seqX.GetTotal(k);

            List <int> listKtupleY = seqY.GetKtupleData(k);
            int        totalY      = seqY.GetTotal(k);

            //Ma
            for (int i = 0; i < listKtupleX.Count; i++)
            {
                //public
                cXi = listKtupleX[i];
                cYi = listKtupleY[i];
                fXi = cXi / totalX;
                fYi = cYi / totalY;

                //Ma
                tmpMa     = Math.Abs(fXi - fYi);
                resultMa += tmpMa;
            }
            //Ma
            SaveResult("Ma_k" + k, seqX, seqY, result, resultMa);

            //result.AddToMatrix("Ma_k" + k, seqX.ID, seqY.ID, resultMa);
            //result.AddToMatrix("Ma_k" + k, seqY.ID, seqX.ID, resultMa);
        }
        private void CalD2(SequenceData seqX, SequenceData seqY, Result result, int k, int m = 0)
        {
            if (result.ContiansKeys("D2_k" + k))
            {
                return;
            }
            //public
            double cXi = 0;
            double cYi = 0;

            //D2
            double resultD2 = 0;
            double tmpXD2   = 0;
            double tmpYD2   = 0;
            double tmpD2    = 0;

            List <int> listKtupleX = seqX.GetKtupleData(k);
            List <int> listKtupleY = seqY.GetKtupleData(k);

            for (int i = 0; i < listKtupleX.Count; i++)
            {
                //public
                cXi = listKtupleX[i];
                cYi = listKtupleY[i];

                //D2
                tmpD2  += cYi * cXi;
                tmpXD2 += cXi * cXi;
                tmpYD2 += cYi * cYi;
            }
            //D2
            tmpXD2   = Math.Sqrt(tmpXD2);
            tmpYD2   = Math.Sqrt(tmpYD2);
            tmpD2    = tmpD2 / (tmpXD2 * tmpYD2);
            resultD2 = 0.5 * (1 - tmpD2);

            SaveResult("D2_k" + k, seqX, seqY, result, resultD2);

            //result.AddToMatrix("D2_k" + k, seqX.ID, seqY.ID, resultD2);
            //result.AddToMatrix("D2_k" + k, seqY.ID, seqX.ID, resultD2);
        }
Exemple #8
0
        private void CalEuMaChD2(SequenceData seqX, SequenceData seqY, Result result, int k, int m = 0)
        {
            //Eu
            double resultEu = 0;
            double tmpEu    = 0;

            //Ma
            double tmpMa    = 0;
            double resultMa = 0;

            //public
            double fXi = 0;
            double fYi = 0;
            double cXi = 0;
            double cYi = 0;

            //Ch
            double resultCh = 0;

            //D2
            double resultD2 = 0;
            double tmpXD2   = 0;
            double tmpYD2   = 0;
            double tmpD2    = 0;

            List <int> listKtupleX = seqX.GetKtupleData(k);
            int        totalX      = seqX.GetTotal(k);

            List <int> listKtupleY = seqY.GetKtupleData(k);
            int        totalY      = seqY.GetTotal(k);

            for (int i = 0; i < listKtupleX.Count; i++)
            {
                //public
                cXi = listKtupleX[i];
                cYi = listKtupleY[i];
                fXi = cXi / totalX;
                fYi = cYi / totalY;

                //D2
                tmpD2  += cYi * cXi;
                tmpXD2 += cXi * cXi;
                tmpYD2 += cYi * cYi;

                //Eu
                tmpEu     = (fXi - fYi) * (fXi - fYi);
                resultEu += tmpEu;

                //Ma
                tmpMa     = Math.Abs(fXi - fYi);
                resultMa += tmpMa;

                //Ch
                resultCh = tmpMa > resultCh ? tmpMa : resultCh;
            }
            //D2
            tmpXD2   = Math.Sqrt(tmpXD2);
            tmpYD2   = Math.Sqrt(tmpYD2);
            tmpD2    = tmpD2 / (tmpXD2 * tmpYD2);
            resultD2 = 0.5 * (1 - tmpD2);

            //Eu
            resultEu = Math.Sqrt(resultEu);

            SaveResult("Eu_k" + k, seqX, seqY, result, resultEu);
            SaveResult("Ch_k" + k, seqX, seqY, result, resultCh);
            SaveResult("Ma_k" + k, seqX, seqY, result, resultMa);
            SaveResult("D2_k" + k, seqX, seqY, result, resultD2);

            //result.AddToMatrix("Eu_k" + k, seqX.ID, seqY.ID, resultEu);
            //result.AddToMatrix("Ch_k" + k, seqX.ID, seqY.ID, resultCh);
            //result.AddToMatrix("Ma_k" + k, seqX.ID, seqY.ID, resultMa);
            //result.AddToMatrix("D2_k" + k, seqX.ID, seqY.ID, resultD2);

            //result.AddToMatrix("Eu_k" + k, seqY.ID, seqX.ID, resultEu);
            //result.AddToMatrix("Ch_k" + k, seqY.ID, seqX.ID, resultCh);
            //result.AddToMatrix("Ma_k" + k, seqY.ID, seqX.ID, resultMa);
            //result.AddToMatrix("D2_k" + k, seqY.ID, seqX.ID, resultD2);
        }