Dear Gavin,the problem is that we are interested in those snps with very low maf. Checking the output I obtained by using snptest, snps with lower maf with respect to all of the other are confirmed in dbsnp as variants with low maf. How can I bypass the problem? Does snptest consider a threshold which I could modify?Exist other similar software on the web you could suggest me?Many thanksGiovanniTo unsubscribe from the list visit this webpage https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=OXSTATGEN&A=12013/11/6 Gavin Band <[log in to unmask]>
HI Giovanni,I think your results are as expected - SNPTEST is unable to fit the model at very rare SNPs.Thanks,Gavin.
On 6 Nov 2013, at 09:40, giovanni acierno wrote:To unsubscribe from the list visit this webpage https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=OXSTATGEN&A=1 <log><output.xlsx>Dear Gavin,attached is the log file deling with my last launch of snptest_v2.5-beta4.I attached here below some lines of my output file.Many thanksGiovanni2013/11/5 Gavin Band <[log in to unmask]>Hi Giovanni,To unsubscribe from the list visit this webpage https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=OXSTATGEN&A=1Re: -total_prob_limit: are you using the latest version (2.5-beta4)? If not please try that version.Please send log files / some lines of output if you want help diagnosing a problem.Thanks,Gavin.On 5 Nov 2013, at 15:08, giovanni acierno wrote:To unsubscribe from the list visit this webpage https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=OXSTATGEN&A=1Dear colleagues,I'm using the second version (V2) of Snptest software. Starting from imputed snps obtained from Impute2 analysis (80.000 imputed snps starting from 1071), I have changed the -total_prob_limit parameter in a range from 0.001 to 1 paired with -method score and -method em respectively. No differences in the output files sum up. Impute2 output corresponds to region 2q35 using hg19 human genome. Out of about 80.000 snps ranging in that region, more than a half have not a related p-value. Is in your experience a recurrent issue? Is the fraction of discarded snps too high?Many thanks,--------------------------------------------------------------Dr. Gavin BandWellcome Trust Centre for Human GeneticsUniversity of Oxford--------------------------------------------------------------
--------------------------------------------------------------Dr. Gavin BandWellcome Trust Centre for Human GeneticsUniversity of Oxford--------------------------------------------------------------
To unsubscribe from the list visit this webpage https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=OXSTATGEN&A=1