StatPheno: Statistic analysis for phenotype data

In JunhuiLi1017/JM4QTN: Joint Mapping for QTN

Description

Statistic analysis for phenotype data, including normality test, analysis of variance and best linear unbiased estimation

Usage

StatPheno(phenoData, defineForm = NULL, effNotation = "G")

Arguments

phenoData	Phenotype data
defineForm	Formula in linear model defined by user
effNotation	Effect nanotation for best linear unbiased estimation

Details

Environment effect, block effect, replication and genotype effect are in order for the fouth column of phenoData

Value

normality_test	The result of shapiro.test
formula	The formula of the linear model
ANOVA	The result of ANOVA
lsmeans	The least squart means for best linear unbiased estimation

Author(s)

Junhui Li

Examples

phenoData <- matrix(NA,180,6)
phenoData <- as.data.frame(phenoData)
colnames(phenoData) <- c("E","B","R","G","Height","Weight")
phenoData[,1] <- as.factor(c(rep(c("A"),60),rep(c("B"),60),rep(c("C"),60)))
phenoData[,2] <- as.factor(c(rep(1,30),rep(2,30)))
phenoData[,3] <- 1
phenoData[,4] <- as.factor(c(1:30))
phenoData[,5] <- rnorm(180,175,5)
phenoData[,6] <- rnorm(180,75,10)
StatPheno(phenoData)

JunhuiLi1017/JM4QTN documentation built on June 4, 2019, 4:10 a.m.