I am having a problem with my code. I am trying to create a GLM model with Zeroinfl since my data set has many 0’s in it. I am trying to study the influence of pond characteristics on pond amphibian diversity but I get an error code that leads to “NA”‘s in my summary table.

i have a df that contains data for the nation and would like to see if the relationship between variables is different at smaller geographies (region, state).
i have tried using the subset argument of glm but this would be very repetitive as i have many regions and many dependent variables.

glm.nb() automatically calculates an init.theta when it is not user specified.

model <- lm(y ~ ., data = A)
Error in model.frame.default(formula = y ~ ., data = A, drop.unused.levels = TRUE) :
variable lengths differ (found for ‘BAP’)