从lm函数获取分组数据的p值(Get the p-values from the lm function for grouped data)

 我正在尝试使用lm()函数和plyr包为我的数据中的每个段拟合一个模型，因为我的数据是按键分组的。 
 
 我已经设法运行模型并获得系数以及R ^ 2和adj平方但我正在努力与p值。 
 
library("plyr")
#Sample data
test_data <- data.frame(key = c("a","a","a","a","a","b","b","b","b","b"),
y = c(100,180,120,60,140,200,220,240,260,280),
x1 = c(50,60,79,85,90,133,140,120,160,170),
x2 =  c(20,18,47,16,15,25,30,25,20,15))

#model
model_1 <- dlply(test_data, .(key), 
    function(test_data) lm(y ~ x1 +     x2,data = test_data))

#coefficients
ldply(model_1, coef)

#adj r-squared
ldply(model_1, function(x) summary(x)$r.squared)
 
 我试过这个让我获得密钥和p值但它没有变量的名称，我需要能够将输出与模型中的系数合并。 
 
#p-values but missing the variable names
ldply(model_1, function(x) summary(x)$coefficients)[,c(1,5)]
 
 我试图使用Do来适应模型，然后从dplyr包中tidy ，这可以很好地处理一个小数据集，因为它实际上返回了我需要的所有东西，但我的实际数据包含超过1,000个不同的段，并且RStudio最终崩溃。 

I am trying to fit a model for each segment in my data using the lm() function in conjunction with the plyr package because my data is grouped by a key.
 
I've managed to run the model and get the coefficients along with the R^2 & adj r-squared but I am struggling with the p-values.
 
library("plyr")
#Sample data
test_data <- data.frame(key = c("a","a","a","a","a","b","b","b","b","b"),
y = c(100,180,120,60,140,200,220,240,260,280),
x1 = c(50,60,79,85,90,133,140,120,160,170),
x2 =  c(20,18,47,16,15,25,30,25,20,15))

#model
model_1 <- dlply(test_data, .(key), 
    function(test_data) lm(y ~ x1 +     x2,data = test_data))

#coefficients
ldply(model_1, coef)

#adj r-squared
ldply(model_1, function(x) summary(x)$r.squared)
 
I've tried this which gets me the key and the p-value but it doesn't have the names of the variables which I need to be able to merge the output with the coefficients from the model later.
 
#p-values but missing the variable names
ldply(model_1, function(x) summary(x)$coefficients)[,c(1,5)]
 
I've tried to fit the models using Do and then tidy from the dplyr package and this works fine with a small data set because it actually returns everything I need but my actual data contains over 1,000 different segments and RStudio end up crashing.

原文：https://stackoverflow.com/questions/36065447