Hi,

I calculated genome-wide XP-EHH scores using 'Selscan' software and normalized the raw scores using the norm function provided in the "Selscan". I am super confused about how I can calculate the p-value for each normalized score.

I have read several papers, and they all just mentioned something like "significance test based on normal distribution" without showing the exact methods.

I am thinking of two solutions:

1. Calculated p-values using 'pnorm' in R
2. Treat the normalized scores as Z scores and compute a p-value from there.

Can anyone give some suggestions or solutions?

I really appreciate any help you can provide.