stat1 {qgen}R Documentation

Statistics for one trait

Description

Calculates the narrow-sense heritability and its confidence intervals for one trait.

Usage

stat1(paraDATA,
      alpha=0.05,
      frommethod="REML",
      Satterthwaite=TRUE,
      tex.table=FALSE,
      file=FALSE,
      path="~/qgen/")

Arguments

paraDATA object of class paraDATA created by the function est (or emp, see below)
alpha two sided error probability
frommethod If the paraDATA–object provides results form different variance–partioning methods you can choose which to take otherwise the REML solutions taken without warning:
    "ANOVAuw"
    ANOVA estimators with unweighted sums of squares
    "ANOVA"
    normal ANOVA estimators
    "REML"
    REML solutions
Satterthwaite Logical flag: Should the approximated degree of freedoms for sire and dam be used (only possible if they are available)?
file Logical flag: Should the object be written to
a file called "stat1.rda"?
tex.table Logical flag: should a LaTeX file be created with a table of all estimates
path the path for the output files

Details

This function calculates the heritability of individual traits. If the paraDATA provides parameters for several characters or several environments or both, then for every character–trait combination heritability is calculated individualy. If the paraDATA has a unbal slot with data about the balancedness of the data set then also the confidence intervals are calculated. If the paraDATA provides parameters estimated with diffent variance partitioning methods then it is possible to choose among them.

Value

an object of the class stat.

Note

Attributes: limes

References

Burdick, R. K. and Graybill, F. A. (1992) Confidence intervals on variance components. Marcel Dekker, New York.

Sen, B., Graybill, F. A. and Ting, N. (1992) Confidence intervals on variance components. Biometrical Journal 3, 259–274.

See Also

(not yet implemented) stat2e, for the same statistics of an index of two environments; statXt, for several traits within one environment.

Examples

## calculates narrow-sense heritability and its confidence intervals
myparameters <- est(sim(the(file=FALSE)), file=FALSE)  # takes the default causal parameters of the()
                                                       # simulates a data set
                                                       # estimates the variance components
stat1(myparameters)               # calculates heritability and
                                  # confidence invals from estimated parameters

[Package qgen version 0.03-02 Index]