Связано ли предупреждающее сообщение "усечь до двоичного кода" с сообщением об ошибке "hessian is singular"?
Когда я выполняю приведенный ниже код, я получаю предупреждающее сообщение;
В truncateToBinary(y): некоторые наблюдения были> 1. Они были усечены до 1.
У меня 111 сайтов и 4 ковариаты сайтов. 2 - значения, варьирующиеся от 0 до 1, а другие 2 - стандартные ковариаты. Когда я запустил код для извлечения моделей из моего немаркированного кадра занятости для моих нестандартизированных ковариат, я получил графики, хотя они были ужасными, т.е. доверительный интервал покрывает весь график. Для моделей, использующих стандартизованные ковариаты, я получил сообщение об ошибке "гессиан сингулярен".
Связаны ли сообщение об ошибке и предыдущее предупреждение?
Я проверил свою матрицу обнаружения (это y), и нет наблюдений, превышающих 1. Относится ли это предупреждающее сообщение к чему-то другому?
Мне трудно поверить, что проблема связана с небольшим размером выборки, когда я запускал тот же фрагмент кода с меньшим набором данных, я не получал этого сообщения об ошибке.
length(WolfDetMat2[,1])
WolfDetMat3<-WolfDetMat2[,-1]
Site.Covs.WolvesB <- read.csv("wolftotcovars18.csv", header = T)
#Site.Covs.Wolves1B<-Site.Covs.WolvesB[,2:5];rownames(Site.Covs.Wolves1B)=Site.Covs.WolvesB[,1]
#head(Site.Covs.Wolves1B)
site.SAL <- read.csv("SALobscovs8.csv", header = T)
site.TAN <- read.csv("TANobscovs8.csv", header = T)
EFF <- read.csv("WolfEFF4.csv", header = T, row.names = 1, na.strings = "NA")
# The Model
WolvesModelB <- unmarkedFrameOccu(y=WolfDetMat3, siteCovs=Site.Covs.WolvesB,
obsCovs=list(site.SAL=site.SAL[,c("occ_1", "occ_2", "occ_3", "occ_4", "occ_5", "occ_6", "occ_7", "occ_8", "occ_9", "occ_10", "occ_11", "occ_12", "occ_13", "occ_14", "occ_15", "occ_16", "occ_17", "occ_18", "occ_19", "occ_20")],
site.TAN=site.TAN[,c("occ_1", "occ_2", "occ_3", "occ_4", "occ_5", "occ_6", "occ_7", "occ_8", "occ_9", "occ_10", "occ_11", "occ_12", "occ_13", "occ_14", "occ_15", "occ_16", "occ_17", "occ_18", "occ_19", "occ_20")],
EFF=EFF[,c("occ_1", "occ_2", "occ_3", "occ_4", "occ_5", "occ_6", "occ_7", "occ_8", "occ_9", "occ_10", "occ_11", "occ_12", "occ_13", "occ_14", "occ_15", "occ_16", "occ_17", "occ_18", "occ_19", "occ_20")]))
m<-occu(~1~1, data=WolvesModelB)
m.site2<-occu(~ site.SAL + site.TAN + EFF ~ 1, data=WolvesModelB)
m.VIL2<-occu(~ site.SAL + site.TAN + EFF ~ Vill_Dist, data=WolvesModelB)
m.RAN2<-occu(~ site.SAL + site.TAN + EFF ~ Rang_Dist, data=WolvesModelB)
m.LEP2<-occu(~ site.SAL + site.TAN + EFF ~ Lep_Det, data=WolvesModelB)
m.PRO2<-occu(~ site.SAL + site.TAN + EFF ~ Prot_Stat, data = WolvesModelB)