IASbreeding软件包
IASbreeding V1.1.0 can calculate kinship matrix, variance component analysis using linear mixed model, breeding value using BLUP. The core algorithm of IASbreeding is built using Rcpp and RcppEigen.
安装IASbreeding软件包
先下载安装R和(或)Rstudio用于执行R语言软件包
install.packages("IASbreeding_1.1.0.zip", repos = NULL)
library(IASbreeding)
IASbreeding软件包使用说明
phe_file <- system.file("extdata", "Phenotype.txt", package = "IASbreeding", mustWork = TRUE)
pheno <- read.table(phe_file,header=TRUE)
ped_file <- system.file("extdata", "Pedigree.txt", package = "IASbreeding", mustWork = TRUE)
ped <- read.table(ped_file,header=TRUE)
newped<-pedcheck(ped)
A<-Akin(newped)
Results <- IASBLUP(BW56 ~ 1 + Sex + Condition,random = "ID", kinship = A, max_iter = 20, data = pheno)
bed_file <- system.file("extdata", "Genotype.bed", package = "IASbreeding", mustWork = TRUE)
file_prefix <- sub("\\.bed$", "", bed_file)
Genotype <-read_bed(file_prefix)
G<-Gkin(Genotype)
Results <- IASBLUP(BW56 ~ 1 + Sex + Condition,random = "ID", kinship = G, max_iter = 20, data = pheno)
- 基于系谱+基因型信息构建H矩阵,运行single-step GBLUP(ssGBLUP)
H<-Hkin(A, G)
Results <- IASBLUP(BW56 ~ 1 + Sex + Condition,random = "ID", kinship = H, max_iter = 20, data = pheno)