- 博客(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
143
原创 判断primary10k是否一致
MATLABH1 = 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
90
原创 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
248
原创 原始数据比较
FIND提取三列数据，数据类型为data.framelibrary(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矩阵,数据类型为matrixlibrary(Hi.
2021-09-18 20:35:53
230
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人