Оценка лямбды для преобразования Йо и Джонсона



У меня есть временной ряд значений осадков в файле csv.Я построил гистограмму данных. Гистограмма смещена влево. Я хотел преобразовать значения так, чтобы они имели нормальное распределение. Я использовал преобразование Йео-Джонсона, доступное в R. преобразованные значения здесь.



Мой вопрос:



В приведенном выше преобразовании я использовал тестовое значение 0,5 для лямбды, которое отлично работает. Есть ли возможность определить оптимальное значение лямбды на основе временной ряд? Буду признателен за любые предложения.



Пока что вот код:



library(car)
dat <- scan("Zamboanga.csv")
hist(dat)
trans <- yjPower(dat,0.5,jacobian.adjusted=TRUE)
hist(trans)


Вот csv-файл .

678   1  

1 ответ:

Сначала найдите оптимальную лямбду, используя функцию boxCox из пакета car для оценки λ по максимальному правдоподобию.

Вы можете построить его следующим образом:

boxCox(your_model, family="yjPower", plotit = TRUE)

пример из резюме

Как сказал Бен Болкер в комментарии, модель здесь может быть чем-то вроде
your_model <- lm(dat~1)

Затем используйте оптимизированную лямбду в существующем коде.

Comments

    Ничего не найдено.