Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

HAPPI.GWAS results not the same as GAPIT3 #1

Open
jrod55 opened this issue Sep 5, 2020 · 0 comments
Open

HAPPI.GWAS results not the same as GAPIT3 #1

jrod55 opened this issue Sep 5, 2020 · 0 comments

Comments

@jrod55
Copy link

jrod55 commented Sep 5, 2020

Hi,

I've run FarmCPU using the demo data provided using the command below, and run the same analysis using GAPIT within R, matching the same options as the Demo_FarmCPU.yaml file. Though the results are similar, there are in fact some differences.

Using HAPPI GWAS:
Rscript HAPPI_GWAS.R -GAPIT Demo_FarmCPU.yaml

Within R:
library(GAPIT3)
phenos = read.table("mdp_traits.txt", header = TRUE)
myGAPIT <- GAPIT(
Y=phenos,
Model.selection = TRUE,
model=c("FarmCPU"),
SNP.MAF = 0.05,
ncpus = 20,
file.G="mdp_genotype_chr",
file.Ext.G="hmp.txt",
file.from=1,
file.to=10,
file.path="./Demo/"
)

Is there a reason for differences in results if presumably the same model should be being fit?

Thank you,
Jonas

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant