Question

我试图在R中计算变量的比值比，不仅是线性的，还有逻辑回归中的二次项。假设模型中有X和X ^ 2。当X取特定值时，我知道如何获得优势比（X的单位变化），但我不知道如何计算此估计的置信区间。我在SAS中找到了这个参考：http://support.sas.com/kb/35/189.html，但是我想在R中做。有什么建议吗？

@BenBolker 这是一个例子：

mydata <-read.csv("http://www.ats.ucla.edu/stat/data/binary.csv")
mydata <- transform(mydata, gpaSquared=gpa^2,greSquared=gre^2)
model <- glm(admit ~ gpa + gpaSquared  + gre , family = binomial(logit), data = mydata)

在此示例中，gpa的优势比取决于gpa的实际值（例如，如果gpa = 4，则单位变化的影响以gpa为单位）。我可以计算gpa = 5和gpa = 4的对数几率，并从中获得优势比，但我不知道如何获得OR的CI。（请忽略在示例中，平方术语不是重要的......）

Answer 1

m <- glm(x~X1^2+X2,data,family=binomial(link="logit"))
summary(m)
confint(m) # 95% CI for the coefficients using profiled log-likelihood
confint.default(m) ## CIs using standard errors
exp(coef(m)) # exponentiated coefficients
exp(confint(m)) # 95% CI for exponentiated coefficients

如何用包含X的平方的逻辑回归估计X与CI的比值比？

1 个答案: