pin: Count error for h2 and corr.

In yzhlinscau/AAfun0s: Added Functions for Echidna

Count error for h2 and corr.

Description

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

Usage

pin(
  object,
  mulp = NULL,
  idx = 1,
  digit = 3,
  all = FALSE,
  signif = FALSE,
  Rres = FALSE
)

## S3 method for class 'esR'
pin(
  object,
  mulp = NULL,
  idx = 1,
  digit = 3,
  all = FALSE,
  signif = FALSE,
  Rres = FALSE
)

pin11(object,formula=NULL,idx=1,signif=FALSE,
                digit=3,all=FALSE,Rres=FALSE)

pin22(...,mulp,signif=FALSE,digit=3)

Arguments

object	Echidna result object in R.
mulp	multi-formula for h2, corr, etc.
idx	trait order to run, 1,2...
digit	Index for decimal number, 3(default).
all	Show variance component and genetic parmameter together, FALSE(default).
signif	Index to output signif levels, FALSE(default) for non-signif.
Rres	Index(TRUE) to restore results, FALSE(default).
formula	formula for h2 or corr.

Details

Count error for h2 and corr value, also outputs significent level.

Value

the result is returned directly.

Author(s)

Yuanzhen Lin <[e-mail]>

Examples

## Not run: 

 library(AFEchidna)

 AFEchidna::read.example(package='AFEchidna', setpath=TRUE)
 file.edit("fm.es0")
 
 res<-echidna(h3~1+Rep,
              random=~Fam,es0.file="fm.es0")
 
 # pin results if using vpredict function
 pin(res)
 
 # run pin function to count genetic parameters
 pin11(res,h2~V1/(V1+V2))
 pin(res,mulp=c(h2~V1/(V1+V2),h2f~V1/(V1+V2/4))

## End(Not run)

yzhlinscau/AAfun0s documentation built on April 18, 2023, 4:11 p.m.