S********a 发帖数: 359 | 1 > unadj.lrm
结果如下
Coef S.E. Wald Z P
Intercept -3.60925 0.132424 -27.26 0.0000
pmavg 0.02205 0.009483 2.33 0.0201
我想取pmavg 的standard error (i.e.,0.009483)
怎么样写R code呢?
> unadj.lrm$SE 好像不对呀 |
m**********r 发帖数: 507 | |
hs 发帖数: 1549 | 3 summary(xxxx)$coef[2,2]
xxxx是 unadj.lrm
【在 S********a 的大作中提到】 : > unadj.lrm : 结果如下 : Coef S.E. Wald Z P : Intercept -3.60925 0.132424 -27.26 0.0000 : pmavg 0.02205 0.009483 2.33 0.0201 : 我想取pmavg 的standard error (i.e.,0.009483) : 怎么样写R code呢? : > unadj.lrm$SE 好像不对呀
|
S********a 发帖数: 359 | 4 > unadj.lrm[2,2]
Error in x[..., drop = drop] : incorrect number of dimensions
咋办?
【在 m**********r 的大作中提到】 : unadj.lrm[2,2] 不行么?
|
S********a 发帖数: 359 | 5 > summary(unadj.lrm)
Error in summary.Design(unadj.lrm) :
> summary(unadj.lrm)$coeff
Error in summary.Design(unadj.lrm) :
adjustment values not defined here or with datadist for pmavg
怎么回事呀?
adjustment values not defined here or with datadist for pmavg
【在 hs 的大作中提到】 : summary(xxxx)$coef[2,2] : xxxx是 unadj.lrm
|
S********a 发帖数: 359 | 6 完整的结果是这样的, 如何取SE呢?谢谢大家
> unadj.lrm
Logistic Regression Model
lrm(formula = combpreecl ~ pmavg, data = mydata2)
Frequencies of Responses
0 1
118691 4360
Obs Max Deriv Model L.R. d.f. P C 123051
7e-10 5.41 1 0.02 0.504
Dxy Gamma Tau-a R2 Brier
0.008 0.015 0.001 0 0.034
Coef S.E. Wald Z P
Intercept -3.60925 0.132424 -27.26 0.0000
pmavg 0.02205 0.009483 2.33 0.0201 |
N****n 发帖数: 1208 | 7 Try this?
http://r.789695.n4.nabble.com/extracting-the-standard-error-in-
123051
【在 S********a 的大作中提到】 : 完整的结果是这样的, 如何取SE呢?谢谢大家 : > unadj.lrm : Logistic Regression Model : lrm(formula = combpreecl ~ pmavg, data = mydata2) : Frequencies of Responses : 0 1 : 118691 4360 : Obs Max Deriv Model L.R. d.f. P C 123051 : 7e-10 5.41 1 0.02 0.504 : Dxy Gamma Tau-a R2 Brier
|
S********a 发帖数: 359 | 8 谢谢,我用在lrm结果上不行
> sqrt(diag(vcov(adj.lrm)))
Error in if (object$family$family %in% c("poisson", "binomial")) 1 else if (
df.r > :
但是用在
> sqrt(diag(vcov(adj.glm))) 就成功了
难道lrm不如glm好用吗?
【在 N****n 的大作中提到】 : Try this? : http://r.789695.n4.nabble.com/extracting-the-standard-error-in- : : 123051
|
N****n 发帖数: 1208 | 9 Results from some of models can be extracted, some are not.
I think that's drawback for R .... still under development, I guess?
(
【在 S********a 的大作中提到】 : 谢谢,我用在lrm结果上不行 : > sqrt(diag(vcov(adj.lrm))) : Error in if (object$family$family %in% c("poisson", "binomial")) 1 else if ( : df.r > : : 但是用在 : > sqrt(diag(vcov(adj.glm))) 就成功了 : 难道lrm不如glm好用吗?
|
i********f 发帖数: 206 | 10 一个比较笨的办法,但是能用
sqrt(diag(yourResult$var))[2]
【在 S********a 的大作中提到】 : 谢谢,我用在lrm结果上不行 : > sqrt(diag(vcov(adj.lrm))) : Error in if (object$family$family %in% c("poisson", "binomial")) 1 else if ( : df.r > : : 但是用在 : > sqrt(diag(vcov(adj.glm))) 就成功了 : 难道lrm不如glm好用吗?
|
B******5 发帖数: 4676 | 11 names(unadj.lrm)或者slot(unadj.lrm)有什么? |
S********a 发帖数: 359 | 12 > sqrt(diag(adj.lrm$var))[12]
pmavg
0.009879157
> sqrt(diag(adj.glm$var))[12]
[1] NA
> summary(adj.glm)$coef[12,2]
[1] 0.009879123
对lrm用你的这个能得出结果,对glm用summary能得出结果。
谢谢!
【在 i********f 的大作中提到】 : 一个比较笨的办法,但是能用 : sqrt(diag(yourResult$var))[2]
|
i********f 发帖数: 206 | |