Как найти существенные коэффициенты

Я могу получить значимые коэффициенты при использовании команды summary. Здесь я вижу, что Sepal.Width и Petal.Length имеют намного более низкие значения P, глядя на звезды. Однако эти звезды исчезают, если я только распечатываю коэффициенты. Как мне вернуть звезды? (Примечание: Я хочу, чтобы звезды на моем выходе.)Как найти существенные коэффициенты

library(plm) 
m1 = plm(Sepal.Length ~ Sepal.Width + Petal.Length + Petal.Width, data=iris, index=c('Species')) 

summary(m1) 
Oneway (individual) effect Within Model 

Call: 
plm(formula = Sepal.Length ~ Sepal.Width + Petal.Length + Petal.Width, 
    data = iris, index = c("Species")) 

Balanced Panel: n=3, T=50, N=150 

Residuals : 
    Min. 1st Qu. Median 3rd Qu.  Max. 
-0.79400 -0.21900 0.00899 0.20300 0.73100 

Coefficients : 
       Estimate Std. Error t-value Pr(>|t|)  
Sepal.Width 0.495889 0.086070 5.7615 4.868e-08 *** 
Petal.Length 0.829244 0.068528 12.1009 < 2.2e-16 *** 
Petal.Width -0.315155 0.151196 -2.0844 0.03889 * 
--- 
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Total Sum of Squares: 38.956 
Residual Sum of Squares: 13.556 
R-Squared:  0.65201 
Adj. R-Squared: 0.62593 
F-statistic: 89.9338 on 3 and 144 DF, p-value: < 2.22e-16

Здесь я пытаюсь просто распечатать коэффициенты. Мои знаменательные звезды исчезли.

summary(m1)$coeff 
        Estimate Std. Error t-value  Pr(>|t|) 
    Sepal.Width 0.4958889 0.08606992 5.761466 4.867516e-08 
    Petal.Length 0.8292439 0.06852765 12.100867 1.073592e-23 
    Petal.Width -0.3151552 0.15119575 -2.084418 3.888826e-02

источник

2016-02-15 Bobe Kryant

вдохновение: [Extract только коэффициенты, значения р являются существенными из логистической модели] (http://stackoverflow.com/questions/16070926/extract-only-coefficients-whose-p -значения-значимые-из-логистической модели) –

@Pascal Спасибо за ссылку. Тем не менее, я хочу, чтобы все мои результаты, я просто хочу звезды рядом с ними. Я просматриваю исходный код, я просто не понимаю, почему звезды уходят. –

Звезды отображаются в методе 'print' объекта _summary_. Когда вы используете 'coef', это целый другой зверь. В сводке объясняются звезды (и точки). Их легко воссоздать, но вопрос в том, почему? У вас есть точные значения p. Лучше не бывает. –

Вы можете добавить звезды в итоговые данные. Ниже код делает это:

library(plm) 
m1 = plm(Sepal.Length ~ Sepal.Width + Petal.Length + Petal.Width, data=iris, index=c('Species')) 

summary(m1) 

modelSum <- data.frame(summary(m1)$coeff) 

modelSum$stars <- ' ' 
modelSum$stars <- ifelse(modelSum$Pr...t.. <= 0.1 & modelSum$Pr...t.. > 0.05 ,'.', modelSum$stars) 
modelSum$stars <- ifelse(modelSum$Pr...t.. <= 0.05 & modelSum$Pr...t.. > 0.01 ,'*', modelSum$stars) 
modelSum$stars <- ifelse(modelSum$Pr...t.. <= 0.01 & modelSum$Pr...t.. > 0.001 ,'**', modelSum$stars) 
modelSum$stars <- ifelse(modelSum$Pr...t.. <= 0.001,'***', modelSum$stars) 
modelSum

источник

2016-02-15 06:46:53

Я ожидаю, что в методе печати для 'summary (plmModel)' были добавлены звезды. Пожалуйста, поддержите, если моя душа поможет вам. :) –

Как найти существенные коэффициенты

ответ

Смежные вопросы