基因数据处理121之SSW的score matrix调整，使得与SparkSW评分一致

其他 2018-09-02 22:42:00 阅读次数: 0
版权声明：本文为博主原创文章，欢迎转载和交流！源码交流：https://github.com/xubo245/ https://blog.csdn.net/bob601450868/article/details/80589856

更多代码请见：https://github.com/xubo245

基因数据处理系列

1.解释

SSW的评分矩阵是128*128的，是按char的int值来进行计算的。而blosum50是蛋白质的，而且不是按ABC顺序来的，所以需要转换。

java中无法运行128*128的String矩阵，所以需要提供blosum的转换

2.代码：

    public static int[][] blosum50() {
        int[][] a = {

        };
        return a;
    }

    public static int[][] blosum50(String[][] strArr) {
        int[][] result = new int[128][128];

        int m = strArr.length;
        int n = strArr[0].length;

        for (int i = 0;i< result.length;i++) {
            for (int j = 0;j<  result[i].length;j++) {
                result[i][j] = Integer.valueOf(strArr[m - 1][n - 1]);
            }
        }
        for (int i = 0 ;i< m;i++) {
             char c = strArr[i][0].charAt(0);
            for (int j = 0 ;j<  result.length;j++) {
                result[Integer.valueOf(c)][j] = Integer.valueOf(strArr[i][strArr[i].length - 1]);
                result[j][Integer.valueOf(c)] = Integer.valueOf(strArr[i][strArr[i].length - 1]);
            }
        }

        //    println(result(0)(0))
        for (int i = 0 ;i< m;i++) {
            char strI = strArr[i][0].charAt(0);
            //       strArr(i)(0).toCharArray

            int valueI = Integer.valueOf(strArr[i][0].charAt(0));
            //      println(strI + "=" + valueI + "=" + strArr(i)(0).charAt(0).toInt)

            for (int j = 1 ;j< n;j++) {

                int valueJ =Integer.valueOf(strArr[j - 1][0].charAt(0));
                int tmp = Integer.valueOf(strArr[i][j]);
                result[valueI][valueJ] = tmp;
                //        arr(value)(j)

            }
        }

        //println(result(0)(0))
        return result;
    }
    public static String[][] blosum50Base() {
        String[][] a = {
                {"A","5","-2","-1","-2","-1","-1","-1","0","-2","-1","-2","-1","-1","-3","-1","1","0","-3","-2","0","-2","-2","-1","-1","-5"},
                {"R","-2","7","-1","-2","-4","1","0","-3","0","-4","-3","3","-2","-3","-3","-1","-1","-3","-1","-3","-1","-3","0","-1","-5"},
                {"N","-1","-1","7","2","-2","0","0","0","1","-3","-4","0","-2","-4","-2","1","0","-4","-2","-3","5","-4","0","-1","-5"},
                {"D","-2","-2","2","8","-4","0","2","-1","-1","-4","-4","-1","-4","-5","-1","0","-1","-5","-3","-4","6","-4","1","-1","-5"},
                {"C","-1","-4","-2","-4","13","-3","-3","-3","-3","-2","-2","-3","-2","-2","-4","-1","-1","-5","-3","-1","-3","-2","-3","-1","-5"},
                {"Q","-1","1","0","0","-3","7","2","-2","1","-3","-2","2","0","-4","-1","0","-1","-1","-1","-3","0","-3","4","-1","-5"},
                {"E","-1","0","0","2","-3","2","6","-3","0","-4","-3","1","-2","-3","-1","-1","-1","-3","-2","-3","1","-3","5","-1","-5"},
                {"G","0","-3","0","-1","-3","-2","-3","8","-2","-4","-4","-2","-3","-4","-2","0","-2","-3","-3","-4","-1","-4","-2","-1","-5"},
                {"H","-2","0","1","-1","-3","1","0","-2","10","-4","-3","0","-1","-1","-2","-1","-2","-3","2","-4","0","-3","0","-1","-5"},
                {"I","-1","-4","-3","-4","-2","-3","-4","-4","-4","5","2","-3","2","0","-3","-3","-1","-3","-1","4","-4","4","-3","-1","-5"},
                {"L","-2","-3","-4","-4","-2","-2","-3","-4","-3","2","5","-3","3","1","-4","-3","-1","-2","-1","1","-4","4","-3","-1","-5"},
                {"K","-1","3","0","-1","-3","2","1","-2","0","-3","-3","6","-2","-4","-1","0","-1","-3","-2","-3","0","-3","1","-1","-5"},
                {"M","-1","-2","-2","-4","-2","0","-2","-3","-1","2","3","-2","7","0","-3","-2","-1","-1","0","1","-3","2","-1","-1","-5"},
                {"F","-3","-3","-4","-5","-2","-4","-3","-4","-1","0","1","-4","0","8","-4","-3","-2","1","4","-1","-4","1","-4","-1","-5"},
                {"P","-1","-3","-2","-1","-4","-1","-1","-2","-2","-3","-4","-1","-3","-4","10","-1","-1","-4","-3","-3","-2","-3","-1","-1","-5"},
                {"S","1","-1","1","0","-1","0","-1","0","-1","-3","-3","0","-2","-3","-1","5","2","-4","-2","-2","0","-3","0","-1","-5"},
                {"T","0","-1","0","-1","-1","-1","-1","-2","-2","-1","-1","-1","-1","-2","-1","2","5","-3","-2","0","0","-1","-1","-1","-5"},
                {"W","-3","-3","-4","-5","-5","-1","-3","-3","-3","-3","-2","-3","-1","1","-4","-4","-3","15","2","-3","-5","-2","-2","-1","-5"},
                {"Y","-2","-1","-2","-3","-3","-1","-2","-3","2","-1","-1","-2","0","4","-3","-2","-2","2","8","-1","-3","-1","-2","-1","-5"},
                {"V","0","-3","-3","-4","-1","-3","-3","-4","-4","4","1","-3","1","-1","-3","-2","0","-3","-1","5","-3","2","-3","-1","-5"},
                {"B","-2","-1","5","6","-3","0","1","-1","0","-4","-4","0","-3","-4","-2","0","0","-5","-3","-3","6","-4","1","-1","-5"},
                {"J","-2","-3","-4","-4","-2","-3","-3","-4","-3","4","4","-3","2","1","-3","-3","-1","-2","-1","2","-4","4","-3","-1","-5"},
                {"Z","-1","0","0","1","-3","4","5","-2","0","-3","-3","1","-1","-4","-1","0","-1","-2","-2","-3","1","-3","5","-1","-5"},
                {"X","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-1","-5"},
                {"*","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","-5","1"}
        };
        return a;
    }

blosum50来源于https://www.ncbi.nlm.nih.gov/IEB/ToolBox/C_DOC/lxr/source/data/BLOSUM50

3.结果：

调整后与SparkSW运行基本一致，但是ref Data比较大时会不同，比如：待分析原因解决

SparkSSW
time:10
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)

记录

hadoop@Master:~/disk2/xubo/project/alignment/SparkSW/SparkSW20161114/alluxio-1.3.0$ tail -f sparkSWSSWtime201611270135.txt 
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
time:4
(1734,UniRef100_Q6GZX4)                                                         
time:5
(1734,UniRef100_Q6GZX4)                                                         
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
time:8
(1734,UniRef100_Q6GZX4)                                                         
time:9
(1734,UniRef100_Q6GZX4)                                                         
time:10
(1734,UniRef100_Q6GZX4)                                                         
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 01:49:22 WARN QueuedThreadPool: 1 threads could not be stopped
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
time:4
(1734,UniRef100_Q6GZX4)                                                         
time:5
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 01:51:26 WARN QueuedThreadPool: 4 threads could not be stopped
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
time:8
(1734,UniRef100_Q6GZX4)                                                         
time:9
(1734,UniRef100_Q6GZX4)                                                         
time:10
(1734,UniRef100_Q6GZX4)                                                         
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
time:4
(1734,UniRef100_Q6GZX4)                                                         
time:5
(1734,UniRef100_Q6GZX4)                                                         
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
time:8
(1734,UniRef100_Q6GZX4)                                                         
time:9
(1734,UniRef100_Q6GZX4)                                                         
time:10
(1734,UniRef100_Q6GZX4)                                                         
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 02:19:43 WARN QueuedThreadPool: 2 threads could not be stopped
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
time:4
(1734,UniRef100_Q6GZX4)                                                         
time:5
(1734,UniRef100_Q6GZX4)                                                         
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
time:8
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 02:25:02 WARN QueuedThreadPool: 2 threads could not be stopped
time:9
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 02:25:48 WARN QueuedThreadPool: 5 threads could not be stopped
time:10
(1734,UniRef100_Q6GZX4)                                                         
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
time:4
(1734,UniRef100_Q6GZX4)                                                         
time:5
(1734,UniRef100_Q6GZX4)                                                         
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
time:8
(1734,UniRef100_Q6GZX4)                                                         
time:9
(1734,UniRef100_Q6GZX4)                                                         
time:10
(1734,UniRef100_Q6GZX4)                                                         
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 03:12:40 WARN QueuedThreadPool: 2 threads could not be stopped
time:2
(1734,UniRef100_Q6GZX4)                                                         
time:3
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 03:15:08 WARN QueuedThreadPool: 4 threads could not be stopped
time:4
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 03:16:24 WARN QueuedThreadPool: 1 threads could not be stopped
time:5
(1734,UniRef100_Q6GZX4)                                                         
time:6
(1734,UniRef100_Q6GZX4)                                                         
time:7
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 03:20:13 WARN QueuedThreadPool: 2 threads could not be stopped
time:8
(1734,UniRef100_Q6GZX4)                                                         
time:9
(1734,UniRef100_Q6GZX4)                                                         
16/11/27 03:22:45 WARN QueuedThreadPool: 4 threads could not be stopped
time:10
(1734,UniRef100_Q6GZX4)                                                         
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:10
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:10
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:10
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:10
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
SparkSW:
time:1
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
time:10
(1734,UniRef100_Q6GZX4)                                                         
(295,UniRef100_Q91FP2)
SSW
time:1
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:2
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:3
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:4
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:5
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:6
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:7
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:8
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
time:9
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
16/11/27 12:57:11 WARN QueuedThreadPool: 1 threads could not be stopped
time:10
(1734,UniRef100_Q6GZX4)                                                         
(335,UniRef100_Q91FP2)
(291,UniRef100_Q196Y1)
SparkSW:
time:1
[Stage 1:====================================>                  (86 + 16) / 128]Exception in thread "main" org.apache.spark.SparkException: Job cancelled because SparkContext was shut down

参考

    【1】https://github.com/xubo245
    【2】http://blog.csdn.net/xubo245/

猜你喜欢

转载自blog.csdn.net/bob601450868/article/details/80589856
今日推荐
《美国对全球网络空间安全与发展的威胁和破坏》报告发布
火速冲上 GitHub 热榜 —— 开源编程语言、框架哪有这么可爱？
北京人形机器人创新中心发布全球首个纯电驱拟人奔跑的全尺寸人形机器人“天工”
LFOSSA 源来如此公开课 | 掌握云原生未来：CNCF 认证全面攻略与备考秘籍
周排行
循环神经网络（rnn）讲解
Tigao教程四：单独的关节运动
金蝶K3WISE15.0-注册套打教程
如何在Mac上配置Kubernetes
Android应用结束自身进程的方法
SpringMVC学习 十三 拦截器栈
中国驻洛杉矶总领馆举行新春招待会
HttpClient get post 发送
11 - three.js 笔记 - 绘制三维字体模型
Mysql递归获取某个父节点下面的所有子节点和子节点上的所有父节点
每日归档
更多
2024-05-01(4)
2024-04-30(1)
2024-04-29(40)
2024-04-28(0)
2024-04-27(56)
2024-04-26(39)
2024-04-25(22)
2024-04-24(36)
2024-04-23(26)
2024-04-22(39)

代码天地 © 2018 codetd.com

境外服务器   最新电视剧2022