corandcov: Provides the Variance-covariance parameters estimates

In emilelatour/lamisc: Latour's Little Helpers

Provides the Variance-covariance parameters estimates

Description

User-written function to output the covariance parameters

Usage

corandcov(glsob, cov = TRUE, ...)

Arguments

glsob	An object of class "gls" representing the linear model fit
cov	Default is TRUE. Not sure how this is used
...	Additional arguments

Value

A list with correlation matrix, variance structure, and covariance matrix.

References

https://stats.idre.ucla.edu/r/examples/alda/r-applied-longitudinal-data-analysis-ch-7/

Examples

## Not run: 
library(dplyr)
library(readr)
library(nlme)

opposites <- readr::read_csv(
  file = "https://stats.idre.ucla.edu/stat/r/examples/alda/data/opposites_pp.txt",
  col_names = TRUE)
dplyr::glimpse(opposites)

unstruct <- gls(opp ~ time * ccog,
                opposites,
                correlation = corSymm(form = ~ 1 |
                                        id),
                weights = varIdent(form = ~ 1 | wave),
                method = "REML")

corandcov(unstruct)

## End(Not run)

emilelatour/lamisc documentation built on May 20, 2024, 2:42 a.m.