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.793408
## 2   G1 Fa 3.052887
## 3   G1 Fa 3.266410
## 4   G1 Fa 3.169346
## 5   G1 Fa 2.872325
## 6   G1 Fa 2.953033
## 7   G1 Fa 2.732439
## 8   G1 Fa 2.863219
## 9   G1 Fa 2.777013
## 10  G1 Fa 3.660082
## 11  G1 Fa 3.287176
## 12  G1 Fa 3.073834
## 13  G1 Fa 3.171644
## 14  G1 Fa 3.248184
## 15  G1 Fa 2.995427
## 16  G1 Fa 3.322549
## 17  G1 Fa 3.218643
## 18  G1 Fa 3.365239
## 19  G1 Fa 2.540294
## 20  G1 Fa 2.797180
## 21  G1 Fa 2.697166
## 22  G1 Fa 3.040030
## 23  G1 Fa 2.883869
## 24  G1 Fa 3.011719
## 25  G1 Fs 3.075270
## 26  G1 Fs 2.871830
## 27  G1 Fs 2.926142
## 28  G1 Fs 3.323535
## 29  G1 Fs 2.567691
## 30  G1 Fs 3.216614
## 31  G1 Fs 2.955244
## 32  G1 Fs 3.035647
## 33  G1 Fs 2.876824
## 34  G1 Fs 2.361977
## 35  G1 Fs 3.540452
## 36  G1 Fs 2.934011
## 37  G1 Fs 2.907562
## 38  G1 Fs 2.815502
## 39  G1 Fs 2.554373
## 40  G1 Fs 2.825487
## 41  G1 Fs 2.883684
## 42  G1 Fs 3.486788
## 43  G1 Fs 3.270348
## 44  G1 Fs 3.107397
## 45  G1 Fs 3.390695
## 46  G1 Fs 2.825622
## 47  G1 Fs 2.822158
## 48  G1 Fs 2.907584
## 49  G2 Fa 3.043960
## 50  G2 Fa 2.816733
## 51  G2 Fa 3.487482
## 52  G2 Fa 2.988910
## 53  G2 Fa 3.067393
## 54  G2 Fa 2.810152
## 55  G2 Fa 3.532661
## 56  G2 Fa 2.970534
## 57  G2 Fa 3.039914
## 58  G2 Fa 2.926793
## 59  G2 Fa 2.895908
## 60  G2 Fa 2.608305
## 61  G2 Fa 2.136715
## 62  G2 Fa 2.848102
## 63  G2 Fa 2.938107
## 64  G2 Fa 2.951347
## 65  G2 Fa 2.934597
## 66  G2 Fa 2.720339
## 67  G2 Fa 3.379810
## 68  G2 Fa 3.153469
## 69  G2 Fa 2.986627
## 70  G2 Fa 2.870662
## 71  G2 Fa 2.923054
## 72  G2 Fa 3.324433
## 73  G2 Fs 3.248232
## 74  G2 Fs 2.653636
## 75  G2 Fs 2.915050
## 76  G2 Fs 2.846837
## 77  G2 Fs 3.119420
## 78  G2 Fs 2.656760
## 79  G2 Fs 2.602913
## 80  G2 Fs 3.260757
## 81  G2 Fs 2.805084
## 82  G2 Fs 3.296860
## 83  G2 Fs 3.072142
## 84  G2 Fs 2.898098
## 85  G2 Fs 2.348458
## 86  G2 Fs 3.307880
## 87  G2 Fs 3.126558
## 88  G2 Fs 3.110467
## 89  G2 Fs 3.458595
## 90  G2 Fs 3.173943
## 91  G2 Fs 3.015438
## 92  G2 Fs 3.041120
## 93  G2 Fs 2.863140
## 94  G2 Fs 3.225992
## 95  G2 Fs 2.992630
## 96  G2 Fs 2.793926
## 97  G3 Fa 2.645237
## 98  G3 Fa 2.975815
## 99  G3 Fa 2.302255
## 100 G3 Fa 3.211765
## 101 G3 Fa 3.084633
## 102 G3 Fa 2.425979
## 103 G3 Fa 2.750995
## 104 G3 Fa 2.769210
## 105 G3 Fa 3.145503
## 106 G3 Fa 2.676121
## 107 G3 Fa 3.130636
## 108 G3 Fa 2.812714
## 109 G3 Fa 2.434484
## 110 G3 Fa 3.033142
## 111 G3 Fa 3.270610
## 112 G3 Fa 3.050890
## 113 G3 Fa 2.728008
## 114 G3 Fa 3.511454
## 115 G3 Fa 3.148624
## 116 G3 Fa 2.965105
## 117 G3 Fa 2.783856
## 118 G3 Fa 2.984447
## 119 G3 Fa 2.775715
## 120 G3 Fa 2.870218
## 121 G3 Fs 4.024147
## 122 G3 Fs 3.665600
## 123 G3 Fs 3.333872
## 124 G3 Fs 2.493231
## 125 G3 Fs 3.310431
## 126 G3 Fs 3.208816
## 127 G3 Fs 2.882196
## 128 G3 Fs 3.392514
## 129 G3 Fs 2.976178
## 130 G3 Fs 2.961146
## 131 G3 Fs 2.648493
## 132 G3 Fs 2.738021
## 133 G3 Fs 2.542976
## 134 G3 Fs 3.241982
## 135 G3 Fs 2.919360
## 136 G3 Fs 3.214831
## 137 G3 Fs 2.392000
## 138 G3 Fs 3.134104
## 139 G3 Fs 3.207237
## 140 G3 Fs 2.629022
## 141 G3 Fs 2.956388
## 142 G3 Fs 2.953316
## 143 G3 Fs 3.016866
## 144 G3 Fs 2.727470

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)
## Warning: package 'ggplot2' was built under R version 4.2.2
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"))