I'm working with the Iris data.One of the variables,PetalWidth,Has a clear bimodal distribution

My understanding is that Multivariate regression Assumes normality for each of the input variables

Can I continue With the variable left alone? Is it necessary To normalize it?

Here's the r code to Access the data

dni3 <- dimnames(iris3)
ii <- data.frame(matrix(aperm(iris3, c(1,3,2)), ncol = 4,
                        dimnames = list(NULL, sub(" L.",".Length",
                                                  sub(" W.",".Width", dni3[[2]])))),
                 Species = gl(3, 50, labels = sub("S", "s", sub("V", "v", dni3[[3]]))))
all.equal(ii, iris) # TRUE

I assume you mean multiple regression?

Anyway, it doesn't assume normality of the input variables. It assumes the response variable is conditionally distributed Gaussian (normal) but doesn't assume anything about the covariates or predictor variables (that said, transforming the covariates so that it's not just a few extreme values dominating the estimated effect often makes sense.)

Let's assume you are modelling petal width and it is bimodal. We might assume that the width of the petal is in some way connected to some of the petals coming from species A and some from species B. Once we account for the effect of species, the bimodality disappears if it was due to species as we essentially subtract each species mean from the data, which moves the two modes of the distribution together to be approximately 0. We then can say that conditional upon the model the response is distributed Gaussian (or rather we hope we can say this, but this is an assumption we need to check). Or in linear regression, that the residuals, which are the data after we have account for the model, are distributed Gaussian.

For example:


mod <- lm(Petal.Width ~ Species, data = iris)
iris <- cbind(iris, Residual = resid(mod))

p1 <- ggplot(iris, aes(x = Petal.Width)) +
  geom_histogram(bins = 10) + labs(title = "Original data")
p2 <- ggplot(iris, aes(x = Residual))+
  geom_histogram(bins = 10) + labs(title = "Residuals")
plot_grid(p1, p2, ncol = 2, align = 'hv', axis = 'lrtb')

Which produces

enter image description here

The residuals aren't really Gaussian, but this is not the right model for a variable constrained at 0 (can't observe zero width or negative width petals). However, the bimodality has gone because it was due to (careful with causality) or related to differences between species.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.