我有一系列线性模型,我想为每个模型报告标准化系数。然而,当我在stargazer中打印模型时,它似乎会自动将标准化系数的显著性星号打印出来,就好像它们是未标准化系数一样。你可以看到下面的差异。
在未标准化值的基础上打印显著性星号是否在统计上是错误的?在stargazer中如何实现这一点?谢谢!
在未标准化值的基础上打印显著性星号是否在统计上是错误的?在stargazer中如何实现这一点?谢谢!
#load libraries
library(stargazer)
library(lm.beta)
#fake data
var1<-rnorm(100, mean=10, sd=5)
var2<-rnorm(100, mean=5, sd=2)
var3<-rnorm(100, mean=2, sd=3)
var4<-rnorm(100, mean=5, sd=1)
df<-data.frame(var1, var2, var3, var4)
#model with unstandardized betas
model1<-lm(var1~var2+var3+var4, data=df)
#Standardized betas
model1.beta<-lm.beta(model1)
#print
stargazer(model1, model1.beta, type='text')