Apply Solve to an reStruct Object

Solve is applied to each pdMat component of a, which results in inverting the positive-definite matrices they represent.

# S3 method for class 'reStruct'
solve(a, b, ...)

Arguments

a

an object inheriting from class "reStruct", representing a random effects structure and consisting of a list of pdMat objects.

b

this argument is only included for consistency with the generic function and is not used in this method function.

...

some methods for this generic require additional arguments. None are used in this method.

Value

an reStruct object similar to a, but with the pdMat components representing the inverses of the matrices represented by the components of a.

Author

José Pinheiro and Douglas Bates bates@stat.wisc.edu

See also

solve.pdMat, reStruct

Examples

rs1 <- reStruct(list(A = pdSymm(diag(1:3), form = ~Score),
  B = pdDiag(2 * diag(4), form = ~Educ)))
solve(rs1)
#> Random effects:
#>  Formula: ~Score | A
#>  Structure: General positive-definite
#>    StdDev    Corr   
#> V1 1.0000000 V1   V2
#> V2 0.7071068 0      
#> V3 0.5773503 0    0 
#> 
#>  Formula: ~Educ | B %in% A
#>  Structure: Diagonal
#>                V1        V2        V3        V4 Residual
#> StdDev: 0.7071068 0.7071068 0.7071068 0.7071068        1
#>