Using bootstrapped regression to test if categorical and/or count data has an effect on the outcome of another variable

[niamhailbhe]

I want to test the effect of Habitat (categorical), caterpillar SR, caterpillar abundance, plant SR and plant abundance on Z. filipendulae (count).

I have decided to use bootstrapped regression as the assumptions of ANCOVA have not been met.

R code below

bootReg<-function(formula=Z. filipendulae~Catepillar.SR+Caterpillar.Abundance+Plant.SR+Plant.Abundance+factor(Habitat), data=AncovaData, i)
{
d <- data[i,]
fit <- lm(formula, data = d)
return(coef(fit))
}

bootResults<-boot(statistic =bootReg,
formula=Z.filipendulae~Caterpillar.SR +
Caterpillar.Abundance+Plant.SR+Plant.Abundance+factor(Habitat),
data=Data,R=2000)

Error

Error in t.star[r, ] <- res[[r]] :
number of items to replace is not a multiple of replacement length

[andresrcs]

I think you have forgot to tell us what your problem is, you are just telling us what have you done so far but not what problems you are having.

Ideally, you should prepare a reproducible example (reprex) illustrating your issue.
Please have a look at this guide, to see how to create one:

FAQ: How to do a minimal reproducible example ( reprex ) for beginners

A minimal reproducible example consists of the following items: A minimal dataset, necessary to reproduce the issue The minimal runnable code necessary to reproduce the issue, which can be run on the given dataset, and including the necessary information on the used packages. Let's quickly go over each one of these with examples: Minimal Dataset (Sample Data) You need to provide a data frame that is small enough to be (reasonably) pasted on a post, but big enough to reproduce your issue. Let's say, as an example, that you are working with the iris data frame head(iris) #> Sepal.Length Sepal.Width Petal.Length Petal.Width Species #> 1 5.1 3.5 1.4 0.…