модель ln(y) = ln(b0) + b1*ln(x). Получить этот ln(b0) я, насколько я понимаю, могу вот так:
dat <- data.frame(x = seq(1,10, 1), y = rnorm(10, mean = 1000))
sl <- summary(lm(log(y) ~ log(x), data = dat))
log(sl$coefficients[[1]])
Внимание, вопрос: остаются ли все остальные данные в модели верными? Ведь нет? И если нет, то можно ли как-то указать трансформацию коэффициента внутри lm()?
Что где когда прям
Обсуждают сегодня