outlierCARBayesST: outlierCARBayesST

In sanchezi/openSilexStatR: Lattice Experiments Data Analysis

Description

a function to detect outlier point in time-course using a CARBayesST modelling

Usage

outlierCARBayesST(modelin, datain, threshold, trait)

Arguments

modelin	a CARBayesST object
datain	input dataframe of parameters used in the CARBayesST modelling
threshold,	numeric threshold for influence with standardized residuals
trait,	character, trait of interest to model (example plantHeight...)

Value

a plot of fitted vs. residuals and a dataset with a column identifying the outlier points

See Also

CARBayesST, ST.CARanova

Examples

 data(plant1)
 library(dplyr)
 mydata<-filter(plant1,!is.na(plant1$thermalTime))
 str(mydata)
 test<-fitCARBayesST(datain=mydata,xvar="thermalTime",
       trait="plantHeight",k=2,
       graphDist=TRUE,burnin=10,n.sample=110,
       formulaModel=as.formula(plantHeight~scenario+genotypeAlias),
       typeModel="anova",verbose=FALSE)
 test2<-outlierCARBayesST(modelin=test[[1]],datain=test[[2]],threshold=4,
        trait="plantHeight")

sanchezi/openSilexStatR documentation built on Sept. 10, 2020, 1:03 p.m.