R qtl package 使用记录

1. qtl安装了很多次，提示R版本不符合；重启了R，安装成功

2. 画出的图如下，更详细的linkage map 图需要借助其它软件

########################

install.packages("qtl")

library(qtl)

source("https://bioconductor.org/biocLite.R")

install.packages("BiocManager")
BiocManager::install("qtl")
if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager")
##########################################

data(fake.bc)
ls()
summary(fake.bc)
plot.missing(fake.bc)
fake.bc <- drop.nullmarkers(fake.bc)
totmar(fake.bc)
fake.bc <- calc.genoprob(fake.bc, step=1, error.prob=0.01)
out.em <- scanone(fake.bc)
out.hk <- scanone(fake.bc, method="hk")
fake.bc <- sim.geno(fake.bc, step=2, n.draws=16, error.prob=0.01)
out.imp <- scanone(fake.bc, method="imp")
summary(out.em)
summary(out.em, threshold = 3)
summary(out.hk, threshold = 3)
summary(out.imp, threshold = 3)

max(out.em) # based on expectation maximization

max(out.hk) # based on Haley-Knott regression

max(out.imp) # based on multiple imputation

plot(out.em, chr=c(2,5))
plot(out.em, out.hk, out.imp, chr=c(2,5))

R qtl package 使用记录相关推荐

R语言数据处理学习记录--使用linkET包完成mental test
R语言数据处理学习记录–使用linkET包完成mental test 注:本文仅作为自己的学习记录以备复习查阅一简单介绍一下这个包依旧是来自于GitHub上,作者给出的介绍:linkET的目标是 ...
使用R/qtl进行QTL分析
使用R/qtl进行QTL分析 QTL分析是进行基因精细定位和克隆的基础,今天教大家使用R包" qtl "进行QTL分析. 在开始分析前,我们需要准备两个输入文件:基因型和表型文件. ...
R语言报错记录The following objects are masked from ‘package:stats’:decompose, spectrum
报错记录
R 语言学习过程全记录 ~
RStudio介绍超详细的教程:https://www.jianshu.com/p/132919ca2ca9 前辈的心得:https://blog.csdn.net/kMD8d5R/article/d ...
R语言基础学习记录4：重要函数
时间: 2018-07-18(学习) 2018-07-22(学习记录) 教程:慕课网 <R语言基础> 讲师:Angelayuan 补充内容: R语言常用函数总结大全.gl()函数学习内容 ...
R语言package安装显示(as ‘lib’ is unspecified)并出现未知路径
想问一个新问题: 安装本地包"dataMAM"package的时候出现这个 > install.packages("G:/Study software/R pack ...
R中package not avalible问题
我们在加载之前安装的R包时,有时会出现加载失败的情况.常规做法是remove.packages(""),然后重新安装,但是有时候会出现安装失败的情况,如下: 1. 针对not wr ...
Ubuntu16.04 配置记录（持续更新）
Ubuntu16.04 配置记录 1.安装中文输入法 https://www.cnblogs.com/darklights/p/7722861.html 2.修改终端命令行路径长度 https://b ...
R+Hadoop大数据方案有哪些坑？
摘要: 为什么有些公司在机器学习业务方面倾向使用 R + Hadoop 方案?因为他们在不懂R和Hadoop的特征应用场景的情况下,恰好抓到了一根免费,开源的稻草.R:R的应用场景不在于无与伦比的统计 ...

R qtl package 使用记录

R qtl package 使用记录相关推荐

最新文章

热门文章