scdrs compute-score output

[NaomiHuntley]

Hello and thank you for this great tool. I am running in to an issue at the scdrs compute-score step.

Each of my traits are a separate .gs file at this point. Here is an example of the first few lines of one of the trait files:

Here is the code I am running:

But my output files are both not as expected:

I would really appreciate any insight that you have on this issue!

Thank you!

[martinjzhang]

Hi @NaomiHuntley ,

.gz file is a compressed file, so head will not work. You can use zcat to read the compressed binary file and then use head to print out the first few lines. Can you try zcat [file_name] | head instead?

[NaomiHuntley]

Hi @martinjzhang oh wow, what a noob mistake! Thank you for that information.

One other issue I am running in to is that since each of my traits is in a separate file, I am running the compute-score command on a loop. However, my output is only one file. Does the command automatically concatenate all the results? Or should I have a separate output file for each of my input traits.

Thank you!

[martinjzhang]

Hi @NaomiHuntley ,

I think you named all your traits ZSTAT in the .gs file. Since scdrs specifies output file names using the trait name (which is ZSTATS for all your traits), all your traits may have the same output file names. So I think the new output results overwrote the old ones.

[NaomiHuntley]

That was the issue! Thank you