Я использую модель без опознавательных знаков, и один из моих ковариат защищен. Этот ковариат имеет 3 уровня 0, 1 и 2. Мне нужно взять усредненную по модели оценку эффекта ковариации и занятости. В настоящее время R читает ковариату как непрерывную, а не категориальную переменную, поэтому мне нужно скорректировать свой код. К сожалению, все методы, которые я могу найти на R для того, чтобы сделать данные категориальными, а не непрерывными, предназначены для «моделей S3», и это «модель S4». Я также пытался сказать R, что данные являются фактором, подобным показанному в приведенном ниже коде.
Site.Covs.Wolves2 <- read.csv("wolftotcovars40.csv", header = T, colClasses = c("factor", "factor", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric"))
Но когда я попытался сделать усредненную по модели оценку, я получил сообщение об ошибке
Some models include more than one instance of the parameter of interest.
This may be due to the presence of interaction/polynomial terms, or variables
with similar names:
see "?modavg" for details on variable specification and "exclude" argument
In addition: Warning message:
In modavg.AICunmarkedFitOccu(M.listTOT3, "Prot_Stat", parm.type = "psi") :
Model names have been supplied automatically in the table
Данные представляют собой 111 сайтов, обозначенных «TN1», «TN2», «TN3» и т. Д. c. хотя это удаляется до формирования кадра занятости, так что это не связано с этим. Я предполагаю, что проблема в том, что есть несколько сайтов с одинаковым защищенным статусом и поэтому имеют одинаковое «значение фактора», но я не уверен, как решить эту проблему.