dput(head(reg,20))

structure(list(sample = c(13L, 14L, 15L, 13L, 14L, 15L, 13L,
14L, 15L, 13L, 14L, 15L, 13L, 14L, 15L, 13L, 14L, 15L, 1L, 2L
), days = c(7L, 7L, 7L, 14L, 14L, 14L, 21L, 21L, 21L, 28L, 28L,
28L, 35L, 35L, 35L, 42L, 42L, 42L, 7L, 7L), treat = structure(c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L), .Label = c("asynchron", "control", "synchron"), class = "factor"),
repl = c(1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L,
1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L), disp = structure(c(2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L), .Label = c("high", "low", "medium"), class = "factor"),
gbiov = c(205088809L, 236043290L, 182258079L, 356710654L,
519036032L, 423556500L, 265916788L, 260495312L, 424304021L,
456717288L, 412530454L, 495216491L, 360012083L, 246070292L,
351623075L, 99075488L, 141565279L, 109303015L, 157719021L,
275901505L), rich = c(15L, 15L, 15L, 12L, 12L, 10L, 10L,
10L, 10L, 8L, 7L, 9L, 9L, 9L, 8L, 6L, 8L, 8L, 15L, 15L),
evenness = c(0.810893754, 0.766130198, 0.82401103, 0.422734873,
0.365797654, 0.406629574, 0.317890623, 0.434262493, 0.475385016,
0.311819454, 0.307131166, 0.321269955, 0.365280751, 0.505205296,
0.43643584, 0.484220445, 0.455471176, 0.423730082, 0.811921742,
0.833306048)), .Names = c("sample", "days", "treat", "repl",
"disp", "gbiov", "rich", "evenness"), row.names = c(NA, 20L), class = "data.frame")

sample   days  treat    repl   disp   gbiov richness evenness

test<-ezANOVA(data=reg, dv=.(gbiov), wid=.(repl), within=.(treat, disp, days))

Warnung: Converting "repl" to factor for ANOVA.
Warnung: "days" will be treated as numeric.
Fehler in ezANOVA_main(data = data, dv = dv, wid = wid, within = within,  :
One or more cells is missing data. Try using ezDesign() to check your data.