turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Specifying a covariance matrix. PROC GLIMMIX?

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

10-29-2009 11:18 AM

I have a dataset based upon observation from related individuals. I have a kinship matrix that estimates the pairwise relationship between every two individuals in the sample. My question is, how do I specify that matrix in a variance component model?

My best guess is PROC GLIMMIX. I understand how to specify the STRUCTURE of a covariance matrix, but that's not what I want to do. I already have the covariance matrix calculated, now I need to use it in a model.

Thanks for any help.

My best guess is PROC GLIMMIX. I understand how to specify the STRUCTURE of a covariance matrix, but that's not what I want to do. I already have the covariance matrix calculated, now I need to use it in a model.

Thanks for any help.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

10-29-2009 08:32 PM

KAZ,

Try Googling "+pedigree +covariance +matrix +proc +mixed" (without the quotes). The first hit should be "Genetic analysis of complex traits using SAS - Google Books Result". If you click on that, you can see a couple of pages of text in which the INBREED procedure is used to construct a kinship matrix. The kinship matrix is then passed to the MIXED procedure.

I would suggest that you purchase the book so that you can get the most from it. It looks like it could be quite useful

Dale

Try Googling "+pedigree +covariance +matrix +proc +mixed" (without the quotes). The first hit should be "Genetic analysis of complex traits using SAS - Google Books Result". If you click on that, you can see a couple of pages of text in which the INBREED procedure is used to construct a kinship matrix. The kinship matrix is then passed to the MIXED procedure.

I would suggest that you purchase the book so that you can get the most from it. It looks like it could be quite useful

Dale

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to Dale

11-03-2009 03:03 PM

Thanks Dale, that is very helpful. I have another related question:

I used PROC INBREED to construct the kinship matrix. Unfortunately I have a very large pedigree, which makes the matrix unwieldy. In reality I only need the kinship coefficient estimated for a small sample of pairwise relationships in the pedigree.

For example, I have observations on N individuals, who are connected via a single pedigree which contains Y individuals. Y >> N.

Can I tell PROC INBREED to use all Y individuals to estimate the kinship coefficients between only the N observed individuals?

Thanks.

I used PROC INBREED to construct the kinship matrix. Unfortunately I have a very large pedigree, which makes the matrix unwieldy. In reality I only need the kinship coefficient estimated for a small sample of pairwise relationships in the pedigree.

For example, I have observations on N individuals, who are connected via a single pedigree which contains Y individuals. Y >> N.

Can I tell PROC INBREED to use all Y individuals to estimate the kinship coefficients between only the N observed individuals?

Thanks.