Comments (5)
The software does not support fitting random effects other than marker effects. Have a look at the BGLR R-package which supports fitting additional random effects.
from bayesr.
All effects are considered "random" in the Bayesian framework. "Fixed" effects are given flat priors (no shrinkage of effects) and estimates are not updated with samples from the posterior distribution (estimates based on information contained in the likelihood only). Informative priors are only implemented for marker effects.
from bayesr.
@syntheke Sorry, you did not answer my question.
I wanna include a design matrix for year (as factor variable) in order to estimate the specific year effects, and another model matrix for blocks in order to estimate one variance component (which means I wanna add a random effect on the intercept). I managed to include the design matrix for fixed effects (year effect) using the covar argument. How can I include a design matrix for block as random effect? Please, is there anyway to do that?
Thank you
from bayesr.
@syntheke Thanks for responding to @leticia-lara
We would like to add an effect with a normal prior where we estimate variance component - so a "random" effect with an "informative" normal prior - not the flat prior.
Is this possible?
Thanks!
from bayesr.
-numit 100000 -burnin 1000 -thin 10
Dear,
Recently I am running BayesR. I have some questions about computation time of bayesR. I'm particularly curious about how long your data took.
About my data, There are 202840 SNPs and 8438 subjects. I include three covariates.
In order to reduce the computation burden, according to the BayesRMannual published in Github (https://github.com/syntheke/bayesR), I used the “500SNPs” strategy by modifying the two options(--msize 500 -mrep 5000). Referring to the manual about time, I infer that my data would take less than six hours. Actually I run the first step about 120 hours, and the first step didn’t finish. I really can't find the reason, so please give me some guidance and suggestions about my code. Thank you very much.
This is the command I used:
##train a model
./bayesRv2 -bfile traindata1 -out traindata1
-numit 50000 -burnin 20000 -seed 333
-blocksize 4 -nthreads 4
-msize 500 -mrep 5000
-covar cov_train1.txt
##prediction
./bayesRv2 -bfile testdata1 -out testdata1
-predict -model traindata1.model -freq traindata1.frq
-param traindata1.param -covar cov_test1.txt
-alpha traindata1.alpha
from bayesr.
Related Issues (16)
- direction of effect estimates HOT 3
- Using -dfvara -3.0 HOT 1
- Printing additional information for all SNPs with -snpout argument HOT 3
- At line 1736 of file baymods.f90; Fortran runtime error: End of file
- BayesR: Some questions about computation time of bayesR. HOT 7
- fam file HOT 2
- allele frequency HOT 3
- Format/Readme file issue HOT 1
- Missing phenotypic values in .fam file
- How do I run BayesR with large memory ? HOT 1
- How to adjust the proportion of SNPs with no additive effect ? HOT 3
- Could I input covariate file into BayesR model? HOT 5
- Shape parameter must be positive HOT 1
- Way to do a weighted phenotype analysis with BayesR?
- End-of-file during read error, discrepancy between PLINK .bim file and simulated .bim file? HOT 13
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from bayesr.