- 博客(4)
- 收藏
- 关注
RSS订阅原创 full2straw
%先生成H1 H1 = primary50k; H1(:,1:2) = ceil(H1(:,1:2)/50000)+1; H1(isnan( H1(:,3)),3) = 0; H1 = double(H1); H1 = sparse(H1(:,1),H1(:,2),H1(:,3)); if sum(sum(triu(H1,1)))==0 H1=tril(H1)+tril(H1,-1)'; else H1=triu(H1)+triu(H1,1).
2021-09-23 14:57:50
162
原创 判断primary10k是否一致
MATLAB H1 = primary10k; linz=(H1(:,1)==0); i=1; while (i<=40505069) if linz(i)==1 H1(i,:)=[]; break end i=i+1; end H1(:,1:2) = ceil(H1(:,1:2)/10000); H1(isnan( H1(:,3)),3) = 0; H1
2021-09-22 23:20:08
103
原创 matlab随笔
MATLAB数据结构之cell，类似于R中的list >> X=[1,2,3;4,5,6;7,8,9] X = 1 2 3 4 5 6 7 8 9 >> Y=zeros(5,1) Y = 0 0 0 0 0 >> Z={X,Y}; >> Z(1,1)%索引数据类型 ans = 1×1 cell 数组 .
2021-09-20 12:13:24
268
原创 原始数据比较
FIND 提取三列数据，数据类型为data.frame library(strawr) noxlink<-straw(norm='KR',fname='D:/xunleidownload/GSE63525_GM12878_insitu_noXlink_combined.hic', chr1loc='1', chr2loc='1', unit='BP',binsize=50000 , matrix = "observed") 转为hic矩阵,数据类型为matrix library(Hi.
2021-09-18 20:35:53
274
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人