育种数据分析之放飞自我分享 http://blog.sciencenet.cn/u/yijiaobai 关注:生物统计,数量遗传,混合线性模型,生物信息,R,Perl,Python,GWAS,GS相关方法,文章及代码

博文

按标题搜索
分类变量logistic回归分析--1
2020-9-23 23:13
1. 二分类logistic回归分析 概念解释 「logistic回归介绍:」 Logistic回归适用于二分类变量(0和1)。模型假设Y服从二项分布,线性模型的拟合形式为: 「什么是优势比(odds radio)」 其中 为Y=1时的优势比。 「什么是对数优势比」 其中 为对数优势比,也称logit。 「什么是 ...
个人分类: 农学统计|15564 次阅读|没有评论
R语言tryCatch使用方法:判断Warning和Error
2020-9-22 19:48
1. 出门掉坑里 因为一些情况,需要判断一行命令运行的状态,然后再做出反应,整体来说: 是否出现warning,出现了怎么处理? 是否出现Error,出现了怎么处理? 没有出现怎么处理? R中判断warning和error状态的函数,我没有找到。一个玩 Java 的同事说, try.....catch 多经典的方法,你怎么不用呢? ...
个人分类: R语言|12874 次阅读|没有评论
GWAS分析中:为何有些SNP效应值大却不显著?
2020-9-16 22:17
之前的GWAS分析中,我查看结果时有过这个疑问,随着不断的学习理解,我已经知道了为什么,这里将我思考的结果分享一下。 1. 错误的理解 某一个SNP,效应值(Effect)越大,就越显著??? 「事实上:」 效应值和显著性是两码事! 2. 用代码说话 GWAS分析中,最简单的就是一般线性模型(GLM),而GLM模 ...
个人分类: 农学统计|5163 次阅读|没有评论
R语言计算玉米测交试验的配合力
2020-9-15 20:52
1. 玉米测交试验 可以认为,玉米侧交试验,是一个NCII的试验,在玉米实际的育种中,由于其测用结合的特性,应用广泛。 试验设计如下: 骨干系:4个 待测系:7个 重复:2次 观测值:产量 「计算目的:」 1 计算待测系的一般配合力(选择最优品系) 2 计算待测系和骨干系的特殊配合力(选出最优组合) ...
个人分类: 农学统计|3510 次阅读|没有评论
R语音读取txt报错:UTF-16不支持 或者<ff><fe><63>\'多字节字符串有错
2020-9-14 22:00
1. 问题描述 一个txt文件,使用R中的 data.table 包中的 fread 函数读取时,报错: dat=fread( test.txt ) Error in fread( test.txt ): Fileisencoded in UTF-16,thisencodingisnot& ...
个人分类: R语言|8062 次阅读|没有评论
线性混合模型学习笔记2:公式解释
2020-9-12 10:11
上一章节中(http://blog.sciencenet.cn/blog-2577109-1195363.html),我们介绍了混合线性模型中的模型假定,这一章节,用一个动物育种的示例,介绍固定因子,随机因子以及相关矩阵的推导,明确什么是G矩阵,什么是R矩阵,以及他们和方差组分的关系。直和和直积也是重要的概念,特别是多个随机因子以及随机因子有互作时, ...
个人分类: 数量遗传学|6929 次阅读|没有评论
冷风骑车与贝多芬的命运交响曲
2020-2-19 22:54
上传感言: 这篇文章是年前骑车上班的路上有感而发,现在由于疫情封闭在家,远程办公状态一直不佳,想起年前在冷风中的感想,要振奋起来,相信疫情过后会有一个更好的发展。要有一股劲在胸中,要相信那些打不倒我们的会让我们更坚强。加油中国! 晨骑车去公司,一出门一阵冷风吹来。降温很明显。想要回去把自行车 ...
个人分类: 随感|2488 次阅读|没有评论
为什么体型比星座更能影响性格?
2019-10-26 22:03
1. 体型分为哪几种? 这个有很多流派,介绍一下克雷奇米的理论,把体型分为四种: 每个体型包括性格特征和易患病症。 矮胖型,具有狂躁气质,高血压,高血糖 瘦高型,具有分裂气质,精神分裂,自我封闭 强壮型,具有黏着气质,冲动,癫痫症 发育异常型,抑郁气质,抑郁症 换句话说,胖子十有八 ...
个人分类: 数量遗传学|2957 次阅读|没有评论
为什么相爱的人不能在一起?
2019-10-24 08:06
一个两天的bug,终于解决了,竟然是很基础的bug,花费了我两天时间,不免感叹:为什么相爱的人不能在一起。 事情是这个样子的: 最近测试一个R包 breedR 的动物模型的功能,用了它的测试数据: 1 library(breedR) 2 data( globulus ) 3 ped-globul ...
个人分类: 农学统计|2179 次阅读|没有评论
进步青年自我麻痹三大件:健身,编程,学英语
2019-10-21 21:39
1. 进步青年自我麻痹三大件 python 学习,从 print(hello world) 开始,到 print(bye bye world) 结束,或者学着,忘着,看书时激励自己,放弃时安慰自己,洗洗睡吧。 为什么? 因为没有目标,学习只是为了学习,是为了麻痹自己堕落状态下焦虑的心,非要做点什么,才配得上自己那进步青年的自我期许。 ...
个人分类: 农学统计|3849 次阅读|没有评论

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-3-29 22:52

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部