Monte Carlo simulation, resampling

目的:

原因:工具能做什么

是什么:Computer simulation that generates are large number of simulated samples of data based on an assumed Data Generating Process (DGP) that characterizes the population from which the simulated samples are drawn.

  • Patterns in those simulated samples are then summarized  and described.
  • Such patterns can be evaluated in terms of substantive theory or in terms of the statistical properties of some

estimator

是什么:Data  Generating Process (DGP)

  • A DGP describes how a values of a variable of interest are  produced in the population.
  • Most DGP’s of interest include a systematic component and a stochastic component.
  • We use statistical analysis to infer characteristics of the DGP by analyzing observable data sampled from the population.
  • In applied statistical work, we never know the DGP – if we did, we wouldn’t need statistical estimates of it.
  • In Monte Carol simulations, we do know the DGP because we create it.


是什么:Resampling

  • Like Monte Carlo simulations, resampling methods use a  computer to generate a large number of simulated samples of data.
  • Also like Monte Carlo simulations, patterns in these simulated samples are then summarized, and the results used

to evaluate substantive theory or statistical estimators.

  • What is different is that the simulated samples are generated by drawning new samples (with replacement) from the sample of data you have
  • In resampling methods, the researcher DOES NOT know or control the DGP, but the goal of learning about the DGP remains the same.

Monte Carlo Simulation of OLS

Know Your Assumptions:




set.seed(123456) # Set the seed for reproducible results
sims= 500 # Set the number of simulations at the top of the script
alpha.1 = numeric(sims) # Empty vector for storing the simulated intercepts
B.1 = numeric(sims) # Empty vector for storing the simulated slopes
a = .2 # True value for the intercept
b =.5 # True value for the slope
n = 1000 # sample size
X = runif(n, -1, 1) # Create a sample of n observations on the variable X.
# Note that this variable is outside the loop, because X
# should be fixed in repeated samples.
for(i in 1:sims)– # Start the loop
Y = a + b*X + rnorm(n, 0, 1) # The true DGP, with N(0, 1) error
model = lm(Y ˜ X) # Estimate OLS Model

alpha.1[i] = model$coef[1] # Put the estimate for the intercept
# in the vector alpha.1
B.1[i] = model$coef[2] # Put the estimate for X in the vector B.1

˝ # End loop



  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Monte Carlo simulation is a powerful tool for modeling complex physical systems. It involves using random numbers to simulate the behavior of a system over time, and can be used to study a wide range of phenomena, from the behavior of particles in a collider to the spread of disease in a population. In the context of the reconstruction toolkit, Monte Carlo simulation can be used to model the behavior of a detector system and the particles that interact with it. This allows researchers to study the performance of the detector under different conditions and optimize its design. To carry out a Monte Carlo simulation with the reconstruction toolkit, several steps are typically involved: 1. Define the geometry of the detector system, including the materials used and the placement of sensors and other components. 2. Specify the particles to be simulated, including their energy, momentum, and other properties. 3. Use the Monte Carlo algorithm to simulate the behavior of the particles as they interact with the detector system. 4. Collect data on the resulting interactions and use this data to reconstruct the original events. 5. Analyze the reconstructed events to study the performance of the detector and make improvements as needed. Overall, Monte Carlo simulation with the reconstruction toolkit is a powerful tool for studying complex physical systems and optimizing their design. It requires expertise in both physics and computer science, but can lead to important insights and breakthroughs in a wide range of fields.
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值