如何计算序列间的种内和种间遗传距离 How to Calculate Max-intra and Min-inter Distance distribution

step 1. Prepare the distance matrix.

Tool: DistanceMatrix In DECIPHER(R package)


library(DECIPHER)

DNAStringSet=readDNAStringSet("S:/xxx.txt")

a= DistanceMatrix(DNAStringSet,

type = "matrix",

includeTerminalGaps = FALSE,

penalizeGapLetterMatches = FALSE,

penalizeGapGapMatches = FALSE,

correction = "none",

processors = NULL,

verbose = TRUE)

write.csv(a, file = "S:/xxxxxxx.csv")


step 2. Search the intra- and inter Distances of all sequences

Tool: Excel

Use MAXIF and MINIF formula to get intra- and inter Distances for each sequence;

 

step 3. Search Max-intra and Min-inter Distances for every species

Tools: Excel, SpeciesIdentify in TaxonDNA

使用数组公式
使用excel数组公式

Use Species Summary Function to get species name list;

Use VLOOKUP formula to search Max-intra and Min-inter distance for every species (Notice that use Species name list as Lookup_value,and sort the previous intra- and inter Distances results by descending and scending to get the final maximum and minimum results;

step 4. Draw scatter diagram.

The x-coordinate is Max-intra distane,The y-coordinateMin-inter distance.


最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容