240 发简信
IP属地:北京
  • @铁线莲 你好,请问stairway plot中的参数L是如何计算的呢,我一直没弄明白

    种群动态历史评估方法之 SFS(Site/Allele Frequency Spectrum)

    SFS的基本理解、介绍的一篇文章:https://theg-cat.com/tag/joint-sfs/[https://theg-cat.com/tag/joint-sfs...

  • 120
    「群体遗传学实战」第一课: 对SNP位点进行注释

    我们用于实战的数据集来自于2019年发表于NG的西瓜文章,它提供了GATK过滤后的SNP数据集用于分析,并且附录提供了完整的样本信息。该SNP数据集包括414个西瓜,2000...

  • 你好,老师!最近在研究你这个脚本的用法,可以我总是一运行就出现问题。
    The drictory is: D:\worksoftware\callCodeml-main\paml
    The treeFile is: D:\worksoftware\callCodeml-main\31.nwk
    0%| | 0/78 [00:00<?, ?it/s]
    2021-08-11 00:26:27 Staring ModelA null of accD_mafft ...
    0%| | 0/78 [00:00<?, ?it/s]
    Traceback (most recent call last):
    File "callCodeml.py", line 191, in <module>
    main()
    File "callCodeml.py", line 178, in main
    run(path_seq, name)
    File "callCodeml.py", line 140, in run
    res_null = call_codeml('null', path_seq, name)
    File "callCodeml.py", line 117, in call_codeml
    b = subprocess.Popen(['codeml', f'{path0}/WorkingDrictory_{start_time}/{name}/{type}/{type}.ctl'], stdout=subprocess.PIPE)
    File "C:\Users\yangz\anaconda3\lib\subprocess.py", line 858, in __init__
    self._execute_child(args, executable, preexec_fn, close_fds,
    File "C:\Users\yangz\anaconda3\lib\subprocess.py", line 1311, in _execute_child
    hp, ht, pid, tid = _winapi.CreateProcess(executable, args,
    FileNotFoundError: [WinError 2] 系统找不到指定的文件。

    能否指教什么问题,找不到指定文件?

    选择压力分析,从序列下载到结果绘图

    选择压力及基因在进化过程中所受到的压力,也称自然选择。 在遗传学中, Ka/Ks或者dN/dS表示的是非同义突变频率(Ka)和同义突变频率(Ks)之间的比例。 一般认为,同义...