Description Usage Arguments Details Value Note Author(s) See Also Examples

Extracts the (co)variance component and fixed effect estimates from an object of class `dmm`

, for the specified set of traits and set of components. Makes tables of component estimates ordered either by trait or by component. Tables include component estimate, its standard error, and its 95 percent confidence limits. If there are class specific components these appear with appropriate labels in the list of components.

1 2 3 |

`object` |
An object of class |

`traitset` |
A vector containing the names of the subset of traits for which tables of (co)variance component estimates are to be constructed. Default is "all" which means all traits present in object |

`componentset` |
A vector containing the names of the subset of (co)variance components for which tables are to be constructed. Default is "all" which means all (co)variance components present in object |

`bytrait` |
Logical flag: should the tables of (co)variance component estimates be constructed with trait varying least rapidly from line to line? If TRUE each subtable contains component estimates for one trait or traitpair and for all components in argument |

`gls` |
Logical flag: should the (co)variance component estimates by GLS-b method be tabled in addition to the (co)variance component estimates by OLS-b method? Default is |

`digits` |
Number of digits for output. This is returned as part of the return value for us
e by the S3 print function |

`...` |
Ellipsis argument. |

This is a long printout with estimates, standard errors and confidence limits, arranged in tables with one estimate per line. For a short printout see function `print.dmm()`

.
In the case of class-specific components, the components which are class-specific are labelled with their class-codes prepended to the variance component name, and components for all classes are included, so the listed components do not sum to phenotypic variance. To list the components in class groups, so that they sum to the class phenotypic variance use function `csummary()`

.

An object of class `gsummary.dmm`

which is a list containing the following
items:

`btables` |
A list of dataframe objects each containing one subtable of estimates of the fixed effects, along with the appropriate standard errors and confidence limits. Based on OLS-b fixed effect estimates. |

`ctables` |
A list of dataframe objects each containing one subtable of estimates of the (co)variance components, along with the appropriate standard errors and confidence limits. Based on OLS-b component estimates. |

`gbtables` |
A list of dataframe objects each containing one subtable of estimates of the fixed effects, along with the appropriate standard errors and confidence limits. Based on GLS-b fixed effect estimates. |

`gctables` |
A list of dataframe objects each containing one subtable of estimates of the (co)variance components, along with the appropriate standard errors and confidence limits. Based on GLS-b component estimates. Only present if argument |

`traits` |
A vector of traitnames as specified in argument |

`components` |
A vector of component names as specified in argument |

`bytrait` |
Logical flag: as specified in argument |

`gls` |
Logical flag: as specified in argument |

`digits` |
A numeric value, as specified in argument |

`call` |
The function call |

There is no provision to constrain the 95 percent confidence limits for component estimates. Hence for small samples, these may vary outside the bounds for the component, that is for components which are variances, they may be negative. Fixed effects are not bounded.
Use `summary()`

if you want to see the components as estimated. Use `csummrry()`

if you want to see the components summing to phenotypic (co)variance, or sorted into class-specific groups.

Neville Jackson

Function `print.summary.dmm()`

.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | ```
# get some data
data(sheep.df)
# prepare it - only need "E" and "A" relationship matrices
sheep.mdf <- mdf(sheep.df,pedcols=c(1:3),factorcols=c(4:6),ycols=c(7:9),
sexcode=c("M","F"),relmat=c("E","A"))
# estimate (co)variance components - individual and maternal
sheep.fit5 <- dmm(sheep.mdf, Ymat ~ 1 + Year + Sex,
components=c("VarE(I)","VarG(Ia)","VarE(M)","VarG(Ma)",
"CovG(Ia,Ma)","CovG(Ma,Ia)"))
# look just at component "VarG(Ma)" across all traits
summary(sheep.fit5,componentset="VarG(Ma)",bytrait=FALSE)
# look just at trait "Cww"
summary(sheep.fit5,traitset="Cww")
# cleanup
rm(sheep.df)
rm(sheep.mdf)
rm(sheep.fit5)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.