基于R语言绘制GGE双标图2

参考资料:

严威凯等: 双标图分析在农作物品种多点试验中的应用【作物学报】

https://cran.r-project.org/web/packages/GGEBiplots/GGEBiplots.pdf

1、如何判断双标图是否充分体现数据中的规律

        在对双标图的解释中,有一个隐含的假设,就是所讨论的 2-D双标图可以充分近似所代表的两
向表数据。

        如何知道上面这个假设是否满足,首先要看双标图的拟合度(Goodness of Fit), 即前两个主成分(PC1 和 PC2)所解释的两向表中总变异的百分数。如果拟合度高, 则双标图能较好地近似实际数据。

        假定待分析数据中有 g 个品种和 e 个环境(地点)。那么百分之百表现这套数据所需要的最大主成分数目是 K =min (g–1, e)。如果数据中品种间和环境间没有任何关系或规律, 每个主成分所解释的变异应当都是总变异的1/K。如果数据中存在一些关系或规律, 那么有些主成分所解释的变异的比例就会大于1/K,而其他的主成分所解释的比例则小于1/K

        按此推理, Yan 和 Tinker提出了“信息比”(IR)的概念, 它是各主成分所解释总变异的份额乘以 K。这样, 对于一个主成分, IR >1 表示含有规律性信息,IR = 1 表示含有独立性信息,IR <1 则表示不含任何规律或重要信息。对于一套数据, 有几个主成分具有 IR ≥1, 就需要几个主成分来充分近似之。如果只需要 1 个主成分, 那么双标图分析就是多余的。如果需要 2 个主成分, 则 2-D 双标图足以表现数据中规律。如果需要 3 个主成分, 则 2-D 双标图不足以表现数据中所有规律,但3-D 双标图则可。如果需要 3 个以上主成分, 说明数据中规律复杂但没有突出规律。

        以GGEBiplots包自带的Ontario案例进行分析,查看本例的主成分信息比情况:

# 查看主成分的特征值、解释总变异的百分数及信息比
K=min(c(nrow(Ontario)-1,ncol(Ontario)))
eigval<-GGE1$eigenvalues
expl<-GGE1$varexpl
IR<-expl*K/100
tab<-data.frame(特征值=eigval,解释总变异百分比=expl,信息比=IR)
tab

        结果显示:只有前 2 个主成分满足 IR >1, 可以认为双标图可以充分体现数据中的规律。

        经验表明,多数情况下,2-DGGE双标图足以充分近似一年内多点试验的数据。但在数据较大(同时有很多品种和很多环境)或涉及到广泛的地理分布,或把不同年份的数据放在一起分析时,情况可能不同。2-D双标图不能充分近似数据时,可以考虑采取以下方法:
        第1种方法是使用3-D双标图。
        第2种方法是把试验点分组并分别进行研究。按照2-D双标图表现出的环境间关系把环境分组,然后按组分别研究,是解决数据较大、关系较复杂的有效途径。
        第3种方法是去掉数据中的低产品种或产量不稳定品种。2-D双标图不能充分表现数据中的规律多是因为数据中的GE关系复杂而得不到充分表现;后者又可能是因为品种主效益较强所致。在这个前提下,去掉一些在所有或多数环境下都低产的品种,可以减少双标图中G的分量,使较多的GE得以在双标图上表现,从而取得对环境更好的区分。

2、如何检验双标图上关系的显著性

        由于双标图本身没有一个判断显著性的尺度,在按照双标图做出重要决定之前需要进行统计检验。因此,原则上应当把双标图当作一个建立假说的工具而不是作出最终结论的工具。具体如何检验,可参考:严威凯等: 双标图分析在农作物品种多点试验中的应用【作物学报】。

3、双标图的定标

        未定标的GGE双标图是不经过数据定标,直接将环境中心化的数据进行特征值分解而建立的。其特点是,在采用聚焦环境的特征值分配(SVP="column")时,图中各环境向量的长度与各环境内品种均值间的标准差(SD,等于环境内表型方差的平方根)呈近似直线关系。近似的程度与双标图的拟合度有关。某试验点内SD的大小可以作为衡量该试验点对品种的表型区分能力的指标。

        用未定标GGE双标图对试验点进行有效评价是有条件的,即各试验点的试验误差相同或相近(称为误差同质性原则);而这个条件在实际中不一定能满足。

        数据定标是把同一环境内各品种的数值除以(或乘以)某个能够表征该环境的量,以便把各环境放在同一尺度下比较。GGEBiplots包中只提供了SD标定,故这里只讨论SD标定。用环境内品种均值间的标准差(SD)对各环境进行定标,即以环境SD去除两向表中相应环境内各品种的数值,可以消除各环境在SD上的差异。在SD-定标的GGE双标图上,在拟合度较高的前提下,各个环境之向量的长度应当是相同或相近的。这一特性有两个用途。第一,由于各环境的向量长度相近,此图更便于直观分析各环境在品种排序上的相似性。第二,可据以判断双标图是否充分表现数据。如果双标图上各环境的向量长度显著不同,则说明该双标图不能充分体现数据中环境间的关系

        与这些优点相伴而来的缺点是, SD定标双标图上环境向量的长度不再是其对品种区分能力的度量

# 查看数据中各试验点的标准差
apply(Ontario,2,sd)
# 对数据进行定标
GGE2<-GGEModel(Ontario,scaling="sd")
# 展示“区分力和代表性”功能图
DiscRep(GGE2)

        SD-定标的GGE 双标图与多元分析中常用的主成分分析模型相一致, 因而在多元数据(如品种-性状两向表, 品种-遗传标记两向表, 等等)分析中有广泛应用。在需要把各环境看得同等重要时, SD-定标的GGE 双标图应当是第一选择

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

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

相关文章

COFFEE AI PARTNER -- 神奇的AI工具,相当我雇佣了一个AI员工,淘汰你的是会使用AI的人

COFFEE AI PARTNER介绍 COFFEE AI PARTNER是由 AI JAVA开发的一款生成式人工智能工具&#xff08;又名AI助手&#xff09;&#xff0c;尝试一下。 首先域名似乎正在备案中&#xff0c;企业邮箱似乎正在采购&#xff0c;目前服务地址是&#xff1a;COFFEE AI PARTNER-官网 官网…

git安装图文

1.下载 通过百度网盘分享的文件&#xff1a;git安装图文 链接&#xff1a;https://pan.baidu.com/s/17ZMiWUIULtrGGba5n-WLeA 提取码&#xff1a;anjm --来自百度网盘超级会员V3的分享 2.安装

【C语言篇】猜数字游戏(赋源码)

文章目录 猜数字游戏前言随机数生成randsrandtime设置随机数生成范围 猜数字游戏的实现 猜数字游戏 前言 在前两篇博客对于分支和循环语句进行了详细的介绍&#xff1a; 分支语句详解 循环语句详解 我们就可以写一写稍微有趣的代码了&#xff0c;比如&#xff1a; 写一个…

二分查找法

&#x1f49d;&#x1f49d;&#x1f49d;欢迎来到我的博客&#xff0c;很高兴能够在这里和您见面&#xff01;希望您在这里可以感受到一份轻松愉快的氛围&#xff0c;不仅可以获得有趣的内容和知识&#xff0c;也可以畅所欲言、分享您的想法和见解。 非常期待和您一起在这个小…

SQL各种注入详解加案例--持续更新

sql注入 联合查询注入案例手工注入判断是否有SQL注入漏洞 sqlmap工具注入 报错注入常用的函数updatexml()函数案例 floor()涉及的函数实现手工注入sqlmap工具注入 盲注布尔盲注案例手工注入脚本sqlmap自动化工具 时间盲注 post注入GET传参和POST传参案例手工注入sqlmap工具 二次…

使用 Python 制作一个属于自己的 AI 搜索引擎

1. 使用到技术 OpenAI KEYSerper KEYBing Search 2. 原理解析 使用Google和Bing的搜搜结果交由OpenAI处理并给出回答。 3. 代码实现 import requests from lxml import etree import os from openai import OpenAI# 从环境变量中加载 API 密钥 os.environ["OPENAI_AP…

MySQL:索引(Index)语句

索引的限制 每个表最多可以有 16 个索引&#xff08;InnoDB 表的限制&#xff09;。 单个索引最多可以包含 16 列。 索引列的最大长度为 767 字节&#xff08;对于 CHAR, VARCHAR, 和 BINARY 类型&#xff09;&#xff0c;3072 字节&#xff08;对于 BLOB 类型&#xff09;。…

浅谈取样器插件之bzm - Free-Form Arrivals Thread Group

浅谈取样器插件之bzm - Free-Form Arrivals Thread Group bzm - Free-Form Arrivals (Ultimate Thread Group) 是一个高级且灵活的线程组插件&#xff0c;专为Apache JMeter设计。它扩展了JMeter的标准线程组功能&#xff0c;允许用户以自由形式定义线程&#xff08;用户&…

SSM项目学习:用xml配置文件或注解开发实现控制反转和依赖注入

什么是SSM SSMSpring(Spring Framework)Spring MVC mybatis Spring Framework系统架构 Spring Framework学习线路 IoC(Inversion of Control)和DI(Dependency Injection) 他们解决的问题&#xff1a;代码耦合度高的问题&#xff0c;需要类自己new对象&#xff0c;修改部分代…

03、DQL(数据查询语句)

目录 1、编写顺序 2、基本查询 3、条件查询 4、聚合函数 5、分组查询 6、排序查询 7、分页查询 8、执行顺序 1、编写顺序 SELECT 字段列表 FROM 表名列表 WHERE 条件列表 GROUP BY 分组字段列表 HAVING 分组后条件列表 ORDER BY 排序字段列表 LIMIT 分页参数2、基本查…

简单的docker学习 第11章 镜像中心

第11章 镜像中心 Docker Hub 与阿里云都是 Docker 的公网镜像中心&#xff0c;用户可以将自己的镜像 push 到公网镜像中心中自己的镜像仓库&#xff0c;并可将仓库设置为私有库&#xff0c;使他人无法看到&#xff0c;更无法 pull&#xff0c;以保证镜像的安全性。不过&#x…

【LeetCode刷题笔记】LCR.27 回文链表

创作不易&#xff0c;本篇文章如果帮助到了你&#xff0c;还请点赞 关注支持一下♡>&#x16966;<)!! 主页专栏有更多知识&#xff0c;如有疑问欢迎大家指正讨论&#xff0c;共同进步&#xff01; 更多算法知识专栏&#xff1a;算法分析&#x1f525; 给大家跳段街舞感谢…

为什么康耐视visionpro的C#二次开发调用的recorddisplay控件偶尔会显示白色的,偶尔又正常了?

recorddisplay控件正常显示 异常显示 原因分析&#xff1a; 没有完全加载recorddisplay控件&#xff0c;有可能是有bug没有完全加载&#xff0c;打断点调试控件是否完全加载。

EMQX服务器安装MQTT测试

cd /usr/local/develop wget https://www.emqx.com/en/downloads/broker/5.7.1/emqx-5.7.1-el7-amd64.tar.gz mkdir -p emqx && tar -zxvf emqx-5.7.1-el7-amd64.tar.gz -C emqx ./emqx/bin/emqx start 重启 ./emqx/bin/emqx restart http://10.8.0.1:18083/ 账号ad…

【Kubernetes】应用的部署(一):金丝雀部署

应用的部署&#xff08;一&#xff09;&#xff1a;金丝雀部署 在项目迭代开发过程中&#xff0c;经常需要对应用进行上线部署。上线部署策略主要有 3 种&#xff1a;金丝雀部署、蓝绿部署 和 滚动部署。 金丝雀部署 也被叫作 灰度部署。金丝雀部署过程&#xff1a;先让一部分…

letcode 分类练习 哈希表 242.有效的字母异位词 349. 两个数组的交集 202. 快乐数 1. 两数之和

letcode 分类练习 哈希表 242.有效的字母异位词 349. 两个数组的交集 202. 快乐数 1. 两数之和 242.有效的字母异位词349.两个数组的交集202.快乐数1.两数之和 242.有效的字母异位词 分别定义两个字母哈希表就可以了 class Solution { public:bool isAnagram(string s, strin…

搭建pxe网络安装环境

实验目的&#xff1a; 搭建pxe网络安装环境实现服务器自动部署 实验原理&#xff1a; PXE 网络安装环境实现服务器自动部署的实验原理为&#xff1a; 待安装的服务器&#xff08;PXE 客户端&#xff09;开机时&#xff0c;BIOS 设置从网络启动&#xff0c;向网络发送请求。…

科普文:JUC系列之ForkJoinPool源码解读ForkJoinWorkerThread

科普文&#xff1a;JUC系列之ForkJoinPool基本使用及原理解读-CSDN博客 科普文&#xff1a;JUC系列之ForkJoinPool源码解读概叙-CSDN博客 科普文&#xff1a;JUC系列之ForkJoinPool源码解读WorkQueue-CSDN博客 科普文&#xff1a;JUC系列之ForkJoinPool源码解读ForkJoinTask…

【第13章】Spring Cloud之Gateway全局异常处理

文章目录 前言一、异常处理1. 响应实体类2. 异常处理类 二、单元测试1. 无可用路由2. 服务不可用 总结 前言 网关作为我们对外服务的入口起着至关重要的作用&#xff0c;我们必须保证网关服务的稳定性&#xff0c;下面来为网关服务增加异常处理机制。 一、异常处理 1. 响应实…

K个一组翻转链表(LeetCode)

题目 给你链表的头节点 &#xff0c;每 个节点一组进行翻转&#xff0c;请你返回修改后的链表。 是一个正整数&#xff0c;它的值小于或等于链表的长度。如果节点总数不是 的整数倍&#xff0c;那么请将最后剩余的节点保持原有顺序。 你不能只是单纯的改变节点内部的值&…