GeneChart R code

This is the R code for the GeneChart assignment used in the paper. The main commands are at the top. Select the Secretary-Genechart.input.csv as stoat.dat. Read in this R code using source(file.choose()). Plot the chart using the command: chart.func(stoat.dat, refpopnames=c("Mainland", "y2005"), include.pops=c("y2006", "y2007", "y2008","y2009")) title(main="..."). This code will then produce a GeneChart assignment using the Mainland and Secretary Island 2005 captures as reference populations, and attempt to assign subsequent Secretary Island captures (from 2006 - 2009) onto this chart. The reference populations and additional included populations can be modified in the command.

本文所附为该论文中用于GeneChart分析任务的R代码。核心命令位于代码顶部。请将Secretary-Genechart.input.csv文件指定为stoat.dat。可通过source(file.choose())命令加载该R代码。执行如下命令即可绘制图表：chart.func(stoat.dat, refpopnames=c("Mainland", "y2005"), include.pops=c("y2006", "y2007", "y2008","y2009")); title(main="...")。该代码将以大陆种群（Mainland）与2005年秘书岛捕获样本作为参考种群，生成GeneChart分析任务，并尝试将2006至2009年后续采集的秘书岛捕获样本映射至该图表中。用户可通过修改该命令参数，自定义参考种群与额外纳入的种群。