A<-c(1,2)
B<-c(3,5,9)

 
joint<-c(A,B)

#Η συνάρτηση choose μας δίνει τους συνδιασμούς των 5 ανά 2
choose(5,2)
[1] 10

> choose(30,10)
[1] 30045015


library(gtools)
#Η συνάρτηση combinations μας δίνει όλους τους δυνδιασμούς από το διάνυσμα joint να επιλέξουμε 2 στοιχεία.
comb<-combinations(5, 2, v=joint, set=TRUE, repeats.allowed=FALSE)
      [,1] [,2]
 [1,]    1    2
 [2,]    1    3
 [3,]    1    5
 [4,]    1    9
 [5,]    2    3
 [6,]    2    5
 [7,]    2    9
 [8,]    3    5
 [9,]    3    9
[10,]    5    9


# με τη συνάρτηση apply εφαρμόζω σε κάθε γραμμή του πίκανα comb τη συνάρτηση της μέσης τιμής
apply(comb,1,mean)
 [1] 1.5 2.0 3.0 5.0 2.5 3.5 5.5 4.0 6.0 7.0


choose(40,20)
[1] 137846528820


# Η συνάρτηση setdiff αφαιρεί τα στοιχεία του ενός διανύσματος από το άλλο
comb3<-matrix(NA,nrow(comb),3)
for (i in 1:nrow(comb)){
comb3[i,]<-setdiff(joint,comb[i,])
}


apply(comb3,1,mean)

Tobs<-abs(apply(comb,1,mean)- apply(comb3,1,mean))

Tobs<-round(Tobs,2)