这个问题是关于如何从S4 lme模型摘要对象中删除字符串的一部分.

假设这是我的模型.

# ----preliminaries-------------------------------------------------------
library("lme4")
library("lattice")
library("minqa")
library("nlme")

##-------------------------------------------------------------------------
data("sleepstudy")
##-------------------------------------------------------------------------
sleepstudy$Days[sleepstudy$Days == "0"] <- "Zero"
sleepstudy$Days[sleepstudy$Days == "1"] <- "One"
sleepstudy$Days[sleepstudy$Days == "2"] <- "Two"
sleepstudy$Days[sleepstudy$Days == "3"] <- "Three"
sleepstudy$Days[sleepstudy$Days == "4"] <- "Four"
sleepstudy$Days[sleepstudy$Days == "5"] <- "Five"
sleepstudy$Days[sleepstudy$Days == "6"] <- "Six"
sleepstudy$Days[sleepstudy$Days == "7"] <- "Seven"

## Model
##-------------------------------------------------------------------------
fm1 <- lme(fixed = Reaction ~ Days ,
           random = ~1 | Days, 
           data = sleepstudy)
summary(fm1)

Linear mixed-effects model fit by REML
  Data: sleepstudy 
       AIC      BIC    logLik
  1855.854 1893.483 -915.9268

Random effects:
 Formula: ~1 | Days
        (Intercept) Residual
StdDev:    30.55409 48.61125

Fixed effects:  Reaction ~ Days 
               Value Std.Error  DF   t-value p-value
(Intercept) 336.6295  32.63178 170 10.316003       0
Days9        14.2217  46.14830   0  0.308174     NaN
DaysFive    -28.1110  46.14830   0 -0.609146     NaN
DaysFour    -47.9801  46.14830   0 -1.039693     NaN
DaysOne     -72.1337  46.14830   0 -1.563086     NaN
DaysSeven   -17.8789  46.14830   0 -0.387423     NaN
DaysSix     -24.4512  46.14830   0 -0.529841     NaN
DaysThree   -53.6375  46.14830   0 -1.162285     NaN
DaysTwo     -71.2676  46.14830   0 -1.544317     NaN
DaysZero    -79.9777  46.14830   0 -1.733058     NaN

我喜欢把fixed effects model summary开头的"101"go 掉

Expecting个我的固定效果总结术语看起来像这样.

Fixed effects:  Reaction ~ Days 
               Value Std.Error  DF   t-value p-value
(Intercept) 336.6295  32.63178 170 10.316003       0
9        14.2217  46.14830   0  0.308174     NaN
Five    -28.1110  46.14830   0 -0.609146     NaN
Four    -47.9801  46.14830   0 -1.039693     NaN
One     -72.1337  46.14830   0 -1.563086     NaN
Seven   -17.8789  46.14830   0 -0.387423     NaN
Six     -24.4512  46.14830   0 -0.529841     NaN
Three   -53.6375  46.14830   0 -1.162285     NaN
Two     -71.2676  46.14830   0 -1.544317     NaN
Zero    -79.9777  46.14830   0 -1.733058     NaN

我已经try

 summary(fm1) <- gsub("Days", "", summary(fm1))

我收到一个错误,说这个对象fm1是类S4,gsub将无法工作.因此,对于如何从示范摘要术语中排除"天数"一词的任何建议,我们将不胜感激.谢谢.

推荐答案

大多数summary()个函数返回的结果可以在打印前修改. 因此,不要运行summary(fm1),而是使用

summ <- summary(fm1)

并在打印前编辑名字. 您可以查看?summary.lme来查找要编辑的组件的名称,或者只是查看names(summ)并猜测. 这对我很有效:

rownames(summ$tTable) <- sub("Days", "", rownames(summ$tTable))
rownames(summ$corFixed) <- sub("Days", "", rownames(summ$corFixed))
colnames(summ$corFixed) <- sub("Days", "", colnames(summ$corFixed))
summ

R相关问答推荐

为什么当我try 在收件箱中使用合并功能时会出现回收错误?

在(g)子中使用asserable字符

有没有一种方法可以从函数中创建一个值的列表,然后将这些值变成R中的直方图?我一直觉得不行

用derrr在R中查找组间的重复项

当月份额减go 当月份额

如何使用STAT_SUMMARY向ggplot2中的密度图添加垂直线

合并DFS列表并将索引提取为新列

使用R中的正则表达式将一列分割为多列

按多列统计频次

变长向量的矢量化和

按组跨多列创建伪变量

按组和连续id计算日期差

使用ifElse语句在ggploy中设置aes y值

Broom.Mixed::Augment不适用于Sample::分析

网络抓取NBA.com

将仪表板中的值框大小更改为Quarto

从矩阵创建系数图

具有由向量定义的可变步长的序列

根据向量对列表元素进行排序

在R中,有没有什么方法可以根据一列中的多个值来过滤行?