確率が2つの分布から描画するR

Question

確率が100000倍の2つの異なる分布から描画しようとしています。残念ながら私はforループで何が間違っているかはわかりませんが、希望する100,000の値の代わりに1つの値をsimulated_dataに加えるだけです。

質問1：これを修正するにはどうすればよいですか？

質問2：リスト内の100,000項目をループする必要がない、はるかに効率的な方法がありますか？ここで

#creating a vector of probabilities 
probabilities <- rep(0.99,100000) 
#creating a vector of booleans 
logicals <- runif(length(probabilities)) < probabilities 

#empty list for my simulated data 
simulated_data <- c() 

#drawing from two different distributions depending on the value in logicals 
for(i in logicals){ 

    if (isTRUE(i)) { 
    simulated_data[i] <- rnorm(n = 1, mean = 0, sd = 1) 
    }else{ 
    simulated_data[i] <- rnorm(n = 1, mean = 0, sd = 10) 
    } 
} 

Answer 1

ここに誰のための素晴らしいソリューションです：

n <- 100000 
prob1 <- 0.99 
prob2 <- 1-prob1 

dist1 <- rnorm(prob1*n, 0, 1) 
dist2 <- rnorm(prob2*n, 0, 10) 

actual_sample <- c(dist1, dist2) 

Answer 2

あなたが確率0.99と0.01で、各要素はSAMPLE1またはsample2のいずれかからランダムに取られ、最終的なサンプルを作成したいと思われます。

正しいアプローチは、それぞれが同じ数の要素を含み、どちらかからランダムに選択する両方のサンプルを生成することです。

正しいアプローチは次のようになります

# Generate both samples 
n = 100000 
sample1 = rnorm(n,0,1) 
sample2 = rnorm(n,0,10) 

# Create the logical vector that will decide whether to take from sample 1 or 2 
s1_s2 = runif(n) < 0.99 

# Create the final sample 
sample = ifelse(s1_s2 , sample1, sample2) 

この場合には、* N SAMPLE2からSAMPLE1と0.01から正確に150個の* n個のサンプルがあることを保証されません。実際：

これは、予想通り0.99 * nに近いですが、正確ではありません。

Answer 3

は、各分布からの値の所望の画分を有するベクターを作成し、値のランダム順列を作成する：

N = 10000 
frac =0.99 
rand_mix = sample(c(rnorm(frac*N, 0, sd=1) , rnorm((1-frac)*N, 0, sd=10))) 

> table(abs(rand_mix) >1.96) 

FALSE TRUE 
9364 636 
> (100000-636)/100000 
[1] 0.99364 

> table(rnorm(10000) >6) 

FALSE 
10000 

割合が固定されています。あなたはおそらくランダム画分（しかし、統計的に150に近い）をwante場合、この方法を試してください。

> table(sample(c(rnorm(10e6), rnorm(10e4, sd=10)), 10e4) > 1.96) 

FALSE TRUE 
97151 2849 

は、と比較する：

> N = 100000 
> frac =0.99 
> rand_mix = sample(c(rnorm(frac*N, 0, sd=1) , rnorm((1-frac)*N, 0, sd=10))) 
> table(rand_mix > 1.96) 

FALSE TRUE 
97117 2883