1

我正在运行一个未标记的占用模型，并且我的协变量之一是受保护状态。这个协变量有 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

数据是标记为“TN1”“TN2”“TN3”等的 111 个站点，尽管在形成占用框架之前已将其删除，因此与此无关。我想问题是有多个站点具有相同的受保护状态，因此具有相同的“因子值”，但我不知道如何解决这个问题。

4

0 回答 0