c++ 生成模拟测序数据代码

统计真实数据不同测序位置碱基的错误率,引入到DNA片段中,从而模拟生成DNA测序数据。

参考 ART: a next-generation sequencing read simulator - PMC

#include <iostream>
#include <random>
#include <string>
#include <cassert>using namespace std;#define SEQ_LEN 36inline double r_prob(){return rand()/(double)RAND_MAX;
}static char rand_base(){short base=(short)ceil(r_prob()*4);//std::cout << "r_prob()*4 " << r_prob()*4 << std::endl;//std::cout << "ceil(r_prob()*4) " << ceil(r_prob()*4) << std::endl;switch(base){case 1:return 'A';case 2:return 'C';case 3:return 'G';case 4:return 'T';default:

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/diannao/6731.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

05_G1垃圾收集器

G1垃圾收集器简介 垃圾优先 Garbage-First&#xff08;G1&#xff09;垃圾收集器面向多处理器机器&#xff0c;适用于大内存场景。它尝试在无需太多配置的情况下实现垃圾收集暂停时间目标&#xff0c;并同时实现高吞吐量。G1旨在通过适用于当前目标应用和环境的功能&#xff0…

3DMax中场景太大如何优化?

如果你在3dMax中进行大型项目时曾因性能低下或崩溃而感到沮丧,那么你就来对地方了。大型3dMax场景由于其复杂性和文件大小而具有挑战性,但使用正确的优化技术,你可以改进工作流并获得更平滑的结果。在今天的文章中,我们将深入探讨如何在大量3ds Max场景中增强性能。 原因:…

5月4(信息差)

&#x1f384; HDMI ARC国产双精度浮点dsp杜比数码7.1声道解码AC3/dts/AAC环绕声光纤、同轴、USB输入解码板KC33C &#x1f30d; 国铁集团回应高铁票价将上涨 https://finance.eastmoney.com/a/202405043066422773.html ✨ 源代码管理平台GitLab发布人工智能编程助手DuoCha…

安装mysql-8.0.19-winx64.zip步骤

1、cmd以管理员身份打开 2、dos切换到mysql的bin目录下 3、doc执行 mysqld --initialize --console&#xff08;注意这里会生成临时密码&#xff09; 4、继续执行&#xff1a; mysqld --install5、启动服务器&#xff1a; net start mysql 完毕&#xff01;&#xff01; …

Python系列五之正则表达式

正则表达式&#xff0c;Regular Expression&#xff0c;可用于在一个目标字符串里对于指定模式的字符进行查找、替换、分割等操作。 比如&#xff0c;判断某个字符串里是否都是数字&#xff0c;或者是否包含指定字符串&#xff0c;又或者更直接的例子是判断电话号码或者邮箱是…

【前端开发---Vue2】史上最详细的Vue2入门教程,从基础到进阶带你彻底掌握Vue(三)

本篇重点分享常见指令修饰符、v-bind指令用于 class 类名 和 style 行内样式 动态控制、v-model在其他表单元素的使用...... 并结合具体案例来让小伙伴们掌握的更透彻&#xff01;喜欢就先关注一下吧~ 声明&#xff1a;图片资源来自于黑马程序员公开学习资料 本人在学习当中&am…

golang学习笔记(协程的基础知识)

golang的协程 协程是一种轻量级的线程&#xff0c;它可以实现并发执行的并行操作。协程是Go语言中的一个核心特性&#xff0c;它使得程序能够以并发的方式运行&#xff0c;并且非常高效。与传统的线程相比&#xff0c;协程的创建和销毁成本非常低&#xff0c;可以方便地启动大…

浏览器中不能使用ES6的扩展语法...报错

浏览器大多数已经支持ES6&#xff08;ECMAScript 2015&#xff09;的扩展语法&#xff08;...&#xff09;&#xff0c;包括Chrome、Firefox、Safari和Edge等。然而&#xff0c;如果你在某些浏览器中遇到无法使用扩展语法的问题&#xff0c;可能是由以下原因导致的&#xff1a;…

致敬每一位辛勤耕耘的程序员

致敬每一位辛勤耕耘的程序员 在数字世界的浩瀚海洋中&#xff0c;有这样一群人&#xff0c;他们默默无闻&#xff0c;却用一行行代码编织出神奇的世界&#xff1b;他们日夜兼程&#xff0c;只为让科技的力量更贴近人们的生活。他们&#xff0c;就是我们最可爱的程序员们。今天…

神经网络之防止过拟合

今天我们来看一下神经网络中防止模型过拟合的方法 在机器学习和深度学习中&#xff0c;过拟合是指模型在训练数据上表现得非常好&#xff0c;但在新的、未见过的数据上表现不佳的现象。这是因为模型过于复杂&#xff0c;以至于它学习了训练数据中的噪声和细节&#xff0c;而不…

一款开源高性能AI应用框架

前言 LobeChat 是一个基于 Next.js 框架构建的 AI 会话应用&#xff0c;旨在提供一个 AI 生产力平台&#xff0c;使用户能够与 AI 进行自然语言交互。 LobeChat应用架构 LobeChat 的整体架构由前端、EdgeRuntime API、Agents 市场、插件市场和独立插件组成。这些组件相互协作&a…

代码随想录算法训练营二刷第14天 | 二叉搜索树完结

235.二叉搜索树的最近公共祖先 要记住 二叉搜索树自带搜索方向 这个特性&#xff01;要找 二叉搜索树的最近公共祖先&#xff0c;那么相当于就是二叉搜索树的搜索&#xff0c;要想到“二叉搜索树的搜索”的迭代法其实是利用了二叉树的特性的&#xff0c;那么&#xff0c;本题也…

P8799 [蓝桥杯 2022 国 B] 齿轮

P8799 [蓝桥杯 2022 国 B] 齿轮 分析 最右边的齿轮的转速是最左边齿轮的q倍 最右边的齿轮的半径是最左边齿轮的q倍 题意即为&#xff1a;查询数组中是否存在两个数&#xff0c;其中一个是另一个的q倍 题目范围&#xff1a;查询次数q:2*10^5&#xff0c;数组范围2*10^5&…

XXL-JOB定时任务

1. xxl-job初识 1.1 xxl-job介绍 xxl-job 是大众点评大佬徐雪里开源的一款分布式任务调度框架&#xff0c;具有简单易用、轻量级、可扩展的特点。相比于Spring Task, Quartz&#xff0c;xxl-job有记录执行日志和运行大盘&#xff0c;方便开发人员和运维人员更好的管理任务。 …

Docker no space left on device

Docker no space left on device 原因&#xff1a;空间不足 解决方式&#xff1a;清理空间,删除未被使用的容器&#xff0c;镜像等资源 docker system prune -a Docker是一种相对使用较简单的容器&#xff0c;我们可以通过以下几种方式获取信息&#xff1a; 1、通过docker ru…

如何进行面向对象分析、面向对象设计和面向对象编程

目录 1.引言 2.案例介绍和难点剖析 3.如何进行面向对象分析 4.如何进行面向对象设计 5.如何进行面向对象编程 6.总结 1.引言 面向对象分析(OOA)、面向对象设计(00D)和面向对象编程(OOP)是面向对象开发的3个主要环节。 在以往的工作中&#xff0c;作者发现&#xff0c;很多…

【JavaEE 初阶(一)】初识线程

❣博主主页: 33的博客❣ ▶️文章专栏分类:JavaEE◀️ &#x1f69a;我的代码仓库: 33的代码仓库&#x1f69a; &#x1faf5;&#x1faf5;&#x1faf5;关注我带你了解更多线程知识 目录 1.前言2.进程3.线程4.线程和进程的区别5.Thread创建线程5.1继承Thread创建线程5.2实现R…

为什么Java的id用long但是做JavaWeb项目时id用String而不是long

为什么Java中id用long? 1.从整数的范围来看 int: 最小值&#xff1a;-2,147,483,648 最大值&#xff1a;2,147,483,647 在Java中&#xff0c;int类型是32位有符号整数&#xff0c;范围是从 -2^31 到 2^31-1。 long: 最小值&#xff1a;-9,…

专项技能训练五《云计算网络技术与应用》实训7-2:使用OpenDaylight界面下发流表

文章目录 使用OpenDaylight的web界面下发流表1. 根据前面实训教程&#xff0c;启动并登录到opendaylight的web控制页面&#xff0c;并使用mininet2. 单击左侧的“Nodes”查看节点信息。其中尤其需要关注Node Id&#xff0c;下发流表的时候会用到Node Id&#xff0c;如下图所示。…

怎么查看java 安装路径

查看Java安装路径的方法因操作系统而异。以下是几种常见操作系统中查看Java安装路径的方法&#xff1a; Windows系统&#xff1a; 使用命令提示符&#xff1a; 打开命令提示符&#xff08;cmd&#xff09;。输入以下命令&#xff1a;for %i in (java.exe) do echo %~$PATH:i。…