Backward RDA model

scope <- list(lower = ~1, upper = ~Si+Al+ K+Ti+Mg+ED+Mn+Sr+Zr+C.N+C.P+S+pH)

ordiR2step(my_RDA,

  •          scope = scope,
    
  •          Pin = 0.05,
    
  •          permutations = how(nperm = 499),
    
  •          trace = TRUE,
    
  •          R2permutations = 1000)
    

Error in reformulate(scope) :
'termlabels' must be a character vector of length at least one
Can someone help me?