это мой первый вопрос на этом форуме, который я считаю очень полезным! У меня проблема, которую я не знаю, как ее интерпретировать или решить:
У меня есть измерения биохимического параметра (называемого «дромс»), который не достигает нормы. Итак, я сделал преобразование BoxCox, следуя этому коду:
sum(is.na(dataframe$droms))
[1] 0
boxcoxnc(dataframe$droms, method='sw', lambda =seq(-5,5,0.01),plot=TRUE, alpha=0.05, verbose=TRUE)
Error in boxcoxnc(data_turbidez$droms, method = "sw", lambda = seq(-5, :
Data must include positive values. Specify shifting parameter, lambda2
Итак, я изменил код на новый:
boxcoxnc(dataframe$droms, method='sw', lambda2=2, lambda =seq(-5,5,0.01),plot=TRUE, alpha=0.05, verbose=TRUE)
lambda.hat : -0.53
statistic : 0.9895904
p.value : 0.1819954
Result : Transformed data are normal.
И я приступил к применению преобразования:
dataframe$droms_normalized <- (dataframe$droms^-0.53-1)/-0.53
Но когда я пытаюсь снова проверить нормальность, глядя на гистограмму, я не получаю ту же гистограмму, что и преобразование boxcox, и получаю вывод из теста Шапиро Вилка:
W = NaN, p-value = NA
Не могли бы вы помочь мне определить, что я делаю неправильно в преобразовании?
Спасибо большое!