IP属地:上海
用原始表型获得相关性系数的矩阵数据的话,则可以使用psych包的corr.test函数。导入数据,注意不要表头不要在里面,不然不能识别,而且不能...
①在E2单元格输入以下数组公式,②按Ctrl+Shift+Enter组合键结束,③然后向右向下填充公式 =INDEX(B,SMALL(IF(2:...
本脚本适用于GAPIT产生的GWAS结果: 汇总文件在tmp中,输出为csv即可。提取的仅仅是显著位点(阈值自己调整)
首先找到/etc/apt/source.list,修改成如下即可 用vim全选删除,直接粘贴,wq保存即可、 之后 安装好宝塔,愉快的玩耍吧
统计分析时候常常格式录入时候并不是R包想要的格式,如果数据量很大,用excel转换及其麻烦,可以考虑使用reshape2包中的melt函数进行转...
现有区间如上,需要统计各区间内的SNP个数 (示例文件中的Chr文件) 公司返还的,或者说从hmp等格式上截取下来的之后SNP位点,例如下图 由...
安装bioconductor这次遇到了很大的问题,一直提示若干其他的问题,发现还是网速的和网络链接的问题导致的,时好时不好,这里在在服务器中安装...
绘制染色体特定片段的LDheatmap使用LDheatmap包genetics转换数据格式
我们GWAS常用的软件GAPIT很方便,但是如果表型数据多,将最后的数据进行汇总和处理将也是一件很麻烦的事情。这里,我尝试编写了一个小脚本,可以...