В моем моделировании независимая переменная - это последовательность, которая равномерно распространяется от 0 до 1. Код ниже показывает, как генерируются x, ошибка u и зависимая переменная y.
Для моей вспомогательной регрессии FGLS мне нужно взять лог x. Однако взятие логарифма 0 дает бесконечность в R, s.t. Я получаю сообщение об ошибке Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : NA/NaN/Inf in 'x' при попытке вычислить вспомогательную регрессию.
x <- seq(from = 0, to = 1, by = 0.01)
N <- length(x)
u <- rnorm(N, mean = 0, sd = x^2+0.5) # heteroskedastic errors
y <- 20 + 3*x + u
reg_ls <- lm(y ~ x)
auxreg <- lm(log(residuals(reg_ls)^2) ~ log(x))
reg_fgls <- lm(log(y) ~ log(x), weights = 1/exp(fitted(auxreg)))
Как вы рекомендуете действовать в такой ситуации, когда независимая переменная содержит нули?