不平衡数据的设计矩阵答案

【问题标题】：Design matrix for unbalanced data不平衡数据的设计矩阵
【发布时间】：2011-12-22 01:16:44
【问题描述】：

我正在尝试提出一个不平衡的设计矩阵/虚拟变量矩阵（在模拟项目的背景下）。目前模拟的不等观测值如下：

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 循环的简单方法？不用说，我的真实数据集有 1000 次观察。非常感谢您对此的任何帮助。

【问题讨论】：

标签： r simulation

【解决方案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

【讨论】：

谢谢。那行得通。我只是将它修改为更通用一点，因为 n_i 非常大。