我想使用 R 中的 Survey 包,使用
svyciprop()
函数来估计各种比例及其置信区间,就像这样
data(api)
dclus1<-svydesign(id=~dnum, fpc=~fpc, data=apiclus1)
svyciprop(~I(ell==0), dclus1, method="li")
我认为我理解R中
I()
函数的原理,但我不明白在这种特殊情况下需要这个函数。事实上,如果我省略 I()
,就会产生相同的结果。
data(api)
dclus1<-svydesign(id=~dnum, fpc=~fpc, data=apiclus1)
svyciprop(~ell==0, dclus1, method="li")
感谢您的帮助并且如果需要更多信息,请告诉我。
公式中的
I()
表示计算表达式 ell==0
并将结果用作预测变量。 当表达式包含 +
、-
等在模型公式中专门解释的运算符时,这是必要的。
在这种情况下没有必要,因为
==
在模型公式中并不意味着特殊的东西。 我假设它是为了一致性而使用的:每当您想在模型公式中构造变量时,请将构造包装在 I()
中,以避免使用被误解的运算符。