我正在构建一个具有按功能分组的多元回归模型,因此接近按每个sub_category分组的100多个模型。我想添加Adj。新数据框中的R sq。和R sq。值具有系数。我能够创建一个只有回归系数但无法添加R sq。值的表。
R代码:
dat <- data.table(main_df)
df_coeff<- dat[,list(intercept=coef(lm(y~x1+x2+x3))[1], coef_price=coef(lmlm(y~x1+x2+x3)))[2],coef_discount=coef(lmlm(y~x1+x2+x3)))[3], coef_clicks=coef(lm(y~x1+x2+x3)))[4], rsq= summary((lm(lm(y~x1+x2+x3)))$r.squared[5]), by=subcatg_id]
我有这个目前的表结构,并希望将r.sq添加到现有表中,以便可以看到每个模型的R.sq。[现有表但没有rsq列] [1]
我目前拥有的此数据帧的R代码:
df_coeff<- dat[,list(intercept=coef(lm(y~x1+x2+x3))[1], coef_price=coef(lmlm(y~x1+x2+x3)))[2],coef_discount=coef(lmlm(y~x1+x2+x3)))[3], coef_clicks=coef(lm(y~x1+x2+x3)))[4], by=subcatg_id]
#Calculating elasticity
df_coeff$exp_price<- exp(df_coeff$coef_price)
df_coeff$exp_discount<- exp(df_coeff$coef_discount)
df_coeff$exp_clicks<- exp(df_coeff$coef_clicks)
df_coeff$exp_elasticity<- df_coeff$exp_price+df_coeff$exp_discount+df_coeff$exp_clicks
df_coeff$elasticity<- apply(cbind(df_coeff$exp_price,df_coeff$exp_discount,df_coeff$exp_clicks),1,sum,na.rm = TRUE)