2011-12-22 2 views
1

시뮬레이션 프로젝트의 컨텍스트에서 불균형 한 디자인 행렬/더미 변수 행렬을 만들려고합니다. 다음과 같이 불평등 관측은 현재 시뮬레이션 :불균형 데이터를위한 설계 매트릭스

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 회의 관측 결과에 포함됩니다. 이 문제에 도움을 주시면 감사하겠습니다.

답변

2

이 하나

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 
+0

감사를 사용해보십시오. 그게 효과가 있었어. n_i가 매우 큰만큼 좀 더 일반적인 것으로 수정했습니다. – user969524