netmeta (version 3.1-1)

invmat: Moore-Penrose Pseudoinverse of a Matrix

Description

Calculates the Moore-Penrose pseudoinverse of a square matrix X.

Usage

invmat(X)

Value

The Moore-Penrose pseudoinverse for matrix X.

Arguments

X

A square matrix.

Author

Gerta Rücker , Guido Schwarzer

Details

This function is used by default in R package netmeta to calculate the Moore-Penrose pseudoinverse L\(^+\) of the Laplacian matrix L (Rücker, 2012):

L\(^+\) = (X - J / n)\(^{-1}\) + J / n with identity matrix J of dimension nxn.

References

Rücker G (2012): Network meta-analysis, electrical networks and graph theory. Research Synthesis Methods, 3, 312--24

See Also

netmeta, solve

Examples

# \donttest{
data(smokingcessation)

pw1 <- pairwise(list(treat1, treat2, treat3),
  event = list(event1, event2, event3), n = list(n1, n2, n3),
  data = smokingcessation, sm = "OR")
net1 <- netmeta(pw1)

invmat(net1$L.matrix.common)

data(Senn2013)

net2 <- netmeta(TE, seTE, treat1, treat2, studlab, data = Senn2013)
L1 <- net2$L.matrix.common
L2 <- invmat(net2$Lplus.matrix.common)
all.equal(round(L1, 10), round(L2, 10))
# }