I'm trying to add a random perturbation term to my response variable (to estimate the impact of multicolinnearity in the beta estimates), but I am not sure on how to accomplish it when the response variable is transformed using a 1/4 power, in R.
Without any transformation in the response variable, according to Faraway's "Linear Models with R," the way to write it is:
g <- lm(hipcenter+10*rnorm(38) ~ ., seatpos).
So if I have a 1/4 power transformation for the response variable, then the code should be:
g <- lm((hipcenter^(1/4))+10*rnorm(38) ~ ., seatpos).
Is this right?