パッケージにnetlm()
というネットワーク回帰を実行しました。すべての結果は正常に見えます。しかし、私は各係数の標準誤差を見つけることができません。返されるモデルオブジェクトには標準エラーはありません。 summary()
もエクスポートしません。それらを得る方法はありますか?`sna`パッケージで` netlm() `を使用した場合の回帰係数の標準誤差の求め方は?
IV0 <- list(SEX, GRADE, YEAR)
M0 <- netlm(K1M, IV0, nullhyp = "qapspp", reps = 2000)
summary.default(M0) ## check structure of returned model object
# Length Class Mode
# coefficients 4 -none- numeric
# fitted.values 812 -none- numeric
# residuals 812 -none- numeric
# qr 4 qr list
# rank 1 -none- numeric
# n 1 -none- numeric
# df.residual 1 -none- numeric
# tstat 4 -none- numeric
# dist 8000 -none- numeric
# pleeq 4 -none- numeric
# pgreq 4 -none- numeric
# pgreqabs 4 -none- numeric
# nullhyp 1 -none- character
# names 4 -none- character
# intercept 1 -none- logical
M0 ## print model object
# OLS Network Model
# Residuals:
# 0% 25% 50% 75% 100%
# -0.3669251 -0.3376203 -0.3066127 0.6623797 0.7340360
# Coefficients:
# Estimate Pr(<=b) Pr(>=b) Pr(>=|b|)
# (intercept) 0.271072656 0.966 0.034 0.0605
# x1 0.009641084 0.602 0.398 0.8360
# x2 -0.031007609 0.169 0.831 0.3440
# x3 0.029304737 0.774 0.226 0.4660
# Residual standard error: 0.47 on 808 degrees of freedom
# Multiple R-squared: 0.002094 Adjusted R-squared: -0.001611
# F-statistic: 0.5651 on 3 and 808 degrees of freedom, p-value: 0.6382
summary(M0) ## print model summary
# Test Diagnostics:
# Null Hypothesis: qapspp
# Replications: 2000
# Coefficient Distribution Summary:
# (intercept) x1 x2 x3
#Min -3.15189 -4.84538 -3.08131 -2.75036
#1stQ -0.65506 -0.99759 -0.65075 -0.68183
#Median 0.01805 -0.09364 0.06947 -0.01831
#Mean 0.02510 -0.03473 0.03247 -0.01179
#3rdQ 0.69636 0.90936 0.72701 0.63639
#Max 3.29170 5.71549 2.71428 2.90498