Tag: qqplot

python – Getting AIC of lognormal distributions using scikit’s GMM

The short answer is that you should not compare aic of models that have different scales of response variables. See post like this or this. So for example you can compare gmm = GaussianMixture(n_components = 1).fit(X=np.expand_dims(data,1)) print(gmm.aic(X=np.expand_dims(data,1))) gmm2 = GaussianMixture(n_components = 2).fit(X=np.expand_dims(data,1)) print(gmm2.aic(X=np.expand_dims(data,1))) 18753.718210403502 17498.87614003799 You can see that the…

Continue Reading python – Getting AIC of lognormal distributions using scikit’s GMM

How to make Q-Q plots for different models at a single chart using R?

How to make Q-Q plots for different models at a single chart using R? 1 Hello everyone! I am trying to make a Q-Q plot to visualise p-values obtained from different models afterr runing GWAS analysis. The key thing is that I want to reflect p-values from the different models…

Continue Reading How to make Q-Q plots for different models at a single chart using R?

Cool Bioinformatics Scripts With Python

You can use this script in two ways read tons of millions of P values from stdin # python zcat pval.txt.gz | qqplot.py -out test -title “QQ plot on the fly” # julia zcat pval.txt.gz | qqplot.jl –out test –title “QQ plot on the fly” warning : If you have…

Continue Reading Cool Bioinformatics Scripts With Python

GitHub – Zilong-Li/BioScripts: Cool Bioinformatics Scripts

You can use make a QQ plot in the following ways. one-liner for reading tons of millions of P values from the pipe # python zcat pval.txt.gz | qqplot.py -out test -title “QQ plot on the fly” # julia (recommand to run it in the REPL) zcat pval.txt.gz | qqplot.jl…

Continue Reading GitHub – Zilong-Li/BioScripts: Cool Bioinformatics Scripts

Cool Bioinformatics Scripts | PythonRepo

You can use this script in two ways read tons of millions of P values from stdin # python zcat pval.txt.gz | qqplot.py -out test -title “QQ plot on the fly” # julia zcat pval.txt.gz | qqplot.jl –out test –title “QQ plot on the fly” warning : If you…

Continue Reading Cool Bioinformatics Scripts | PythonRepo

Dealing with P value inflation of SNPs identified by qqplot

Dealing with P value inflation of SNPs identified by qqplot 0 I’m attempting to call SNPs that changed in allele frequency in a population. The population was sequenced at two times (before and after a drought) using pooled-sequencing of the whole genome. Fisher exact tests identified 7 million SNPs, which…

Continue Reading Dealing with P value inflation of SNPs identified by qqplot