Я использую функцию tableby из пакета arsenal для создания сводных таблиц. Для большей части статистики, которую мне нужно сгенерировать, этот пакет дает мне именно тот формат, который мне задают, за исключением одного. Мне нужно получить в той же ячейке что-то вроде этого:
n (%) [95% CI процента]
Пока я использую countpct функция, которая дает мне "n (%)" и binomCI, которая дает мне пропорцию с 95% -ным доверительным интервалом, но она удваивает число строк в моей финальной таблице, поэтому она не идеальна ...
Как я могу сделать так, чтобы все было на одной линии?
Я пытался выяснить, смогу ли я создать другую функцию из исходных, но я не совсем понимаю их синтаксис ...
Спасибо за вашу помощь.
РЕДАКТИРОВАТЬ: Вот воспроизводимый пример.
Код для исходных функций можно найти здесь .
Так вот, что у меня сейчас:
data<-NULL
data$Visit2<-c(rep("Responder",121),rep("Not Responder",29),rep("Responder",4),rep("Not Responder",47))
data$Group<-c(rep("Tx",150),rep("No Tx",51))
data<-as.data.frame(data)
library(arsenal)
my_controls <- tableby.control(test = F,total = F, cat.stats = c("countpct" ,"binomCI"), conf.level = 0.95)
summary(tableby(Group ~ Visit2,
data = data,
control = my_controls),
digits=2, digits.p=3, digits.pct=1)
# Results :
| | No Tx (N=51) | Tx (N=150) |
|:-------------------------------|:-----------------:|:-----------------:|
|**Visit2** | | |
| Not Responder | 47 (92.2%) | 29 (19.3%) |
| Responder | 4 (7.8%) | 121 (80.7%) |
| Not Responder | 0.92 (0.81, 0.98) | 0.19 (0.13, 0.27) |
| Responder | 0.08 (0.02, 0.19) | 0.81 (0.73, 0.87) |
And this is what I want :
| | No Tx (N=51) | Tx (N=150) |
|:----------------|:-------------------------:|:------------------------:|
|**Visit2** | |
| Not Responder | 47 (92.2%) [81.1, 97.8] | 29 (19.3%) [13.3, 26.6] |
| Responder | 4 (7.8%) [2.2, 18.9] | 121 (80.7%) [73.4, 86.7] |
|