R：この条件付き2項確率変数を生成するには？

Question

次のモデルからXiを生成します。私はbinomial distributionからデータを生成する方法を知っていますが、私ができることは以下の通りですが、同時にデータを生成するにはどうすればY = 1またはY =-1の条件を満たすことができますか？

set.seed(100) 
n1<-50 
n2<-50 
p1<- 0.4 
p2<- 0.3 

y <- c(rep(1,50),rep(-1,50)) 
xi1<- c(rbinom(n1,1,p1), rbinom(n2,1,p2)) 

Answer 1

rbinom()はベクトル化されているので、我々はあなたがよりよい密接?rbinomを読んだ

n <- ifelse(y == 1, n1, n2) 
p <- ifelse(y == 1, p1, p2) 
rbinom(n, 1, p) 

を行うことができます。