首页
学习
活动
专区
圈层
工具
发布

多个性状的BLUP计算综合育种值

最近有老师问如何计算多个育种值的综合育种值,因为选择育种时,需要考虑多个性状,而这些性状单位不一样,怎么合并计算?...这里面其实就是包括两部分内容: 1,多个性状的育种值进行标准化 2,多个性状的育种值赋予权重 为何要进行标准化? 所谓标准化,就是将数据变为平均数为0,方差为1的数据集。...将所有性状都标准化,就意味着这些性状的尺度都是一样的,是一个总体,可以进行加减乘除的操作。然后可以对不同的性质设置不同的权重,才可以进行综合育种值的计算。...最后将权重对应的性状相乘,相加变为综合育种值。...) d2$y3 = (d2$y3-mean(d2$y3))/(sd(d2$y3)) d2$index = d2$y1*0.3 + d2$y2*0.3 + d2$y3*0.4 head(d2) 两者结果是一致的

1.5K10
  • 您找到你想要的搜索结果了吗?
    是的
    没有找到

    如何合并多个体细胞突变检测工具的输出结果

    简介 在肿瘤基因组数据分析方法中,很多文章会使用多个软件 call 突变,然后进行合并,常用见的方法就是多个软件得到的结果中,一个突变在两个或两个以上的软件检测到即保留,认为其可信度较高。...主要是针对体细胞突变检测工具 Mutect2 和 Strelka2 的结果。 vcf 文件格式 既然是要对 vcf 文件进行合并,就要深入了解 vcf 文件格式。...表示 REF T 参考基因组中所对应的碱基 ALT G 样品中该位置突变后的碱基 QUAL ....其体细胞突变检测的结果分为 SNV 和 INDELs 两个 vcf 文件,需要先进行合并: zcat 6.strelka/case1_biorep_B/results/variants/somatic.snvs.vcf.gz...可以看到合并后的 vcf 文件 header 部分有多出来 bcftools 的处理信息。 这种方法比较简单粗暴,实际上也有其他工具可以实现,欢迎大家在评论区提出自己的方法。

    30900

    dotnet 修复 ILLinkTasksAssembly 特性的值的计算结果无效

    提示 元素 UsingTask 中“AssemblyFile”特性的值“$(ILLinkTasksAssembly)”的计算结果“”无效。...如果发现自己的设备上不存在 Microsoft.NET.ILLink.Tasks 这个文件夹,那么请将 dotnet sdk 卸载重新安装,或者安装更新版本的 sdk 然后查看自己的环境变量,是否有设置特定版本的...dotnet sdk 如果有设置,就删除此项值或者修改为更新版本 创建空白项目,找找自己本地的 ILLinkTasksAssembly 定义的做法是在空项目的 csproj 里面添加如下代码,用来输出...本文的方法能修复的是在构建和加载项目提示如下内容 error : 元素 中“AssemblyFile”特性的值“$(ILLinkTasksAssembly)”的计算结果“”无效...TaskName="ILLink" AssemblyFile="$(ILLinkTasksAssembly)" /> 上面提示内容就是 $(ILLinkTasksAssembly) 没有定义,拿到空值。

    1.3K20

    js中带有参数的函数作为值传入后调用问题

    ❝小闫语录:你可以菜,但是就这么菜下去是不是有点过分了 ❞ 每天不是在写 bug,就是在解 bug 的路上~更多精彩文章请关注公众号『Pythonnote』或者『全栈技术精选』 1.无参数函数作为参数传入调用...当根据实际情况,函数需要作为参数传入时,一般采用如下方式直接调用即可: function fuc1() { console.log(1); } function fuc2(a) { a();...} fuc2(fuc1); // 1 2.有参数函数作为参数传入调用 一般的函数都有参数,那么这种情况如何传参呢?...console.log(param); } function fuc2(a, b) { a(b); } fuc2(fuc1, "欢迎关注微信公众号:全栈技术精选"); 3.有参数函数作为事件方法 现在要将传入的函数作为点击事件的处理程序...❝因为在你写 fuc1("我是小闫同学啊") 时,默认就调用了此函数,都不需要点击。 ❞ 如何才能达到在点击时才弹出窗口呢?

    11.7K40

    【Kotlin 协程】Flow 异步流 ① ( 以异步返回返回多个返回值 | 同步调用返回多个值的弊端 | 尝试在 sequence 中调用挂起函数返回多个返回值 | 协程中调用挂起函数返回集合 )

    文章目录 一、以异步返回返回多个返回值 二、同步调用返回多个值的弊端 三、尝试在 sequence 中调用挂起函数返回多个返回值 四、协程中调用挂起函数返回集合 一、以异步返回返回多个返回值 ----...| 协程的 suspend 挂起函数 ) 博客 ; 如果要 以异步的方式 返回多个元素的返回值 , 可以使用如下方案 : 集合 序列 Suspend 挂起函数 Flow 异步流 二、同步调用返回多个值的弊端...System.out: 4 2022-12-22 12:33:04.703 15427-15427/kim.hsl.coroutine I/System.out: 5 三、尝试在 sequence 中调用挂起函数返回多个返回值...---- 如果要 以异步方式 返回多个返回值 , 可以在协程中调用挂起函数返回集合 , 但是该方案只能一次性返回多个返回值 , 不能持续不断的 先后 返回 多个 返回值 ; 代码示例 : package...super.onCreate(savedInstanceState) setContentView(R.layout.activity_main) // 携程中调用挂起函数返回多个值

    13K30

    aardio的whttp库调用post()后如何获取header中的cookie值

    目前的whttp库调用get和post后无法通过readHeader()函数读取返回的header。...因为readHeader函数必须在请求完成之前调用才能获取到header,而一鹤写的库里面,只有请求的method=”head”时才调用这个函数, 其他如post、get方法都不会调用。...控制这个逻辑的代码在whttp库的down函数里面,大概573行: if( method == “HEAD” || noReceiveData ){ this.readHeader(); this.endRequest...(); return true; } 只要把    this.readHeader(); 这行代码移动到这个判断语句外面, 就可以在post()之后再调用readHeader来获取返回的http头了...其实whttp是可以自动保存cookie的,那为什么我非要把它读出来呢? 是为了在多线程中共用cookie,才必须把这个header读出来。 本人和一鹤沟通, 希望把这个库这样改一下, 被拒绝。

    1.6K40

    明明结果是对的,为什么被合并查询后得到的结果却出错了?| Power Query躲坑

    最近,有位朋友在一个实际工作问题中,在表2使用合并查询从表1的结果中匹配最高(阶段)项,眼看着表1的结果是对的,但表2里却得到了错误的返回结果,具体情况如图所示: 为什么会这样?...我们先来看表1的处理情况。 为了合并查询得到最高阶段项,对表1进行降序排序: 然后通过删除重复项保留最高阶段数据: 从表1的结果来看,的确保留了最高阶段的数据。...然后,在表2里使用合并查询获取表1中的结果并展开: 咦!!! 表1的处理结果明明是阶段4(报价),为什么合并查询得到的结果却是阶段2(售前)? 这难道是Power Query的Bug吗?...实际上,经过这么多年的发展,Power Query已经相当成熟了,合并查询这种基本的功能,不大可能存在这样的Bug或低级错误的。...这里的问题根源其实是表1的处理问题,我以往发布的多篇文章案例中,在涉及Power Query中使用排序的问题时会强调,Power Query的排序需要增加添加索引或Table.Buffer的步骤,使排序的结果真正

    4.2K20

    Python 源代码里的算法——如何合并多个有序列表并使得结果依然有序?

    摄影:产品经理 朝闻道,晚上喝酒 去年的一篇文章《一日一技:在 Python 里面如何合并多个有序列表并使得结果依然有序?》,我很自不量力地提到了“多个有序列表”。...但实际上,那篇文章仅仅是合并两个有序列表而已。真正要合并多个有序列表并使结果依然有序,会难得多。...与 E 合并得到最终结果。...先把 A、B 列表的元素合并,得到 X;再把 C、D列表的元素合并得到 Y、然后 XY 合并得到 Z;最后把 Z 与 E 合并得到最终结果。...接下来,第359行,把堆顶列表的下标为0的项替换为新的值(原来所在列表的第二个值)。 第360行,把当前堆顶的列表替换为新的列表。

    2.5K10

    这种作者只提供了每组多个重复样本合并后数据的怎么办

    数据解释和验证:最后,对分析结果进行解释,并可能需要通过实验(如qPCR)来验证关键基因的表达变化。...早期分享的关于mfuzz 分析的帖子: Mfuzz做转录变化的时间趋势分析后对每个趋势分组挑一个代表性基因 使用Mfuzz包做时间序列分析 关闭利用Mfuzz包对转录变化的时间趋势进行分析(学徒作业参考...fread(fs[1],data.table = F) head(tmp) gid <- fread(fs[1],data.table = F)[,2] head(gid) # 批量读取并按照列合并...rownames(rpkm) colnames(rpkm_exp) <- colnames(rpkm) head(rpkm_exp) 表达矩阵如下: 学徒作业 拿到 GSE65031 数据集预处理后的表达矩阵...,进行 Mfuzz 分析,并与 文章中的结果进行比较,文献《Whole transcriptome analysis reveals a role for OGG1-initiated DNA repair

    36911

    【Leetcode -1171.从链表中删去总和值为零的连续节点 -1669.合并两个链表】

    Leetcode -1171.从链表中删去总和值为零的连续节点 题目:给你一个链表的头节点 head,请你编写代码,反复删去链表中由 总和 值为 0 的连续节点组成的序列,直到不存在这样的序列为止。...删除完毕后,请你返回最终结果链表的头节点。 你可以返回任何满足题目要求的答案。 (注意,下面示例中的所有序列,都是对 ListNode 对象序列化的表示。)...对于链表中的每个节点,节点的值: - 1000 <= node.val <= 1000....//prev迭代 prev = prev->next; } return dummy->next; } Leetcode -1669.合并两个链表...请你将 list1 中下标从 a 到 b 的全部节点都删除,并将list2 接在被删除节点的位置。 请你返回结果链表的头指针。

    61810

    ​cytoscape的十大插件之五--Centiscape(计算多个中心值)

    下面是cytoscape讲师的笔记 一、初识 据文献介绍,Centiscape是目前唯一能同时计算多个中心值的Cytoscape插件。...Betweenness (Cspb(v),中间度) 定义: 代表节点中心值的指标,与 Stress相似,但含有更多的意义 首先计算网络中任意两点的最短路径数目总和,然后计算其中经过v点的最短路径的数量与总和的比...Closeness (Cclo(v),紧密度) 定义: 代表节点中心值的指标 若计算节点v的离心率,先求出其与网络任意节点的最短距离,后求其总和,最后求倒数 作用: 一般与网络的平均紧密度作比较 紧密度高...,因此是个权重的结果 高质心值的节点,意味着它有很多个邻近点(不单指第一邻近点) 意义: 相对于平均值,高质心值的点更有可能参与协调其他高度连接蛋白的活性,共同致力于某一特定细胞活性的调控(如细胞黏附,...结果 具体数值结果在网络下方呈现 网络右边还有具体调控参数 右方是默认求出各参数平均值,≥阈值的节点会在网络下方的表格出现 可划动调控阈值范围,也可以选择≤阈值的节点,点解 update 后,网络下方结果表格会发生变化

    8.9K63

    【实战】将多个不规则多级表头的工作表合并为一个规范的一维表数据结果表

    最近在项目里,有个临时的小需求,需要将一些行列交叉结构的表格进行汇总合并,转换成规范的一维表数据结构进行后续的分析使用。...从一开始想到的使用VBA拼接字符串方式,完成PowerQuery的M语言查询字符串,然后转换成使用插件方式来实现相同功能更顺手,最后发现,在当前工作薄里使用PowerQuery来获取当前工作薄的其他工作表内容...,也是可行的,并且不需要转换智能表就可以把数据抽取至PowerQuery内。...再最后,发现PowerQuery直接就支持了这种多工作表合并,只要自定义函数时,定义的参数合适,直接使用自定义函数返回一个表结果,就可以展开后得到多行记录的纵向合并(类似原生PowerQuery在处理同一文件夹的多个文件纵向合并的效果...整个实现的过程,也并非一步到位,借着在知识星球里发表,经过各星友一起讨论启发,逐渐完善起来最终的结果。探索是曲折的,但众人一起合力时,就会有出乎意料的精彩结果出来。

    3.6K20

    Excel VBA解读(140): 从调用单元格中获取先前计算的值

    学习Excel技术,关注微信公众号: excelperfect 如果有一个依赖于一些计算慢的资源的用户定义函数,可能希望该用户定义函数在大多数情况下只返回其占用的单元格中最后一次计算得到的值,并且只偶尔使用计算慢的资源...假设要给用户定义函数传递一个计算慢的资源的参数,并让一个开关告诉它何时使用计算慢的资源。...有几种方法可以获得先前为用户定义函数计算的值,它们各有优缺点。...xlerrNA); //coerce the caller ref xloResult.Coerce(xloCaller); return xloResult.Ret(); } 上述代码没有搞懂,有兴趣的朋友研究后可以分享一下...小结 有几种方法可以从VBA用户定义函数的最后一次计算中获取先前的值,但最好的解决方案需要使用C++ XLL。

    10.3K20

    微信公众号网页开发,获取 调用微信公众号接口后 的返回值,PHP

    大家好,又见面了,我是你们的朋友全栈君。...在微信公众号网页开发中,我们通常会用到许多微信团队提供的接口,我们在运行单个php调用接口的时候,返回值直接出现在php的echo中,我们根据这个值来判断代码正确与否,但是当php文件和微信公众号连接起来的时候...,我们无法获取返回值,下面就展示一下如何将返回值的信息写入txt文件中。...;//打开一个文件,如果没有则创建,a+表示文件可读可写 $txt1 = "\r\n";//换行 $txt2 = "返回值信息";//返回值信息 $txt3 = date("Y/m/d G:i:s");...$txt3; fwrite($myfile, $txt); fclose($myfile); } fanhui($txt2); 这样我们就可以通过txt的文件来获取到调用接口的返回值了。

    1.6K20

    【C语言刷题每日一题】——计算11-12+13-14+15 …… + 199 - 1100 的值,打印出结果

    第一种实现方式——每次判断是奇数项还是偶数项对应不同的计算 解题思路 首先需要创建一个浮点型变量来存储表达式的结果,这里暂且称为sum 使用一个for循环来实现除数的递增,从1到100 进入for循环.../ i; } } printf("%lf\n", sum); return 0; } 第二种实现方式——使用一个标记变量来区分每次+还是- 解题思路 首先需要创建一个浮点型变量来存储表达式的结果...,这里暂且称为sum 创建一个标记变量flag,初始为1 使用一个for循环来实现除数的递增,从1到100 进入for循环,每次计算表达式时,额外乘以flag变量,乘完之后改变flag的正负号,进入下一轮...100; i++) { sum = sum + flag * 1.0 / i; flag = -flag; } printf("%lf\n", sum); return 0; } 计算结果...总结 两种解题方式在这个问题的实现上不会有太大差距,但在计算的数据足够庞大时,第二种方式在效率上会更出色一些,因为它比第一种方式每次都少了一条if判断语句

    34410
    领券