出力データのサブセットで線形回帰を使用してR-二乗

Question

lag,psit、およびvarのデータセットが変数としてあります。私はPSITがVAR lm(psit~var, df)によって予測された線形回帰を実行したいlagの各サブセットに対して

set.seed(560) 
df<-data.frame(lag= rep(1:40, each=228), psit= rep(rnorm(228, 20, 
10)),var=rnorm(9120, 50, 10)) 

。私は、R-二乗値、F-統計量、およびモデルのp-値をデータフレームに出力したいと思います。

output<-data.frame(lag= rep(1:40, each=1), rsqrd= rep(rnorm(40, .5, 
0.01)),fstat=rnorm(40,5, 0.05), pvalue=rnorm(40,0.01, 0.1)) 

私が試してみました：

library(broom) 
rsqrd <- 
    models %>% 
    ungroup() %>% group_by(lag) %>% 
    do(glance(.$model[[1]])) 

glanceがかかる：

models <- 
    df %>% 
    group_by(lag) %>% 
    do(model = lm(psit ~ var,data = .)) 


rsqrd<- 
    model %>% 
    do(data.frame(
    lag = .$lag, 
    summary(.)$r.squared, summary(.)$fstatistic,summary(.)$pvalue)) 

Answer 1

あなたは線形モデルのシリーズから必要な要約統計量を取得するためにほうきライブラリのglance機能を使用することができます単一線形モデルを作成し、要約統計量の1行のデータフレームを返します。