Hi,
>>A = rep(rep(c(0,1),4),2)
>>B = rep(rep(c(0,1),each=4),2)
>>C = rep(rep(rep(c(0,1),2),each=2),2)
>>X = data.frame(A,B,C,rnorm(16))
>>
>>
> X = data.frame(A,B,C,D=rnorm(16))
> X
A B C D
1 0 0 0 0.829987
2 1 0 0 1.319203
3 0 0 1 0.223752
4 1 0 1 0.017353
5 0 1 0 1.288284
6 1 1 0 1.384687
7 0 1 1 0.186565
8 1 1 1 -0.002799
9 0 0 0 -0.875205
10 1 0 0 1.872563
11 0 0 1 -0.317170
12 1 0 1 -0.766612
13 0 1 0 -1.423495
14 1 1 0 -0.042385
15 0 1 1 -1.191448
16 1 1 1 0.731018
> Q<-unique(X[1:3])
> Q
A B C
1 0 0 0
2 1 0 0
3 0 0 1
4 1 0 1
5 0 1 0
6 1 1 0
7 0 1 1
8 1 1 1
> for(i in 1:nrow(Q)){
+ print(X[apply(X[,1:3],1,function(x){all(x==Q[i,])}),])
+ }
A B C D
1 0 0 0 0.8300
9 0 0 0 -0.8752
A B C D
2 1 0 0 1.319
10 1 0 0 1.873
A B C D
3 0 0 1 0.2238
11 0 0 1 -0.3172
A B C D
4 1 0 1 0.01735
12 1 0 1 -0.76661
A B C D
5 0 1 0 1.288
13 0 1 0 -1.423
A B C D
6 1 1 0 1.38469
14 1 1 0 -0.04238
A B C D
7 0 1 1 0.1866
15 0 1 1 -1.1914
A B C D
8 1 1 1 -0.002799
16 1 1 1 0.731018
>
HTH
--
Ferdinand Alimadhi
Programmer / Analyst
Harvard University
The Institute for Quantitative Social Science
(617) 496-0187
[EMAIL PROTECTED]
www.iq.harvard.edu
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
