Comments (7)
syntheke
commented on August 13, 2024
Hi Jiefang,
looking at your command list, how did you perform "prediction" when fitting the model failed ( "first step" didn't complete")?
First you may confirm that the software is running as expected on your system using the examples provided.
Without access to the data its hard to provide guidance, I would try the following
- exclude covariates
- omit msize mrep
- reduce #SNPs
Cheers
from bayesr.
jiefangDuan
commented on August 13, 2024
Dear,
I am very grateful for your prompt reply and appreciate your valuable
suggestions.
About the data I use is UKBiobank that is imputed with mean .
About the command line I provided to you, is all my commands. I sent it
to you for completeness. Indeed, my fitting model failed, so the prediction
step was never performed. As for the suggestion you mentioned, before I
sent you my initial email, I have tried the operation of excluding
covariates, which also took a long time and was not completed. Of course,
the premise of these operations is that I used the "500SNPs" strategy. And
I never try to omit the options of msize and mrep because of the size of my
data. Next,I will try it. In terms of reducing the number of SNPs, it is
true that it can work with the case that includes 100,000 SNPs (the number
of SNPs have halved ) and doesn't include covariates. But reducing the
number of SNPs may be a last resort for me. Next, I'll try the suggestion
you mentioned and see the performance!
I don't have much experience in bioinformatics, so please feel free to
speak up if you have reminders.
Thank you very much! Best wishes for you!
Kind Regards,
Jiefang
syntheke ***@***.***> 于2021年12月2日周四 21:11写道：
… Closed #15 <#15>.
—
You are receiving this because you authored the thread.
Reply to this email directly, view it on GitHub
<#15 (comment)>, or
unsubscribe
<https://github.com/notifications/unsubscribe-auth/ATFZBC6VYTPLYQ3K35TBEATUO5V7BANCNFSM5JGVFNYA>
.
Triage notifications on the go with GitHub Mobile for iOS
<https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675>
or Android
<https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>.
from bayesr.
syntheke
commented on August 13, 2024
Hi Jiefang,
I'm not familiar with the Biobank data, but bayesR only accepts 0,1,2 as genotypes, missing genotypes (NA) are allowed.
Your data should run within a day, without using the reduced update.
The *.hyp file provides current estimates of the parameters during the run of the program and can be used to monitor the computations. You could set -burnin to zero and -numit to a low number first and check the the &.hyp file.
Cheers
from bayesr.
jiefangDuan
commented on August 13, 2024
Hi syntheke,
I'm sorry to bother you again.
According to your advice, I have tried the idea that setting -burnin to 0
and -numit to 500. It works and only takes 8 minutes. the defalut of
-nummit and -burnin can't finish within a day for my data. Now, I've
learned about the knowledge of MCMC, and I still don't know how to set the
parameters of -numit and -burnin for my data properly. All I know is that
the computing time should increase linearly with the number of iterations.
Can you give me some guidance on *how I can determine the suitable
parameter settings for my data* *about the options of -numit and -burnin*?
I will appreciate your advice! Looking forward to hearing from you!
Kind regards,
Jiefang
from bayesr.
syntheke
commented on August 13, 2024
Hi Jiefang,
did you look at the output files generated for the run with "-numit 500"?
from bayesr.
jiefangDuan
commented on August 13, 2024
Hi syntheke,
Sorry to reply to you late.
Indeed, I check the output files generated for the run with "-numit 500".
But I don't know how many I should set to -numit and -burnin is appropriate
for my application based on the output files.
Please forgive me for my limited English expression ability.
The following are the output files and their corresponding details.
[image: aaaaa.png]
[image: bbb.png]
1. traindata1.hyp
[image: 11111.png]
[image: 22222.png]
2. traindata1.log
[image: trainlof.png]
3. traindata1.model
[image: model.png]
4. traindata1.param
5. testdata1.log
[image: testlog.png]
syntheke ***@***.***> 于2021年12月6日周一 10:05写道：
… Hi Jiefang,
did you look at the output files generated for the run with "-numit 500"?
—
You are receiving this because you authored the thread.
Reply to this email directly, view it on GitHub
<#15 (comment)>, or
unsubscribe
<https://github.com/notifications/unsubscribe-auth/ATFZBCYVA23J3BIEW2DVANLUPQK7HANCNFSM5JGVFNYA>
.
Triage notifications on the go with GitHub Mobile for iOS
<https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675>
or Android
<https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>.
from bayesr.
syntheke
commented on August 13, 2024
Hi Jiefang,
I seems that your files were not uploaded correctly, I can only see the name of the files but no content.
Before continuing the conversation you need to
A. Test that the software runs on your system as expected.
Did you run the examples that are provided with the software and do results agree with those in the example folders?
B. Make sure that your data is formatted as required by bayesR
Running only a few hunded iterations and looking at the output files can sometimes highlight issues with the input data or model specifications. Obviously, if you're not familiar with bayesR you don't always know what to look for. However, you could try to run bayesR using a subset of randomly selected SNPs (as suggested earlier), for example 50k, and compare the derived model with your expectations or results from previous analysis (e.g. GREML).
Cheers
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
- How include random effects in the model HOT 5
- 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.