在R中,如果我写:
reg1 <- lm(y ~ x, data = ds)回归模型信息存储在具有列表结构的reg1中。
我可以写:
value.of.intercept <- reg1$coefficients[1] 将截取的值保存到名为"value.of.intercept“的变量中
我该怎么写才能将截取的统计显著性值保存到变量"p.value.of.intercept“中?
我可以在总结中看到它的价值。如果我写道:
summary(reg1)截取的统计显著性值紧接在此文本下:Pr(>|t|)
发布于 2015-07-29 06:27:52
您可以设置summary对象的子集。您可以使用
coef(summary(reg1))为了获得带注释的系数,您可以使用以下命令获取截取的p值
coef(summary(reg1))[1, 4]发布于 2015-07-29 18:22:06
如果你使用扫把包中的tidy函数,你会得到一个数据帧中的所有系数。
library(broom)
lm_model <- lm(mpg ~ ., data = mtcars)
lm_model_coefficients <- tidy(lm_model)
lm_model_coefficients
term estimate std.error statistic p.value
1 (Intercept) 12.30337416 18.71788443 0.6573058 0.51812440
2 cyl -0.11144048 1.04502336 -0.1066392 0.91608738
3 disp 0.01333524 0.01785750 0.7467585 0.46348865
4 hp -0.02148212 0.02176858 -0.9868407 0.33495531
5 drat 0.78711097 1.63537307 0.4813036 0.63527790
6 wt -3.71530393 1.89441430 -1.9611887 0.06325215
7 qsec 0.82104075 0.73084480 1.1234133 0.27394127
8 vs 0.31776281 2.10450861 0.1509915 0.88142347
9 am 2.52022689 2.05665055 1.2254035 0.23398971
10 gear 0.65541302 1.49325996 0.4389142 0.66520643
11 carb -0.19941925 0.82875250 -0.2406258 0.81217871https://stackoverflow.com/questions/31688168
复制相似问题