library(FrF2)

########################################################################################
# Full factorial

#-- 2 fattori

Runs<-2^2

Full2F<-FrF2(nruns=Runs, nfactors=2, randomize = TRUE, seed=7)

summary(Full2F)

#Per vedere tutto il contenuto dell'oggetto creato:

attributes(Full2F)

#Per estrarre solo la tabellina con il disegno sperimentale:

Disegno<-attributes(Full2F)$desnum

#Verifica Design Bilanciato

apply(Design,2,sum) # vuol dire: applica la funzione somma all'oggetto Design, per colonna

#Verifica Ortogonalità

sum(apply(Design,1,prod))

#---
#Se cambio il seed, cambia la sequenza degli esperimenti:

FrF2(nruns=Runs, nfactors=2, randomize = TRUE, seed=365)


########
#-- 3 fattori (completare spazi vuoti)

Runs<-2^3

Full3F<-FrF2(nruns=Runs, nfactors= 3, randomize = TRUE, seed=7)

summary(Full3F)


########
#-- 6 fattori (completare spazi vuoti)

Runs<-2^6

Full6F<-FrF2(nruns=Runs, nfactors=  6, randomize = TRUE, seed=7)

summary(Full6F)

########################################################################################
# Fractional factorial


#-- 3 fattori - fractional 1/2

Runs<-(2^3)/2

Fractional3F1.2<-FrF2(nruns=Runs, nfactors=3, randomize = TRUE, seed=7)

summary(Fractional3F1.2)

#Estrazione della tabellina con il disegno sperimentale:

DesignFr<-attributes(Fractional3F1.2)$desnum

#Verifica Design Bilanciato

apply(DesignFr,2,sum) # vuol dire: applica la funzione somma all'oggetto Design, per colonna

#Verifica Ortogonalità

sum(apply(DesignFr,1,prod))


#-- 6 fattori - fractional 1/2 (completare spazi vuoti)

Runs<-(2^6)/2 

Fractional6F1.2<-FrF2(nruns=Runs, nfactors= 6  , randomize = TRUE, seed=7)

summary(Fractional6F1.2)


#-- 6 fattori - fractional 1/4 (completare spazi vuoti)

Runs<-(2^6)/4 

Fractional6F1.4<-FrF2(nruns=Runs, nfactors= 6  , randomize = TRUE, seed=7)

summary(Fractional6F1.4)


#-- 6 fattori - fractional 1/8 (completare spazi vuoti)

Runs<-(2^6)/8

Fractional6F1.8<-FrF2(nruns=Runs, nfactors= 6 , randomize = TRUE, seed=7)

summary(Fractional6F1.8)

################################################################################
### INSTALLAZIONE MOLTO LUNGA!!!

library(agricolae)

# Per vedere l'elenco dei tipi di design sperimentale ottenibili:

designs <- ls("package:agricolae", pattern = "design")
designs