不平衡データ用の設計行列

Question

私は、（シミュレーションプロジェクトのコンテキストで）アンバランスな設計行列/ダミー変数行列を考案しようとしています。次のように不均等な観測が現在、シミュレートされています

set.seed(1) 
n_i <- as.matrix(round(runif(3,min=0, max=5))) 
n_i 
     [,1] 
[1,] 1 
[2,] 2 
[3,] 3 

私は上記のシミュレートされたサンプルサイズに基づいて、次の形式の行列を作成したいと思います：

des <- cbind(c(1,rep(0,5)),c(0,rep(1,2),rep(0,3)), c(rep(0,3),rep(1,3))) 

は、これを行うのが簡単な方法です複雑なfor-loopsは必要ありませんか？言うまでもなく、私の実際のデータセットは、1000年の観測データに基づいています。 これについて助けていただければ幸いです。

Answer 1

この1

Trt <- factor(x=rep(x=c(1, 2, 3), times=c(1, 2, 3)), labels=c("Trt1", "Trt2", "Trt3")) 
model.matrix(~Trt-1, data=Trt) 



    TrtTrt1 TrtTrt2 TrtTrt3 
1  1  0  0 
2  0  1  0 
3  0  1  0 
4  0  0  1 
5  0  0  1 
6  0  0  1