如何将nlme对象转换为表格

问题描述 投票：0回答：1

我有以下数据集

df3 = structure(list(A = c("51", "51", "51", "51", "51", "51", "51", 
"51", "51", "51", "51", "51", "51", "51", "51", "51", "51", "51", 
"51", "51", "51", "51", "51", "51", "51", "51", "51", "51", "51", 
"51", "52", "52", "52", "52", "52", "52", "52", "52", "52", "52", 
"52", "52", "52", "52", "52", "52", "52", "52", "52", "52", "52", 
"52", "52", "52", "52", "52", "52", "52", "52", "52", "53", "53", 
"53", "53", "53", "53", "53", "53", "53", "53", "53", "53", "53", 
"53", "53", "53", "53", "53", "53", "53", "53", "53", "53", "53", 
"53", "53", "53", "53", "53", "53", "54", "54", "54", "54", "54", 
"54", "54", "54", "54", "54", "54", "54", "54", "54", "54", "54", 
"54", "54", "54", "54", "54", "54", "54", "54", "54", "54", "54", 
"54", "54", "54", "56", "56", "56", "56", "56", "56", "56", "56", 
"56", "56", "56", "56", "56", "56", "56", "56", "56", "56", "56", 
"56", "56", "56", "56", "56", "56", "56", "56", "56", "56", "56", 
"57", "57", "57", "57", "57", "57", "57", "57", "57", "57", "57", 
"57", "57", "57", "57", "57", "57", "57", "57", "57", "57", "57", 
"57", "57", "57", "57", "57", "57", "57", "57"), B = c(84.5, 
84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 
84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 
84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 72, 72, 72, 72, 72, 
72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 
72, 72, 72, 72, 72, 72, 72, 72, 72, 62, 62, 62, 62, 62, 62, 62, 
62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 
62, 62, 62, 62, 62, 62, 62, 88, 88, 88, 88, 88, 88, 88, 88, 88, 
88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 
88, 88, 88, 88, 88, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 
96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 
96, 96, 96, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 
67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 
67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5), 
    C = c("0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
    "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
    "0", "0", "0", "0", "0", "0", "0", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "0", "0", "0", "0", "0", 
    "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
    "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
    "0"), D = c(73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 
    73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 
    73, 73, 73, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 
    71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 
    71, 71, 71, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 
    72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 
    72, 72, 72, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 
    71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 
    71, 71, 71, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 
    73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 
    73, 73, 73, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 
    74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 
    74, 74, 74), E = c(136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815)), row.names = c(NA, -180L), class = c("tbl_df", 
"tbl", "data.frame"))

我运行了以下代码来获取混合模型

obj = (nlme::lme(E ~ B + C + D, random = ~1|A, df3) %>% 
      summary())$tTable 

flextable(obj)

但是我无法获得显示结果的表格。有谁知道如何使用您的任何软件包（flextable、gtsummary、kableExtra 等）获取该表吗？ 当然我想获取所有列（包括变量名）

                   Value  Std.Error   DF     t-value      p-value
(Intercept) -12.68019966 3.10808245 1827  -4.0797501 4.702378e-05
B            -0.24611480 0.02288007   59 -10.7567325 1.557930e-15
C1            0.36800234 3.57315266   59   0.1029909 9.183194e-01
D             0.08818415 0.02749166   59   3.2076696 2.163542e-03

谢谢

r kable gtsummary flextable nlme
1个回答
2
投票

如果将 

summary()$tTable
转换为数据框（或小标题）

df3 = structure(list(A = c("51", "51", "51", "51", "51", "51", "51", 
"51", "51", "51", "51", "51", "51", "51", "51", "51", "51", "51", 
"51", "51", "51", "51", "51", "51", "51", "51", "51", "51", "51", 
"51", "52", "52", "52", "52", "52", "52", "52", "52", "52", "52", 
"52", "52", "52", "52", "52", "52", "52", "52", "52", "52", "52", 
"52", "52", "52", "52", "52", "52", "52", "52", "52", "53", "53", 
"53", "53", "53", "53", "53", "53", "53", "53", "53", "53", "53", 
"53", "53", "53", "53", "53", "53", "53", "53", "53", "53", "53", 
"53", "53", "53", "53", "53", "53", "54", "54", "54", "54", "54", 
"54", "54", "54", "54", "54", "54", "54", "54", "54", "54", "54", 
"54", "54", "54", "54", "54", "54", "54", "54", "54", "54", "54", 
"54", "54", "54", "56", "56", "56", "56", "56", "56", "56", "56", 
"56", "56", "56", "56", "56", "56", "56", "56", "56", "56", "56", 
"56", "56", "56", "56", "56", "56", "56", "56", "56", "56", "56", 
"57", "57", "57", "57", "57", "57", "57", "57", "57", "57", "57", 
"57", "57", "57", "57", "57", "57", "57", "57", "57", "57", "57", 
"57", "57", "57", "57", "57", "57", "57", "57"), B = c(84.5, 
84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 
84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 
84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 84.5, 72, 72, 72, 72, 72, 
72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 
72, 72, 72, 72, 72, 72, 72, 72, 72, 62, 62, 62, 62, 62, 62, 62, 
62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 62, 
62, 62, 62, 62, 62, 62, 62, 88, 88, 88, 88, 88, 88, 88, 88, 88, 
88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 
88, 88, 88, 88, 88, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 
96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 
96, 96, 96, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 
67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 
67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5, 67.5), 
    C = c("0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
    "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
    "0", "0", "0", "0", "0", "0", "0", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", "1", 
    "1", "1", "1", "1", "1", "1", "1", "0", "0", "0", "0", "0", 
    "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
    "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
    "0"), D = c(73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 
    73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 
    73, 73, 73, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 
    71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 
    71, 71, 71, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 
    72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 
    72, 72, 72, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 
    71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 
    71, 71, 71, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 
    73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 
    73, 73, 73, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 
    74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 
    74, 74, 74), E = c(136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 136.061064886813, 
    136.061064886813, 136.061064886813, 136.061064886813, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 91.138451650252, 91.138451650252, 91.138451650252, 
    91.138451650252, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 144.422871899338, 
    144.422871899338, 144.422871899338, 144.422871899338, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 105.603501603159, 105.603501603159, 105.603501603159, 
    105.603501603159, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 128.444375466774, 
    128.444375466774, 128.444375466774, 128.444375466774, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815, 57.1855792612815, 57.1855792612815, 57.1855792612815, 
    57.1855792612815)), row.names = c(NA, -180L), class = c("tbl_df", 
"tbl", "data.frame"))  

library(dplyr)
library(knitr)
mod <- nlme::lme(E ~ B + C + D, random = ~1|A, df3)
summary(mod)$tTable %>% as_tibble(rownames="variable") %>% knitr::kable(digits=3)
变量 价值 标准误差 DF t 值 p 值
（拦截） 93.900 0.382 174 245.532 0
B 0.697 0.006 2 123.619 0
C1 -75.549 0.633 2 -119.272 0
D -1.499 0.006 2 -268.038 0

创建于 2023-09-04，使用 reprex v2.0.2

最新问题
© www.soinside.com 2019 - 2024. All rights reserved.