[R] Relationship between covariance and inverse covariance matrices
Bert Gunter
gunter.berton at gene.com
Thu Dec 8 23:49:36 CET 2011
What does this have to do with R? Is this homework?
I suggest you post to some sort of math list. Perhaps others will have
more specific suggestions where.
-- Bert
On Thu, Dec 8, 2011 at 2:42 PM, Vivian Shih <vivs at ucla.edu> wrote:
> Hi,
>
> I've been trying to figure out a special set of covariance matrices that
> causes some symmetric zero elements in the inverse covariance matrix but am
> having trouble figuring out if that is possible.
>
> Say, for example, matrix a is a 4x4 covariance matrix with equal variance
> and zero covariance elements, i.e.
>
> [,1] [,2] [,3] [,4]
> [1,] 4 0 0 0
> [2,] 0 4 0 0
> [3,] 0 0 4 0
> [4,] 0 0 0 4
>
> Now if we let a[1,2]=a[2,1]=3, then the inverse covariance matrix will
> have nonzero elements on the diagonals as well as for elements [1,2] and
> [2,1]. If we further let a[3,4]=a[4,3]=0.5 then the indices of the nonzero
> elements of the covariance matrix also matches those indices of the inverse.
>
> The problem is, if any of the nonzero off-diagonal indices match, then
> the inverse covariance matrix will have non-matching nonzero elements. For
> example, if a[1,2]=a[2,1]=3 as before but now we'll let a[2,3]=a[3,2]=0.5,
> then a would be:
>
> [,1] [,2] [,3] [,4]
> [1,] 4 3.0 0.0 0
> [2,] 3 4.0 0.5 0
> [3,] 0 0.5 4.0 0
> [4,] 0 0.0 0.0 4
>
> The inverse covariance matrix is now:
> [,1] [,2] [,3] [,4]
> [1,] 0.58333333 -0.44444444 0.05555556 0.00
> [2,] -0.44444444 0.59259259 -0.07407407 0.00
> [3,] 0.05555556 -0.07407407 0.25925926 0.00
> [4,] 0.00000000 0.00000000 0.00000000 0.25
>
> If we let a[1,2] and a[2,3] be nonzero, then the inverse will create a
> nonzero [1,3]. Does that happen all the time? I've tried to write out the
> algebraic system of linear equations for a and a-inverse but couldn't come
> up with anything.
>
> Let me know if I'm not clear on anything. Basically I'd just like to see
> what type of covariance matrices will produce an inverse covariance matrix
> with some zero elements.
>
>
>
>
>
> Thanks,
> Vivian
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
--
Bert Gunter
Genentech Nonclinical Biostatistics
Internal Contact Info:
Phone: 467-7374
Website:
http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm
More information about the R-help
mailing list