Diseño de factorial completo en parcelas divididas F1 = fenotipo (banano) G1, G2, G3 F2 = Fertilización aerea, terrestre rto = peso de racimos

df=data.frame(F1=gl(3,48,144,c("G1", "G2", "G3")),
           F2 =gl(2,24,144,c("Fa","Fs")),
           rto=rnorm(144,3,0.3))
df
##     F1 F2      rto
## 1   G1 Fa 2.669912
## 2   G1 Fa 3.140361
## 3   G1 Fa 3.518024
## 4   G1 Fa 3.485215
## 5   G1 Fa 3.340479
## 6   G1 Fa 2.969980
## 7   G1 Fa 2.892411
## 8   G1 Fa 3.149653
## 9   G1 Fa 2.690410
## 10  G1 Fa 2.904199
## 11  G1 Fa 2.857982
## 12  G1 Fa 2.800308
## 13  G1 Fa 2.892949
## 14  G1 Fa 3.272095
## 15  G1 Fa 3.443228
## 16  G1 Fa 2.689810
## 17  G1 Fa 3.123594
## 18  G1 Fa 2.978549
## 19  G1 Fa 3.217690
## 20  G1 Fa 2.728782
## 21  G1 Fa 3.076851
## 22  G1 Fa 3.107514
## 23  G1 Fa 3.185944
## 24  G1 Fa 2.247613
## 25  G1 Fs 2.835793
## 26  G1 Fs 3.182322
## 27  G1 Fs 3.209232
## 28  G1 Fs 3.139713
## 29  G1 Fs 3.292466
## 30  G1 Fs 3.502019
## 31  G1 Fs 3.275316
## 32  G1 Fs 2.661605
## 33  G1 Fs 2.838835
## 34  G1 Fs 2.863728
## 35  G1 Fs 2.920500
## 36  G1 Fs 3.117117
## 37  G1 Fs 2.775655
## 38  G1 Fs 3.618535
## 39  G1 Fs 2.829905
## 40  G1 Fs 3.163559
## 41  G1 Fs 3.471944
## 42  G1 Fs 3.655987
## 43  G1 Fs 2.910571
## 44  G1 Fs 2.643041
## 45  G1 Fs 3.420271
## 46  G1 Fs 2.455840
## 47  G1 Fs 3.103784
## 48  G1 Fs 3.110934
## 49  G2 Fa 2.911778
## 50  G2 Fa 3.058870
## 51  G2 Fa 2.916051
## 52  G2 Fa 2.953607
## 53  G2 Fa 2.594278
## 54  G2 Fa 3.068266
## 55  G2 Fa 2.934806
## 56  G2 Fa 3.231439
## 57  G2 Fa 3.051148
## 58  G2 Fa 2.563202
## 59  G2 Fa 2.464490
## 60  G2 Fa 2.506190
## 61  G2 Fa 2.892703
## 62  G2 Fa 2.913931
## 63  G2 Fa 2.839928
## 64  G2 Fa 2.587202
## 65  G2 Fa 2.969633
## 66  G2 Fa 3.207205
## 67  G2 Fa 3.367616
## 68  G2 Fa 3.059054
## 69  G2 Fa 2.989130
## 70  G2 Fa 3.142592
## 71  G2 Fa 3.113372
## 72  G2 Fa 2.841086
## 73  G2 Fs 2.807487
## 74  G2 Fs 3.259099
## 75  G2 Fs 2.672358
## 76  G2 Fs 2.878697
## 77  G2 Fs 3.137912
## 78  G2 Fs 2.761476
## 79  G2 Fs 3.049596
## 80  G2 Fs 3.270572
## 81  G2 Fs 2.808517
## 82  G2 Fs 3.369989
## 83  G2 Fs 2.773614
## 84  G2 Fs 2.904388
## 85  G2 Fs 2.813534
## 86  G2 Fs 3.271039
## 87  G2 Fs 2.721590
## 88  G2 Fs 2.875843
## 89  G2 Fs 2.604702
## 90  G2 Fs 3.742127
## 91  G2 Fs 3.053621
## 92  G2 Fs 2.907748
## 93  G2 Fs 2.845240
## 94  G2 Fs 2.767210
## 95  G2 Fs 3.610125
## 96  G2 Fs 2.647971
## 97  G3 Fa 2.910861
## 98  G3 Fa 2.491055
## 99  G3 Fa 3.025870
## 100 G3 Fa 2.825915
## 101 G3 Fa 2.447013
## 102 G3 Fa 2.896025
## 103 G3 Fa 3.069293
## 104 G3 Fa 2.671603
## 105 G3 Fa 3.114424
## 106 G3 Fa 3.074042
## 107 G3 Fa 3.273742
## 108 G3 Fa 3.460566
## 109 G3 Fa 2.133249
## 110 G3 Fa 2.650105
## 111 G3 Fa 2.966768
## 112 G3 Fa 2.887896
## 113 G3 Fa 3.058996
## 114 G3 Fa 3.003867
## 115 G3 Fa 3.040171
## 116 G3 Fa 3.304394
## 117 G3 Fa 3.326141
## 118 G3 Fa 3.588591
## 119 G3 Fa 3.152094
## 120 G3 Fa 3.051220
## 121 G3 Fs 2.750603
## 122 G3 Fs 3.086440
## 123 G3 Fs 2.707333
## 124 G3 Fs 2.822003
## 125 G3 Fs 2.389534
## 126 G3 Fs 2.635696
## 127 G3 Fs 2.597349
## 128 G3 Fs 2.618821
## 129 G3 Fs 3.047567
## 130 G3 Fs 3.520053
## 131 G3 Fs 3.159852
## 132 G3 Fs 3.367996
## 133 G3 Fs 2.458325
## 134 G3 Fs 2.778209
## 135 G3 Fs 2.687566
## 136 G3 Fs 3.235726
## 137 G3 Fs 3.117048
## 138 G3 Fs 3.247701
## 139 G3 Fs 2.736795
## 140 G3 Fs 3.312276
## 141 G3 Fs 2.710664
## 142 G3 Fs 3.374332
## 143 G3 Fs 3.144872
## 144 G3 Fs 2.586841

Arbol colapsable

library(collapsibleTree)
collapsibleTree(df,c("F1","F2","rto"))
xy=expand.grid(x=seq(0,110,10),
               y=0,110,10)

df$x=xy$x
df$y=xy$y

df$trt=interaction(df$F1,df$F2)
df$trt=sample(df$trt,144)
library(ggplot2)
ggplot(df,aes(x=x,y=y,fill=trt))+
  geom_tile()

Factorial incompleto en arreglo completamente al azar F1 = fertilización Orgánica F2 = Dosis de la F.O. resp = grados Brix

df = data.frame(brix = runif(24,5,10),
                F1 = gl(3,8,24,c("GALL", "POLL", "MIX")),
                F2 = gl(4,2,24,c("d0","d5","d10","d15")))
collapsibleTree(df,c("F1","F2","brix"))