我有一个高度偏斜的结果变量TnI,我使用log2对其进行了转换。我使用了出色的 rms 包来绘制 OLS 预测。是否可以对 log2(TnI( 进行指数化以获得预测因子与 TnI 而不是 log2(TnI( 的图?非常感谢,安娜玛丽
dd <- datadist(df); options(datadist="dd")
m1 <- ols(log2(tni) ~ age + ischaemia, data=df, x=TRUE, y=TRUE)
plot(Predict(m1)
是的;变换xx = log_2(x)
的逆是x = 2^xx
。更一般地说,使用基b
获取x
日志的反面是b^x
。
例如:
> x <- c(4,8)
> xx <- log2(x)
> 2^xx
[1] 4 8
对于Predict()
函数/方法,您需要执行
predvals <- 2^Predict(m1)$yhat
或类似的东西,从 Predict()
返回的数据框中提取预测值。如果存在lower
和upper
的分量(如果您请求置信区间(,则可以以相同的方式提取和反向转换这些分量。