Multiple GRM for HE regression

jaleal
New Member

Posts: 5

Multiple GRM for HE regression Feb 25, 2016 0:47:37 GMT

Quote

Post by jaleal on Feb 25, 2016 0:47:37 GMT

Hello, is it possible to use --HEreg with --mgrm or does it not make sense to do HE regression with multiple grms?

I tried this:

gcta64 --HEreg --mgrm GRMLIST.txt --pheno PHENFILE.phen --out OUTFILE

Where GRMLIST.txt is a text file with the prefixes of the grm's I want to use. It seems to be treating GRMLIST.txt as the prefix rather than a file containing the prefixes. I get the error

Error: can not open the file [GRMLIST.txt.grm.id] to read.

Thanks,
Jaleal

Jian Yang Administrator Posts: 362	Multiple GRM for HE regression Feb 25, 2016 4:32:32 GMT Quote Select Post Deselect Post Link to Post Member Give Gift Back to Top Post by Jian Yang on Feb 25, 2016 4:32:32 GMT This hasn't implemented in GCTA but you can run it in R. A R function to read the GRM is available at gcta.freeforums.net/post/576/thread.

jaleal New Member Posts: 5	Multiple GRM for HE regression Feb 25, 2016 15:39:54 GMT Quote Select Post Deselect Post Link to Post Member Give Gift Back to Top Post by jaleal on Feb 25, 2016 15:39:54 GMT Yup, makes sense. Thanks

Post by jaleal on Feb 25, 2016 0:47:37 GMT

Post by Jian Yang on Feb 25, 2016 4:32:32 GMT

Post by jaleal on Feb 25, 2016 15:39:54 GMT