原创 R qtl package 使用記錄
1. qtl安裝了很多次,提示R版本不符合;重啓了R,安裝成功 2. 畫出的圖如下,更詳細的linkage map 圖需要藉助其它軟件 ######################## install.packages("qtl"
原创 package is not available (for R version XXX)
install.packages('package_name', dependencies=TRUE, repos='http://cran.rstudio.com/') R 安裝包的提示 錯誤,常用的是添加repo: https://
原创 Fst, pi, TajimaD plink 計算
Fst, pi, TajimaD 計算 首先使用vcftools 計算出來各個各個值 vcftools --vcf test.vcf --window-pi 3000 --out test vcftools --vcf test.
原创 awk 多條件判斷
awk '{if (NR==2){print $0/4}else if(NR==4){print $0/4} else {print}}' xxx.txt
原创 /lib/ld-linux.so.2: bad ELF interpreter
原因是64位系統使用了32位的軟件,解決方法網上有很多 就是重新安裝 glibc的 i686 ,然後export LD_LIBRARY_PATH; 不過這個用一次就會無效;查詢其它網頁發現root安裝這個時候有個依賴也會一同安裝:
原创 grep -rn 查找及批量替換
grep -rn可以關鍵詞查找符合條件的文件的行;去重文件名,然後xargs替換模式 grep -rn "xxx" . |cut -d ':' -f 1 |sort |uniq |xargs sed -i 's/xxx /yyy/g'
原创 GATK Somatic 流程學習1
1. 主要是參考鏈接https://software.broadinstitute.org/gatk/documentation/article?id=9183#1 的步驟,閱讀了解相應的介紹 2. https://console.clo
原创 perl lookup issue 解決
1. 折騰了大半天: 首先重新 編譯安裝Perl,使其支持多線程: 2.tar vxf perl-5.22.0.tar.bz2 3 cd perl-5.22.0 4 ./Configure -des -Dprefix=~/.plenv
原创 bioinformatics 雜誌list
一個 很好的投稿參考,生信雜誌信息網頁 http://bioinfolab.miamioh.edu/bioinfolab/journals.php
原创 一個比較好用的RNA seq在線分析工具iDEP
iDEP 是一個比較好用的在線RNAseq分析工具,上傳raw counts 文件可以輸出heatmap, pca等,還可以做target 預測,功能強大 http://bioinformatics.sdstate.edu/idep/ h
原创 一個比較代碼差異的好網站diffchecker
一個比較代碼差異的好網站 https://www.diffchecker.com/
原创 好用的桌面日曆軟件desktopcal
一款好用的桌面日曆軟件,可以寫日程計劃: http://chs.desktopcal.com/chs/
原创 bedtools getfasta 使用
bedtools getfasta , 可以很方便的根據 bed 提取 序列
原创 GWAS 分析Fst畫圖
GWAS 分析 Fst https://rpubs.com/rossibarra/61700 有詳盡的步驟可以參考畫圖 PS 當然這個網站 https://rpubs.com/ 也有很多很好的參考資料
原创 從VCF文件畫PCA圖
1. 使用plink 獲取 SNP_test.eigenvec 文件 plink --maf 0.05 --allow-extra-chr --vcf SNP.vcf --pca header tabs -out SNP_test eig