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 | |
