pin: Count error for h2 in breedR

In yzhlinscau/breedRPlus: Augments the Use of 'breedR' in Fitting Mixed Models

Description

pin This function counts standard error(se) for heritability(h2) and corr value and also outputs significent level for corr value in breedR package.

Usage

## S3 method for class 'remlf90'
pin(object, formula, signif = FALSE, digit = 5,
  vres = FALSE)

Arguments

object	an object of breedR result
formula	formula for h2 (or corr, not now)
signif	Index to output signif levels, F(default) for non-signif.
digit	a number to control decimal point,5(default).
vres	Index(T) to return results in vectors, F(default) for direct results.

Details

counts standard error(se) for heritability(h2)

Value

the result is returned directly.

Author(s)

Yuanzhen Lin <[e-mail]>

References

breedRPlus website:https://github.com/yzhlinscau/

Examples

library(breedR)
library(breedRPlus)

data(globulus)
res.animal <- remlf90(fixed = phe_X ~ 1,
                      random = ~ gg,
                      genetic = list(model = 'add_animal',
                      pedigree = globulus[, 1:3],
                      id = 'self'),
                      data = globulus)
                      
pin(res.animal, h2~V2/(V1+V2+V3))

yzhlinscau/breedRPlus documentation built on May 10, 2019, 1:54 a.m.