微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

有没有办法可以从此输出摘要中提取值?

如何解决有没有办法可以从此输出摘要中提取值?

如何提取 ('Estimate','std.Error','Pr (> | z |)' 的值>) 从输出结果中提取并将它们放在下面数据示例中的表格格式中?

require(GJRM)
set.seed(123)
x1 <- sample(1:100,size = 20)
bid1 <- sample(c(5,10,20,30),size = 20,replace = T)
bid2 <- sample(c(5,replace = T)
ans1 <- sample(c(1,0),replace = T)
ans2 <- sample(c(1,replace = T)
df <- cbind(x1,bid1,bid2,ans1,ans2)
df <- as.data.frame(df)
treat.eq <- ans1 ~ bid1 + x1
out.eq <- ans2 ~ bid2 + x1
f.list <- list(treat.eq,out.eq)
mr <- c("probit","probit")
## Model
bvp <- gjrm(f.list,data=df,Model="B",margins= mr)
summary(bvp)

解决方法

数据存在于 model$tableP1model$tableP2 中:

library(GJRM)
bvp <- gjrm(f.list,data=df,Model="B",margins= mr)
model <- summary(bvp)
model$tableP1

#               Estimate Std. Error   z value  Pr(>|z|)
#(Intercept)  0.83797602 0.79822170  1.049804 0.2938084
#bid1        -0.04682910 0.03137645 -1.492492 0.1355702
#x1          -0.01065357 0.00993613 -1.072205 0.2836278

model$tableP2

#                Estimate Std. Error    z value  Pr(>|z|)
#(Intercept)  0.434060433 0.64213719  0.6759621 0.4990647
#bid2        -0.012464492 0.02930301 -0.4253656 0.6705702
#x1          -0.005763137 0.01025844 -0.5617947 0.5742559

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。