‘|’ not meaningful for factors - running ANCOVA in WRS2 - ANCOVA works but NOT with random factors!

Hi!

I haven't been able to run an ANCOVA with random effects.

I get error message "'|' not meaningful for factors"

Please help!

Thank you!

.
.
.

my code

anco2 <-  ancova(Above_Biomass ~                   #attempting to add random effects
                      Proportion_myco_col +
                      Proportion_myco_col*InvasionTrt +
                      InvasionTrt,      
                      data = short) 

ancorand <-  ancova(Above_Biomass ~                   #attempting to add random effects
                      Proportion_myco_col +
                      Proportion_myco_col*InvasionTrt +
                      InvasionTrt +
                      (1|DF_type) +
                      (1|SiteID) +
                      (1|block),       
                    data = short)

MRE

structure(list(block = c(1L, 1L, 2L, 1L, 1L, 1L, 1L, 3L, 3L,
2L, 1L, 2L, 1L, 1L, 3L, 2L, 2L, 3L, 2L, 1L, 2L, 3L, 2L, 2L, 1L,
2L, 2L, 1L, 1L, 1L, 3L, 3L, 1L, 3L, 2L, 2L, 2L, 1L, 2L, 3L, 3L,
2L, 3L, 3L, 1L, 3L, 3L, 1L, 3L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L),
SiteID = c(1L, 2L, 2L, 2L, 1L, 2L, 3L, 3L, 2L, 1L, 1L, 2L,
1L, 3L, 1L, 2L, 3L, 2L, 3L, 3L, 2L, 1L, 1L, 2L, 3L, 1L, 2L,
2L, 2L, 1L, 3L, 2L, 3L, 1L, 3L, 3L, 3L, 2L, 3L, 3L, 2L, 1L,
1L, 3L, 1L, 1L, 3L, 3L, 2L, 2L, 3L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 1L, 1L, 1L,
1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L), DF_type = structure(c(1L, 2L, 2L, 3L,
3L, 2L, 3L, 3L, 2L, 1L, 2L, 3L, 2L, 2L, 3L, 3L, 1L, 3L, 2L,
1L, 3L, 2L, 3L, 1L, 2L, 1L, 2L, 3L, 1L, 1L, 1L, 2L, 3L, 3L,
3L, 2L, 3L, 1L, 1L, 3L, 1L, 2L, 1L, 2L, 3L, 2L, 2L, 1L, 3L,
1L, 1L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L,
1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 2L, 3L, 1L, 2L, 3L, 1L, 2L,
3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 3L, 1L, 2L, 3L, 1L, 3L, 1L,
3L, 1L, 3L, 1L, 3L, 1L, 2L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L,
3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L,
3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L,
3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L), .Label = c("primadonna",
"street_fighter_1", "street_fighter_2"), class = "factor"),
InvasionTrt = structure(c(1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L,
1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L,
2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L,
2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L,
2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L,
1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L,
2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L,
2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L,
1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L), .Label = c("broom",
"uninvaded"), class = "factor"), Above_Biomass = c(27.1,
78.2, 179, 70.7, 106.8, 176.4, 122.4, NA, 124.8, 85.6, 87.5,
70, 120.4, 381.7, 41.8, 54.1, 46.2, 138, 59.5, 107.6, 208.1,
62.5, 86.9, 46.3, 75.3, 276.6, 94.8, 78.1, 172.3, 81.5, 66.8,
102.5, 37.8, 472.4, 56.8, 317.8, 115.8, 118.1, 230.2, 149.1,
286.5, 150, NA, NA, 59.7, 235.6, NA, 139, 232.8, 120.9, NA,
NA, NA, 42.5, 59.5, 52, NA, NA, NA, 27.9, NA, 82.9, NA, NA,
NA, 68.3, 64.4, 90.2, NA, NA, NA, 66.7, NA, NA, NA, 89.6,
62.1, 121.8, NA, NA, NA, 52.1, 22.3, 65.6, NA, NA, NA, NA,
NA, 8.9, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 53.8, 75.8,
NA, NA, NA, 67.3, 46.6, 127.7, NA, NA, NA, 247.5, 201.1,
NA, NA, NA, NA, 95.65, 53.6, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA), Proportion_myco_col = c(0.027777778,
0.027777778, 0.041666667, 0.055555556, 0.083333333, 0.090909091,
0.097222222, 0.097222222, 0.152777778, 0.166666667, 0.171875,
0.180555556, 0.180555556, 0.180555556, 0.194444444, 0.208333333,
0.263888889, 0.263888889, 0.277777778, 0.333333333, 0.333333333,
0.361111111, 0.388888889, 0.416666667, 0.416666667, 0.416666667,
0.416666667, 0.424242424, 0.472222222, 0.486111111, 0.486111111,
0.5, 0.513888889, 0.513888889, 0.527777778, 0.527777778,
0.541666667, 0.555555556, 0.569444444, 0.583333333, 0.611111111,
0.638888889, 0.652777778, 0.652777778, 0.71875, 0.736111111,
0.75, 0.777777778, 0.819444444, 0.833333333, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, 0.107142857, NA, NA, NA, 0.111111111, 0.222222222,
0.625, NA, NA, NA, 0.736111111, 0.708333333, NA, NA, NA,
NA, 0.416666667, 0.341463415, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA), block1 = c(1,
1, 2, 1, 1, 1, 1, 3, 3, 2, 1, 2, 1, 1, 3, 2, 2, 3, 2, 1,
2, 3, 2, 2, 1, 2, 2, 1, 1, 1, 3, 3, 1, 3, 2, 2, 2, 1, 2,
3, 3, 2, 3, 3, 1, 3, 3, 1, 3, 2, 1, 2, 2, 2, 2, 2, 2, 2,
2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
2, 2, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5), site = c(1,
2, 2, 2, 1, 2, 3, 3, 2, 1, 1, 2, 1, 3, 1, 2, 3, 2, 3, 3,
2, 1, 1, 2, 3, 1, 2, 2, 2, 1, 3, 2, 3, 1, 3, 3, 3, 2, 3,
3, 2, 1, 1, 3, 1, 1, 3, 3, 2, 2, 3, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3,
3, 3, 3, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3,
3, 3, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3,
3, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2,
2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3), geno = c(1,
2, 2, 3, 3, 2, 3, 3, 2, 1, 2, 3, 2, 2, 3, 3, 1, 3, 2, 1,
3, 2, 3, 1, 2, 1, 2, 3, 1, 1, 1, 2, 3, 3, 3, 2, 3, 1, 1,
3, 1, 2, 1, 2, 3, 2, 2, 1, 3, 1, 1, 3, 1, 2, 3, 1, 2, 3,
1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 2, 3, 1, 2, 3,
1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 3, 1, 2, 3, 1, 3, 1, 3, 1,
3, 1, 3, 1, 2, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2,
3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3,
1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3), inv = structure(c(1L,
1L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L,
2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L,
1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L,
2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,
1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L,
2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L,
1L), .Label = c("broom", "uninvaded"), class = "factor"),
myco = c(0.027777778, 0.027777778, 0.041666667, 0.055555556,
0.083333333, 0.090909091, 0.097222222, 0.097222222, 0.152777778,
0.166666667, 0.171875, 0.180555556, 0.180555556, 0.180555556,
0.194444444, 0.208333333, 0.263888889, 0.263888889, 0.277777778,
0.333333333, 0.333333333, 0.361111111, 0.388888889, 0.416666667,
0.416666667, 0.416666667, 0.416666667, 0.424242424, 0.472222222,
0.486111111, 0.486111111, 0.5, 0.513888889, 0.513888889,
0.527777778, 0.527777778, 0.541666667, 0.555555556, 0.569444444,
0.583333333, 0.611111111, 0.638888889, 0.652777778, 0.652777778,
0.71875, 0.736111111, 0.75, 0.777777778, 0.819444444, 0.833333333,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, 0.107142857, NA, NA, NA, 0.111111111,
0.222222222, 0.625, NA, NA, NA, 0.736111111, 0.708333333,
NA, NA, NA, NA, 0.416666667, 0.341463415, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA)), row.names = c(1L,
2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L,
16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L,
29L, 30L, 31L, 32L, 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L,
42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 100L, 101L, 102L,
103L, 104L, 105L, 106L, 107L, 108L, 109L, 110L, 111L, 112L, 113L,
114L, 115L, 116L, 117L, 118L, 119L, 120L, 121L, 122L, 123L, 124L,
125L, 126L, 127L, 128L, 129L, 130L, 131L, 132L, 133L, 134L, 135L,
136L, 137L, 138L, 139L, 140L, 141L, 142L, 143L, 144L, 145L, 146L,
147L, 148L, 149L, 150L, 250L, 251L, 252L, 253L, 254L, 255L, 256L,
257L, 258L, 259L, 260L, 261L, 262L, 263L, 264L, 265L, 266L, 267L,
268L, 269L, 270L, 271L, 272L, 273L, 274L, 275L, 276L, 277L, 278L,
279L, 280L, 281L, 282L, 283L, 284L, 285L, 286L, 287L, 288L, 289L,
290L, 291L, 292L, 293L, 294L, 295L, 296L, 297L, 298L, 299L, 300L
), class = "data.frame")

This topic was automatically closed 21 days after the last reply. New replies are no longer allowed.