R语言混合线性模型中BLUE值和BLUP值实例分析
发表于:2024-10-02 作者:千家信息网编辑
千家信息网最后更新 2024年10月02日,本文小编为大家详细介绍"R语言混合线性模型中BLUE值和BLUP值实例分析",内容详细,步骤清晰,细节处理妥当,希望这篇"R语言混合线性模型中BLUE值和BLUP值实例分析"文章能帮助大家解决疑惑,下
千家信息网最后更新 2024年10月02日R语言混合线性模型中BLUE值和BLUP值实例分析
本文小编为大家详细介绍"R语言混合线性模型中BLUE值和BLUP值实例分析",内容详细,步骤清晰,细节处理妥当,希望这篇"R语言混合线性模型中BLUE值和BLUP值实例分析"文章能帮助大家解决疑惑,下面跟着小编的思路慢慢深入,一起来学习新知识吧。
混合线性模型中BLUE值 VS BLUP值
最佳线性无偏预测(best linear unbiased prediction, 简称BLUP),又音译为"布拉普"[1],是统计学上用于线性混合模型对随机效应进行预测的一种方法。最佳线性无偏预测由C.R. Henderson提出。随机效应的最佳线性无偏预测(BLUP)等同于固定效应的最佳线性无偏估计(best linear unbiased estimates, BLUE)(参见高斯-马尔可夫定理)。因为对固定效应使用估计一词,而对随机效应使用预测,这两个术语基本是等同的。BLUP被大量使用于动物育种。--wiki
BLUP值,相当于是对混合线性模型中随机因子的预测;
BLUE值,相当于是对混合线性模型中固定因子的估算
predict means:预测均值,固定因子和随机因子都可以预测均值,它的尺度和表型值尺度一致
将处理作为固定因子
将处理作为固定因子setwd("D:\\02 ASReml\\blue VS blup")library(asreml)library(tidyverse)dat <- read.csv("MaizeRILs.csv",head=T)for (i in 1:4) dat[,i] <- as.factor(dat[,i])as1 <- asreml(height ~ location/rep + location*RIL,data=dat)ASReml: Tue May 08 11:07:55 2018 LogLik S2 DF wall cpu -723.8797 64.8862 244 11:07:55 0.1 -723.8797 64.8862 244 11:07:55 0.0Finished on: Tue May 08 11:07:55 2018LogLikelihood Converged #计算品种的BLUE值ablue <- coef(as1)$fixedblue1 <- ablue[grep("^RIL_RIL*",rownames(ablue)),] %>% as.data.frame()head(blue1)#计算品种的预测均值(predict means)pv1 <- predict(as1,"RIL")$predictions$pvalsASReml: Tue May 08 11:13:33 2018 LogLik S2 DF wall cpu -723.8797 64.8862 244 11:13:33 0.0 -723.8797 64.8862 244 11:13:33 0.0Finished on: Tue May 08 11:13:33 2018LogLikelihood Converged head(pv1) # 类似SAS中的lsmeans#运行模型:因素作为随机因子as2 <- asreml(height ~ 1,random = ~location/rep + location*RIL,data=dat)ASReml: Tue May 08 11:13:34 2018 LogLik S2 DF wall cpu -1646.5302 233.5135 495 11:13:34 0.0 -1569.0397 137.9186 495 11:13:34 0.0 -1507.3257 94.6888 495 11:13:34 0.0 -1471.3354 74.5149 495 11:13:34 0.0 -1462.9209 67.6142 495 11:13:34 0.0 -1461.7649 65.3553 495 11:13:34 0.0 -1461.7228 64.9069 495 11:13:34 0.0 -1461.7228 64.8863 495 11:13:34 0.0 -1461.7228 64.8862 495 11:13:34 0.0Finished on: Tue May 08 11:13:34 2018LogLikelihood Converged blup <- coef(as2)$randomblup2 <- blup[grep("^RIL_RIL-*",rownames(blup)),] %>% as.data.frame()head(blup2)#预测均值pv2 <- predict(as2,"RIL")$predictions$pvalsASReml: Tue May 08 11:13:34 2018 LogLik S2 DF wall cpu -1461.7228 64.8862 495 11:13:34 0.0 -1461.7228 64.8862 495 11:13:34 0.0 -1461.7228 64.8862 495 11:13:34 0.0 -1461.7228 64.8862 495 11:13:34 0.0Finished on: Tue May 08 11:13:34 2018LogLikelihood Converged head(pv2)#计算遗传力summary(as2)$varcompstr(dat)'data.frame': 496 obs. of 9 variables: $ location: Factor w/ 4 levels "ARC","CLY","PPAC",..: 1 1 2 2 3 3 4 4 1 1 ... $ rep : Factor w/ 2 levels "1","2": 1 2 1 2 1 2 1 2 1 2 ... $ block : Factor w/ 8 levels "1","2","3","4",..: 4 6 5 4 8 5 1 4 1 2 ... $ plot : Factor w/ 122 levels "1","2","3","4",..: 28 47 36 92 64 40 7 27 6 9 ... $ RIL : Factor w/ 62 levels "RIL-1","RIL-11",..: 1 1 1 1 1 1 1 1 2 2 ... $ pollen : int 73 74 71 73 97 95 72 72 69 69 ... $ silking : int 77 79 74 77 101 100 78 78 71 72 ... $ ASI : int 4 5 3 4 4 5 6 6 2 3 ... $ height : num 182 169 213 203 156 ...VSNR::pin(as2,h3 ~ V3/(V3 + V4/4 + V5/(2*4)))将数据保存到excel中library(openxlsx)write.xlsx(blue1,"blue.xlsx")write.xlsx(blup2,"blup.xlsx")write.xlsx(pv1,"pm1.xlsx")write.xlsx(pv2,"pm2.xlsx")结果解析RIL是基因型pm2-random是RIL作为随机因子的预测均值pm1-fixed是RIL作为固定因子时的预测均值blue是RIL作为固定因子的BLUE值blup是RIL作为随机因子的BLUP值pm2-blup 是随机因子的预测均值 减去 随机因子的BLUP值,可以看到得到的是一个常数(均值)pm1-mu-random 是固定因子的预测均值 减去 固定依着你的BLUE值, 可以看到不是一个常数blup/blue_effect=heritibility 是BLUP值 除以 BLUE效应值,得到的是遗传力常数备注:blue_effect是用固定因子的预测均值 减去 整体均值
读到这里,这篇"R语言混合线性模型中BLUE值和BLUP值实例分析"文章已经介绍完毕,想要掌握这篇文章的知识点还需要大家自己动手实践使用过才能领会,如果想了解更多相关内容的文章,欢迎关注行业资讯频道。
因子
线性
均值
模型
混合
效应
无偏
实例
实例分析
语言
分析
常数
文章
处理
内容
品种
尺度
遗传
妥当
一致
数据库的安全要保护哪些东西
数据库安全各自的含义是什么
生产安全数据库录入
数据库的安全性及管理
数据库安全策略包含哪些
海淀数据库安全审计系统
建立农村房屋安全信息数据库
易用的数据库客户端支持安全管理
连接数据库失败ssl安全错误
数据库的锁怎样保障安全
高校网络安全配置命令
幼儿园小班网络安全教育图文
网络技术具体内容
网络安全民法
nba数据库进攻效率
linux系统当服务器
中山专业软件开发批发价格
数据库竞争态势
戴尔服务器双路主板和cpu
怎么查看服务器内存和cpu
网络技术对社会的贡献
沧州汇德网络技术
大专学历的软件开发人员
佳都科技是互联网公司吗
fm2017数据库
摸排校园网络安全情况说明
厦门移动医疗健康软件开发
翔优网络安全绘画
灵娱网络技术有限公司招聘
网络安全管理员报考网站
联想服务器管理端口
乐天华网络技术工作怎么样
墨小密南京网络技术
数据库技术与应用期末总结
南宁聚通付网络技术有限公司
河南金蝶财务软件开发
2022年四月份网络安全宣传
一本通 网络技术
网络安全法'
招聘软件开发乌鲁木齐