монотонно увеличивающаяся ошибка в точках оценки fhat при отображении результатов kde.local.test из пакета ks - PullRequest
Новая
       9

монотонно увеличивающаяся ошибка в точках оценки fhat при отображении результатов kde.local.test из пакета ks

0 голосов
/ 28 сентября 2018

У меня проблемы с отображением результатов теста функции kde.local.test (пакет ks из R) для локальных различий в плотности ядра между двумя образцами с двухмерной оценкой плотности ядра

Я подготовил дваПримеры.Пример 1 работает, пример 2 - нет.Почему?

Пример № 1: не работает 

data("iris")  
dat <- iris[order(iris$Sepal.Width, iris$Sepal.Length),]  
x1 <- iris[1:75,c("Sepal.Length", "Sepal.Width")]  
x2 <- iris[76:150,c("Sepal.Length", "Sepal.Width")]  
H1samse <- Hpi(x1, pre = "scale", nstage = 2, pilot = "samse") 
H2samse <- Hpi(x2, pre = "scale", nstage = 2, pilot = "samse")  
loct <- kde.local.test(x1,x2, H1 = H1samse, H2 = H2samse)  
plot(loct)

В результате выдается следующая ошибка: 

Error in .filled.contour(fhat$eval.points[[1]], fhat$eval.points[[2]],  : 
  invalid contour values: have to increase strictly monotonously

Пример# 2: отлично работает 

data("crabs")  
x1 <- crabs[crabs$sex=="M", c("FL", "RW")]  
x2 <- crabs[crabs$sex=="F", c("FL", "RW")]  
H1samse <- Hpi(x1, pre = "sphere", nstage = 2, pilot = "samse")  
H2samse <- Hpi(x2, pre = "sphere", nstage = 2, pilot = "samse")  
loct <- kde.local.test(x1,x2, H1 = H1samse, H2 = H2samse)  
plot(loct)  
H1scv<- Hscv(x1)  
H2scv <- Hscv(x2)  
loct <- kde.local.test(x1,x2, H1 = H1scv, H2 = H2scv)  
plot(loct)
...