首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

文献笔记五十四:全基因组关联分析鉴定拟南芥中控制种子大小的调节因子

--maf 0.05 --recode --out 172sample_maf_filter 给每个snp位点添加ID bcftools view -h 172sample_maf_filter.recode.vcf...> head bcftools view -H 172sample_maf_filter.recode.vcf | grep "^#" -v | awk '{$3=$1"_"$2;print $0}'...beagle.25Nov19.28d.jar gt=172sample_maf_filter_snpID.vcf out=172sample_out ne=172 基因型填充需要很长时间,这里先用没有填充的做接下来的分析...SNP_Density.pvalue.jpg 这个rMVP包好厉害,一条命令MVP.Report(df)所有图就全部生成了 文章开头提到的两篇简书文章还用rMVP这个R语言做了全基因组关联分析,另外找时间来重复这部分内容...还有一部分内容是用R语言的glmnet做岭回归分析分析,这部分内容暂时还没有思路!

1.9K20
您找到你想要的搜索结果了吗?
是的
没有找到

R tips: rlang中的expression操作符

R中,library函数的表现有点特殊,传给它的参数变量不是类似于常规R表达式的即时执行,而是像是被‘冻结’了一样。...这是一个rlang中定义的一个操作符函数。 它的本质原理在于:R在运行代码时,会先将代码解析成叫做expression的中间态,然后再执行expression即可获取代码结果。...操作在tidyverse系列中很常见,比如可以将因子变量的水平值重新编码的函数fct_recode: ### 定义一个因子变量 test_factor <- factor(letters[1:5])...var_name =" 但是会报错,原因是因为在<em>R</em>中=操作符要求比较严格,如果是引号括起来就<em>没有</em>问题了,但是括起来的时候,!!...var_name := 1))) #$test #[1] 1 列表可以使用list2函数,但是如果是向量的话,rlang<em>包</em>是<em>没有</em>c2函数的,这个时候可以先用list2处理,然后unlist函数转换为向量

1.4K10

R数据科学|第十一章内容介绍

使用forcats处理因子 因子在 R 中用于处理分类变量。从历史上看,因子远比字符串更容易处理。因此,R 基础中的很多函数都自动将字符串转换为因子。...准备工作 使用forcats来处理因子,这个提供了能够处理因子的工具,其中还包括了处理因子的大量辅助函数。...x2 <- c("Dec", "Apr", "Jam", "Mar") (2) 其对月份的排序没有意义。...修改因子水平 可以使用fct_recode()函数,它可以对每个水平进行修改或重新编码。该函数会让没有明确提及的水平保持原样,如果不小心修改了一个不存在的水平,那么它也会给出警告。...我们将其修改为较为详细的排比结构: gss_cat %>% mutate(partyid = fct_recode(partyid, "Republican, strong" = "Strong republican

57320

Linux常用命令笔记

php5-intl php-pear php5-imagick php5-imap php5-mcrypt php5-memcache php5-ming php5-ps php5-pspell php5-recode...3274 tar z- 用gzip对存档压缩或解压 x-从存档展开文件 v-详细显示处理的文件 f-指定存档或设备 chmod + 更新源: sudo apt-get update 更新已安装的sudo...apt-get upgrade 账号和密码文件 ls -al /etc/passwd /etc/shadow -rw-rr– 1 root root 1611 2011-12-16 05:30 /etc.../passwd # 保存的是账号 -rw-r—– 1 root shadow 947 2011-12-16 05:30 /etc/shadow # 账号的密码等信息 /etc/passwd root:x...密码过期的宽恕时间:如果在5中设置的日期过后,用户仍然没有修改密码,则该用户还可以继续使用的天数 账号失效日期,过了这个日期账号就不能用了 保留的 以上面最后用户letuknowit为例(15400对应的日期为

95210

【学习】《R实战》读书笔记(第四章)

PPV课R语言读书会以“学习、分享、进步”为宗旨,通过成员协作完成R语言专业书籍的精读和分享,达到学习和研究R语言的目的。...一个工作例子 案例:研究男女管理者领导他们组织的差异。 数据集(数值5、4、3、2和1分别表示强烈反对、反对、中立、支持、强烈支持)如图1所示。...图3:R逻辑运算符号 基于逻辑运算符号关系来改变变量中内容,逻辑值为真,则可以修改。 拓展可以了解within()函数和car中的recode()函数。 重命名变量 重命名变量名,可以采用方法。...方法二:reshape中的rename()函数。 方法三:names()函数。 缺失值 数据集往往是不完整,因为各种缘由存有缺失值。 使用is.na()函数检查缺失值。.../ 2 《R in action》第一部分第四章 本栏目文章由PPV课R语言读书会提供,转载请注明来自PPV课R语言读书会。

84050

R语言从入门到精通:Day5

同时,有些操作符前后有空格,有些又没有,到底该怎么写才对?除此之外,R语言的换行和缩进似乎也不像python那么严格。那我们平时自己写代码的时候应该遵循什么样的规则呢?...这个文档是Google’s R Style Guide,众所周知,医生在临床上工作有临床指南,那程序员工作也有“指南”,上面的这个文档就是我们写R代码的指南。...话不多说,我们进入这次课程的主要内容,上节课中完成了R语言中的数据导入工作,那么这次课程的主要内容就是给大家介绍一些数据管理和操作的基本函数或语句。...其中,car中的recode()和doBy中的recodevar()、R语言中自带的cut(),这三个函数都是很受欢迎的变量重编码函数。...*plyr是一个集合了很多数据集操作函数的R,大家可以查看其帮助文档进一步学习。 ?

1.6K30
领券