| Attachment | Size |
|---|---|
 code_univar.R | 8.05 KB |
I have 75 metabolites. I am running a univariate ACE model analysis for each one to decompose the phenotypic variation into estimated genetic and environmental contributions.
Is it possible to adjust the code attached to run the model for each of these simultaneously to get the output (h, c and e), or is it only possible to run the model for each metabolite individually?
Thank-you
Printer-friendly version
Please clarify what you're trying to do here. Do you want to analyze each metabolite marginally, or analyze all the metabolites jointly (and estimate genetic and environmental sources of covariance among them)? That is, are you trying to automate doing 75 monophenotype analyses, or do you want to run 1 analysis of 75 phenotypes?
Thank-you for your reply.
I am trying to automate doing 75 monophenotype analyses, but I'm struggling to devise correct code at the moment.
Any suggestions would be greatly appreciated.