R语言的regmedint包c_cond语句怎么使用？

dxy_bqvnton8

2023-01-30

最近在用R语言regmedint包做中介分析，有些地方不太懂，就是CVAR,C_COND语句，R语言官网给的代码只有一个协变量，C_COND是3，不清楚想加入多个协变量应该怎么加？c_cond是根据什么设置的？下面是自己尝试加协变量跑出来的，出现报错，但不知道怎么改，望各位大佬指点，做毕业论文比较急！

regmedint_obj1 <- regmedint(data = be,

## Variables

yvar = "sleeptime.group7",

avar = "AZI",

mvar = "log5HT",

cvar = c("age.1","SEX"),

## Values at which effects are evaluated

a0 = 0,

a1 = 1,

m_cde = 1,

c_cond = 4,

## Model types

mreg = "linear",

yreg = "logistic",

## Additional specification

interaction = TRUE,

casecontrol = FALSE)

Error: length(c_cond) not equal to length(cvar)

2 个回答

土井挞克树

2023-02-01

有帮助

你的思路基本是对的协变量部位用括号括起来就可以了

loveliufudan

2023-01-30

有帮助

在R语言的regmedint包中，CVAR表示协变量，C_COND表示协变量的条件。如果你要加入多个协变量，可以将协变量写为一个向量，如：cvar = c("age.1","SEX")。C_COND的长度必须与协变量的长度相同，因此，如果你有两个协变量，C_COND的长度应该为2，如：c_cond = c(1, 2)。该值的设置取决于协变量的类型，详情请参见regmedint的官方文档。

报错的原因是：变量cvar的数量为2，而c_cond的数量为1。

c_cond参数是表示cvar（协变量）的状态，其值应与cvar的数量相等。

如果要添加多个协变量，可以将c_cond改为一个长度为cvar长度的数组，每个值表示相应的cvar的状态。例如：

c_cond = c(3, 2)