A Optimizing Pseudoexperiment Sample Size for Accurate Results

  • A
  • Thread starter Thread starter ChrisVer
  • Start date Start date
ChrisVer
Science Advisor
Messages
3,372
Reaction score
465
How can someone decide what's the best number of produced pseudoexperiments in his set up?
In particular I have a value N which I want to vary wrt 2 nuisance parameters with relative uncertainties \delta_1,\delta_2. I am producing n pseudoexperiments (samples) in each calculating the mean and the standard deviation of
N_i =N_i^0 \Big[1 + \delta_1 \mathcal{N}_1(0,1) + \delta_2 \mathcal{N}_2(0,1) \Big]
How can I decide whether the n-trials I am choosing is optimal?
I have reached the following conclusion after some thinking but I am not sure... the sample's relative uncertainty that is \sqrt{\text{Var}(N)}/\bar{N} should be as small as possible... is that a correct way?
 
Physics news on Phys.org
The answer is "as many as you can". There is no point where another pseudo-experiment makes things worse and not better.
 
Vanadium 50 said:
The answer is "as many as you can". There is no point where another pseudo-experiment makes things worse and not better.

In general I've seen plots where they show the "observed" let's say value (that they vary in each PE), and the sampled values for different/increasing n-trials (eg 100,1000,10000,1000000)... so, the sampled values get closer to the observed value, but also they get more statistically constrained... I was wondering if by looking at something like this allows someone to decide with what n they should go.
 
As much as your computers can reasonably compute. The uncertainty from the limited number of pseudoexperiments goes down, which is great.
 
Thread 'Why is there such a difference between the total cross-section data? (simulation vs. experiment)'
Well, I'm simulating a neutron-proton scattering phase shift. The equation that I solve numerically is the Phase function method and is $$ \frac{d}{dr}[\delta_{i+1}] = \frac{2\mu}{\hbar^2}\frac{V(r)}{k^2}\sin(kr + \delta_i)$$ ##\delta_i## is the phase shift for triplet and singlet state, ##\mu## is the reduced mass for neutron-proton, ##k=\sqrt{2\mu E_{cm}/\hbar^2}## is the wave number and ##V(r)## is the potential of interaction like Yukawa, Wood-Saxon, Square well potential, etc. I first...
Toponium is a hadron which is the bound state of a valance top quark and a valance antitop quark. Oversimplified presentations often state that top quarks don't form hadrons, because they decay to bottom quarks extremely rapidly after they are created, leaving no time to form a hadron. And, the vast majority of the time, this is true. But, the lifetime of a top quark is only an average lifetime. Sometimes it decays faster and sometimes it decays slower. In the highly improbable case that...
I'm following this paper by Kitaev on SL(2,R) representations and I'm having a problem in the normalization of the continuous eigenfunctions (eqs. (67)-(70)), which satisfy \langle f_s | f_{s'} \rangle = \int_{0}^{1} \frac{2}{(1-u)^2} f_s(u)^* f_{s'}(u) \, du. \tag{67} The singular contribution of the integral arises at the endpoint u=1 of the integral, and in the limit u \to 1, the function f_s(u) takes on the form f_s(u) \approx a_s (1-u)^{1/2 + i s} + a_s^* (1-u)^{1/2 - i s}. \tag{70}...

Similar threads

Replies
18
Views
3K
Replies
24
Views
5K
Replies
7
Views
2K
Replies
28
Views
6K
4
Replies
175
Views
25K
Replies
16
Views
6K
Back
Top