Printing additional information for all SNPs with -snpout argument about bayesr HOT 3 CLOSED

Information for each saved MCMC sample is given in the snpout file. For example, if a row contains information for 100 SNPs the remaining 10400 SNPs have zero effects. However, the effect size of a SNP is given as variance explained (beta^2).
If you want betas you have to change the line in subroutine output_snploc from "write(ce,'(e15.6)') g(i)**2" to "write(ce,'(e15.6)') g(i)"

Cheers

from bayesr.

The SNP effects are estimated by the sample means of the posterior estimates and are given in the 'beta' column in file simout.param.

Cheers

from bayesr.

Thank you.
However, how can I obtain 1000 samples? Basically I have 10500 markers and I would like to obtain a matrix with 1000 x 10500, were each row is one sample from the posterior distribution.

Cheers

from bayesr.