lh3 / psmc

Implementation of the Pairwise Sequentially Markovian Coalescent (PSMC) model
Other
146 stars 60 forks source link

Confidence interval for PSMC #20

Open Homap opened 4 years ago

Homap commented 4 years ago

Hello,

I have a question about calculating confidence interval for the PSMC method. I have done the bootstrapping step. To calculate the confidence interval, I am using the maximum error equation for a t distribution:

E = t(0.025)*s/sqrt(n)

to get the confidence interval as follows:

resampling_mean - E < population_mean < resampling_mean + E

In order to find the critical value of t, I need to have the degrees of freedom which would be my sample size - 1. I actually have a problem in understanding what the sample size is really in PSMC analysis.

I was wondering what is the common way of calculating CI for PSMC analysis?

Thank you, Homa