Last data update: 2014.03.03

R: Variance-covariance matrix
vcov.clmeR Documentation

Variance-covariance matrix

Description

Extracts variance-covariance matrix for objects of class clme.

Variance-covariance matrix

Usage

## S3 method for class 'clme'
vcov(object, ...)

## S3 method for class 'summary.clme'
vcov(object, ...)

Arguments

object

object of class clme.

...

space for additional arguments

Value

Numeric matrix.

See Also

CLME-package clme

Examples

data( rat.blood )
cons <- list(order = "simple", decreasing = FALSE, node = 1 )
clme.out <- clme(mcv ~ time + temp + sex + (1|id), data = rat.blood ,
                 constraints = cons, seed = 42, nsim = 0)

vcov( clme.out )

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(CLME)
Loading required package: shiny
Loading required package: lme4
Loading required package: Matrix

Attaching package: 'CLME'

The following object is masked from 'package:lme4':

    VarCorr

The following object is masked from 'package:stats':

    sigma

> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/CLME/vcov.clme.Rd_%03d_medium.png", width=480, height=480)
> ### Name: vcov.clme
> ### Title: Variance-covariance matrix
> ### Aliases: vcov.clme vcov.summary.clme
> 
> ### ** Examples
> 
> data( rat.blood )
> cons <- list(order = "simple", decreasing = FALSE, node = 1 )
> clme.out <- clme(mcv ~ time + temp + sex + (1|id), data = rat.blood ,
+                  constraints = cons, seed = 42, nsim = 0)
> 
> vcov( clme.out )
               time0 Hour time6 Hour time24 Hour time48 Hour time72 Hour
time0 Hour   3.593126e-01  0.3165738   0.3165738   0.3165738   0.3165738
time6 Hour   3.165738e-01  0.3432856   0.3219162   0.3219162   0.3219162
time24 Hour  3.165738e-01  0.3219162   0.3432856   0.3219162   0.3219162
time48 Hour  3.165738e-01  0.3219162   0.3219162   0.3432856   0.3219162
time72 Hour  3.165738e-01  0.3219162   0.3219162   0.3219162   0.3432856
tempRT       8.554270e-16 -0.0106847  -0.0106847  -0.0106847  -0.0106847
sexMale     -3.213226e-01 -0.3213226  -0.3213226  -0.3213226  -0.3213226
                   tempRT       sexMale
time0 Hour   8.506160e-16 -3.213226e-01
time6 Hour  -1.068470e-02 -3.213226e-01
time24 Hour -1.068470e-02 -3.213226e-01
time48 Hour -1.068470e-02 -3.213226e-01
time72 Hour -1.068470e-02 -3.213226e-01
tempRT       2.136939e-02 -1.155694e-15
sexMale     -1.149595e-15  6.426451e-01
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>