[GenABEL-dev] Precision and scientific notation in ProbABEL
L.C. Karssen
lennart at karssen.org
Wed Aug 14 10:03:45 CEST 2013
Thanks for the input Nicola, I've set the output to 6 significant
digits, although it seems that the cout function (used for printing)
sometimes "eats" a trailing zero.
Best,
Lennart.
On 30-07-13 10:33, Nicola Pirastu wrote:
> Dear Lennart,
>
> I think that switching to scientific notation is not really necessary and could lead to a little of loss in precision unless of course you still
> use 6 significant digits which will translate in just a reduction of 0 in the values.
> So if for example we were to choose scientific notation with 3 significant digits, although this would not affect very much the final results we could be asked to submit more and
> would not be able to comply.
> So to summarize I think that if it does not have any effect on performance of ProbABEL 6 significant digits without scientific notation is fine.
>
> Best
>
> Nicola
>
>
> Dr. Nicola Pirastu PhD
> Research Fellow
> Medical Sciences, Chirurgical and Health Department
> University of Trieste
> Medical Genetics
> IRCCS Burlo Garofolo
> Via dell'Istria 65/1
> 34137 Italy
> tel. +390403785539
>
> Il giorno 30/lug/2013, alle ore 09:32, "L.C. Karssen" <lennart at karssen.org> ha scritto:
>
>> Dear list,
>>
>> I'm finalising version 0.4.0 of ProbABEL and there are two things I'd
>> like your opinion on:
>>
>> 1) with what precision should we print the betas, standard errors and
>> Chi^2 values to the output files?
>>
>> 2) Should we use scientific notation in the output (for betas, standard
>> errors and Chi^2)?
>>
>> In ProbABEL v0.3.0 and earlier output was simply sent to cout without
>> any explicit formatting. In practice this lead usually to 6 significant
>> digits, but sometimes less. My proposal is to fix the precision at 6
>> significant digits.
>>
>> Regarding item 2): most of the betas I see are in the range between 0
>> and 10, although in case of no effect beta's can be of the order of
>> 1e-2, 1e-3. All in all, I don't think switching to scientific notation
>> will improve the output.
>>
>>
>> What are your opinions?
>>
>>
>> Thanks,
>>
>> Lennart.
>> --
>> -----------------------------------------------------------------
>> L.C. Karssen
>> Utrecht
>> The Netherlands
>>
>> lennart at karssen.org
>> http://blog.karssen.org
>>
>> Stuur mij aub geen Word of Powerpoint bestanden!
>> Zie http://www.gnu.org/philosophy/no-word-attachments.nl.html
>> ------------------------------------------------------------------
>>
>> _______________________________________________
>> genabel-devel mailing list
>> genabel-devel at lists.r-forge.r-project.org
>> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/genabel-devel
>
> AVVISO DI RISERVATEZZA Informazioni riservate possono essere contenute nel messaggio o nei suoi allegati. Se non siete i destinatari indicati nel messaggio, o responsabili per la sua consegna alla persona, o se avete ricevuto il messaggio per errore, siete pregati di non trascriverlo, copiarlo o inviarlo a nessuno. In tal caso vi invitiamo a cancellare il messaggio ed i suoi allegati. Grazie. CONFIDENTIALITY NOTICE Confidential information may be contained in this message or in its attachments. If you are not the addressee indicated in this message, or responsible for message delivering to that person, or if you have received this message in error, you may not transcribe, copy or deliver this message to anyone. In that case, you should delete this message and its attachments. Thank you.
>
--
-----------------------------------------------------------------
L.C. Karssen
Utrecht
The Netherlands
lennart at karssen.org
http://blog.karssen.org
Stuur mij aub geen Word of Powerpoint bestanden!
Zie http://www.gnu.org/philosophy/no-word-attachments.nl.html
------------------------------------------------------------------
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 230 bytes
Desc: OpenPGP digital signature
URL: <http://lists.r-forge.r-project.org/pipermail/genabel-devel/attachments/20130814/3552c5e3/attachment.sig>
More information about the genabel-devel
mailing list