RGE (version 1.0)

RGEgibbs: Response from Genotype to Environment

Description

This function performance the GIBBS sampler for analyze reginals trials.

Usage

RGEgibbs(data, gen_c, env_c, blk_c, y_c, prior.g = NULL,
        prior.vg =NULL,prior.b = NULL, prior.dfb = NULL,
        prior.sv = NULL, prior.dfsv = NULL, prior.se = NULL, 
        prior.dfse = NULL,
        burnin = 10, thin = 5, niter = 50, saveAt = 10)

Arguments

data

data.frame

gen_c

Number of the column from de data.frame with the genotypes information.

env_c

Number of the column from de data.frame with the environment information.

blk_c

Number of the column from de data.frame with the block information.

y_c

Number of the column from de data.frame with the phenotype information.

prior.g

Vector with prior information of the means of genotypes

prior.vg

Vector with prior information of the variances of the means of genotypes

prior.b

Vector with prior information of the variances of block within environment

prior.dfb

Vector with prior information of hyperparameter degree of credibility of the variances of block within environment.

prior.sv

Vector with prior information of the stability's variances

prior.dfsv

Vector with prior information of hyperparameter degree of credibility of the stability's variances.

prior.se

Vector with prior information of the error's variances

prior.dfse

Vector with prior information of hyperparameter degree of credibility of the variances of error.

burnin

Number of iteration to be consider as burn-in period. This period is not saved in the final result.

thin

The thinning interval between consecutive observations. This interval is not saved in the final result.

niter

Numbers of iterations to be saved

saveAt

Save object outtS4 with samples of the posterior distribution on the work directory each "saveAt" iteration

Value

Matrix with samples of the posterior distribution

References

Cotes Torres, J. M., Gonzalez Jaimes, E. P., & Cotes Torres, A. (2016). Seleccion de Genotipos con Alta Respuesta y Estabilidad Fenotipica en Pruebas Regionales: Recuperando el Concepto Biologico. Revista Facultad De Ciencias Basicas, 8(2), 226-243.

Examples

# NOT RUN {
##data(datos)
##m<-RGEgibbs(data=datos,gen_c=5,env_c=1,blk_c=4,y_c=9,
##thin=5,burnin=100,niter=10000,saveAt=1000)
# }