「被験者背景シミュレーション」の編集履歴（バックアップ）一覧はこちら

「被験者背景シミュレーション」(2007/02/01 (木) 21:36:05) の最新版変更点

追加された行は緑色になります。

削除された行は赤色になります。

 # n = 19 (nonsmoker)
 # n = 49 (former smoker)
 # n = 32 (smoker)
 smoke.num <- c(rep(1, 19), rep(2, 49), rep(3, 32))
 label.smoke <- c("Nonsmoker","Former smoker","Smoker")
 smoke <- factor(smoke.num, levels=1:3, labels=label.smoke)
 
 mu <- mean(log(smoke.num))
 sg <- sqrt(var(log(smoke.num)))
 
 # Empirical distribution
 Pi <- cumsum(as.numeric(table(smoke)) / length(smoke))
 Pi <- Pi[-length(Pi)]
 
 # Continuous critical value
 CrV <- exp(mu + sg * qnorm(Pi))
 
 # Simulation
 set.seed(27011)
 nsim <- 1000
 smoke.sim.continuous <- exp(rnorm(nsim, mean=mu, sd=sg))
 smoke.sim <- cut(smoke.sim.continuous, breaks=c(0,CrV,max(smoke.sim.continuous)), labels=label.smoke)
 
 table(smoke.sim)
 
 # Plot
 hist(smoke.sim.continuous, freq=F, xlab='Smoke', main='')
 lines(density(smoke.sim.continuous))
 abline(v=CrV, col=4)