r - lm 输出的逐列 cbind

Question

我想对两个 lm 输出进行逐列比较。类似这样的事情，除了我希望看到估计值旁边的估计值嵌套到估计值，Std Error 旁边的 Std Error 等。

    x=runif(6);y=runif(6);z=runif(6)
    x2=runif(6);y2=runif(6);z2=runif(6)
    lm1=lm(x~z+y)
    lm2=lm(x2~z2+y2)
    s=summary(lm1)[["coefficients"]]
    s2=summary(lm2)[["coefficients"]]
    > cbind(s,s2)
               Estimate Std. Error    t value  Pr(>|t|)   Estimate Std. Error   t value   Pr(>|t|)
(Intercept)  0.67359508  0.8137390  0.8277778 0.4684978  0.6944551  0.2398032  2.895938 0.06270911
z           -0.19911923  1.0491558 -0.1897899 0.8615892  1.2509814  0.5678805  2.202896 0.11486078
y           -0.07097291  0.4569228 -0.1553280 0.8864254 -0.8926959  0.4984554 -1.790924 0.17123009

Answer 1

试试这个:

# column bind
res <- cbind(s,s2)

# make unique column names
colnames(res) <- make.unique(colnames(res))

# sort column names
ix <- order(colnames(res))

# result
res[,ix]

# output

#               Estimate Estimate.1  Pr(>|t|) Pr(>|t|).1 Std. Error Std. Error.1    t value t value.1
# (Intercept)  1.7701957 0.09783934 0.5505218  0.7042956  2.6399847    0.2342630  0.6705325 0.4176475
# z            0.2393442 0.45691026 0.8111899  0.6096530  0.9180588    0.8041508  0.2607068 0.5681898
# y           -1.7020430 0.12906057 0.6185858  0.8798162  3.0758497    0.7846930 -0.5533570 0.1644727