r - R上的双变量概率问题

Question

我尝试使用 GJRM 包在 R 上做一个 bivariateprobit，但没有成功。我有两个方程，像往常一样有两个不同的因变量，但每个方程都有不同的解释变量集。

这是我已经编码的内容：

EQ1 <- emploi10 ~ SM + femme + jeune + jeune_ad + ad_trent + ad_quarant + ad_vieux + dip_manq + dip_sans + dip_sup_court + dip_sup_long + CSP_manquant + salarie_prive + salarie_privé_cadre_dg + salarie_public_cadre + salarie_public + exp_mini + exp_moy  + chom_un + chom_deux + indep +public + enf_06 + abs_couple_2006_2010 + couple_2006 + couple_2010 + percep_bas + percep_moy + mal_chroniques + limitations + pb_faible + pb_élevé + pb_tr_élevé
EQ2 <- SM ~ élevé_seul + violences + femme + jeune + jeune_ad + ad_trent + ad_quarant + ad_vieux + dip_manq + dip_sans + dip_sup_court + dip_sup_long + CSP_manquant + salarie_prive + salarie_privé_cadre_dg + salarie_public_cadre + salarie_public + exp_mini + exp_moy  + chom_un + chom_deux + indep +public + enf_06 + abs_couple_2006_2010 + couple_2006 + couple_2010 + percep_bas + percep_moy + mal_chroniques + limitations + pb_faible + pb_élevé + pb_tr_élevé
f.list <- list(EQ1, EQ2)
mr <- c("probit", "probit")
library(GJRM)
bvp <- gjrm(f.list,data = data2, Model="B", margins=mr)

但是，在这里我收到了以下 R 消息：变量“data2”的类型（列表）无效。我不明白有什么问题。我尝试编写data=as.data.frame(data2)但收到相同的错误消息。

任何人都可以提供帮助，因为我完全按照我的老师所做的吗？

非常感谢

阿切尔

标签： rlogistic-regression