R/Python 中的 Stata 图
Stata plot in R/Python
R 中的以下绘图使用了哪个函数？
至少对我来说它看起来像一个预定义的函数。
编辑：好的，根据 Claudio，它似乎是 Stata。
新问题：
python/R 中是否有任何类似的东西可以得到这个输出？
系数怎么计算？这是什么系数？
您可以使用 pystata 包，它允许您在 Python 环境中 运行 stata。您可以从 here
下载
我们可以在 R 中使用 lm() 获得相同的信息，并结合 broom 和 kable 获得表格：
model <- lm(mpg ~ disp + qsec, data=mtcars)
summary(model) # show results
coefficients(model) # model coefficients
confint(model, level=0.95) # CIs for model parameters
fitted(model) # predicted values
residuals(model) # residuals
anova(model) # anova table
vcov(model) # covariance matrix for model parameters
influence(model) # regression diagnostics
Combining with:
library(broom) # for tidy()
library(knitr) # for kable()
out <- tidy(model)
out1 <- tidy(anova(model))
kable(out1)
kable(out)
输出：
|term | df| sumsq| meansq| statistic| p.value|
|:---------|--:|---------:|---------:|----------:|---------:|
|disp | 1| 808.88850| 808.88850| 74.8166388| 0.0000000|
|qsec | 1| 3.62193| 3.62193| 0.3350037| 0.5671961|
|Residuals | 29| 313.53676| 10.81161| NA| NA|
> kable(out)
|term | estimate| std.error| statistic| p.value|
|:-----------|----------:|---------:|----------:|---------:|
|(Intercept) | 25.5045079| 7.1840940| 3.5501356| 0.0013359|
|disp | -0.0398877| 0.0052882| -7.5428272| 0.0000000|
|qsec | 0.2122880| 0.3667758| 0.5787951| 0.5671961|
R 中的以下绘图使用了哪个函数？ 至少对我来说它看起来像一个预定义的函数。
编辑：好的，根据 Claudio，它似乎是 Stata。
新问题： python/R 中是否有任何类似的东西可以得到这个输出？ 系数怎么计算？这是什么系数？
您可以使用 pystata 包，它允许您在 Python 环境中 运行 stata。您可以从 here
我们可以在 R 中使用 lm() 获得相同的信息，并结合 broom 和 kable 获得表格：
model <- lm(mpg ~ disp + qsec, data=mtcars)
summary(model) # show results
coefficients(model) # model coefficients
confint(model, level=0.95) # CIs for model parameters
fitted(model) # predicted values
residuals(model) # residuals
anova(model) # anova table
vcov(model) # covariance matrix for model parameters
influence(model) # regression diagnostics
Combining with:
library(broom) # for tidy()
library(knitr) # for kable()
out <- tidy(model)
out1 <- tidy(anova(model))
kable(out1)
kable(out)
输出：
|term | df| sumsq| meansq| statistic| p.value|
|:---------|--:|---------:|---------:|----------:|---------:|
|disp | 1| 808.88850| 808.88850| 74.8166388| 0.0000000|
|qsec | 1| 3.62193| 3.62193| 0.3350037| 0.5671961|
|Residuals | 29| 313.53676| 10.81161| NA| NA|
> kable(out)
|term | estimate| std.error| statistic| p.value|
|:-----------|----------:|---------:|----------:|---------:|
|(Intercept) | 25.5045079| 7.1840940| 3.5501356| 0.0013359|
|disp | -0.0398877| 0.0052882| -7.5428272| 0.0000000|
|qsec | 0.2122880| 0.3667758| 0.5787951| 0.5671961|