如何解决在 r (frm) 中使用分数响应模型中的因子
我一直在尝试使用 frm
包将分数响应模型拟合到我的数据中,但无法使其与我的因子变量一起使用。我能够使用因子并在 STATA 软件中拟合模型,所以我知道在这种类型的模型中是可能的。但是,为了进一步分析,我需要模型在 R 中,因此我无法中继 STATA 模型。
到目前为止,我的代码如下所示:
library(frm)
data <- read.csv("data/model_data_shanon_only.csv",fileEncoding="UTF-8-BOM")
data[,'crop19_group']<- as.factor(data[,'crop19_group'])
data[,'field_id']<-as.factor(data[,'field_id'])
data[,'crop_group']<- as.factor(data[,'crop_group'])
data[,'river_type']<- as.factor(data[,'river_type'])
cluster<- data$field_id
y<-data$kde
x<- (data[,c("fild_sz","ptch_sz","ds_t_pt","rivr_ln","ds_t_rv","SH_Index_750"
"crop_group","crop19_group","river_type")])
model<- frm(y,x,linkfrac = "logit",var.cluster = cluster,dfc = TRUE,var.type = "cluster")
我收到此错误:
Error in `[[<-.data.frame`(`*tmp*`,i,value = c(13001L,13001L,: replacement has 100350 rows,data has 11150
当我使用相同的代码但省略因子变量并仅使用数字变量时,模型运行良好。 数据不包含 NA 变量,所以我知道这不是问题。此外,所有因子数据都被转换为数字而不是名称。
有人对如何合并这些因素有任何建议吗?
感谢您的时间。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。