The function summaryByVisit() produces a table with descriptive statistics for continues and categorical variables at different time points (visits). It is largely based on the function gtsummary::tbl_summary() and gtsummary::tbl_strata_nested_stack. The changes as compared to these functions are:

A hierarchical summary table is created by variable (lab-value), visit and visitgroup.

Setup and data

To demonstrate the various functionalities of the function, we will create a small dataset. We have three hypotetical lab-values: LDH, ANC and Lymphocytes and two categorical variables, grade and stage. Each lab-value has been measured at 10 time points (visits), which can further be divided into three visit groups (baseline, treatment and follow-up). In addition, we have a grouping variable called ‘arm’.

data<-NULL
visit <- c(paste0(rep("Visit ", 10), rbind(c(1:10))),
           paste0(rep("Visit ", 10), rbind(c(1:10))),
           paste0(rep("Visit ", 10), rbind(c(1:10))))
data <- as.data.frame(cbind( visit, rnorm(30)))
data<-as.data.frame(rbind(data, data, data, data, data))
data$visitgroup<- ifelse(data$visit %in% c("Visit 1", "Visit 2"), "Baseline", ifelse(data$visit %in% c("Visit 3", "Visit 4"), "Treatment", "Follow-up"))
data$visitgroup<-factor(data$visitgroup, levels = c("Baseline", "Treatment", "Follow-up"))
data$LDH<-rnorm(150)
data$Lymphocytes<-rnorm(150)
data$ANC<-rnorm(150)
data$LDH[3]<-NA
data$stage <- c(rep("II", 75), rep("I", 75))
data$stage <- ifelse(data$Lymphocytes>1.5, NA, data$stage)
data$grade <- c(rep("IV", 50), rep("III", 50), rep("II", 25), rep("I", 25))
data$arm<- c(rep("Arm A", 70), rep("Arm B", 80))

Basic table

Now, we use summarySCI::summaryByVisit to create a summary table for median and range of each lab-value by visit.

summaryByVisit(data = data,
         vars = c("LDH", "Lymphocytes", "stage"),
         visit = "visit", 
         add_n = TRUE)

Visit	N1	N = 152
LDH
Visit 1	15	0.3 (-1.0, 2.3)
Visit 2	15	0.2 (-1.6, 2.6)
Visit 3	14	0.4 (-1.9, 1.9)
Visit 4	15	-0.3 (-1.4, 1.5)
Visit 5	15	0.0 (-2.0, 0.9)
Visit 6	15	-0.2 (-1.5, 1.8)
Visit 7	15	-0.1 (-1.8, 0.8)
Visit 8	15	-0.3 (-2.8, 1.5)
Visit 9	15	-0.2 (-2.6, 2.7)
Visit 10	15	-0.1 (-2.3, 1.6)
Lymphocytes
Visit 1	15	0.3 (-2.2, 1.7)
Visit 2	15	0.1 (-1.8, 1.2)
Visit 3	15	0.5 (-2.2, 2.2)
Visit 4	15	-0.4 (-2.8, 2.5)
Visit 5	15	-0.5 (-1.4, 0.6)
Visit 6	15	-0.2 (-1.7, 1.9)
Visit 7	15	0.1 (-0.8, 1.8)
Visit 8	15	0.1 (-2.4, 1.2)
Visit 9	15	-0.2 (-2.8, 1.3)
Visit 10	15	-0.2 (-1.6, 2.2)
stage
Visit 1	14
I		6.0 (40.0%)
II		8.0 (53.3%)
Missing		1.0 (6.7%)
Visit 2	15
I		7.0 (46.7%)
II		8.0 (53.3%)
Visit 3	13
I		6.0 (40.0%)
II		7.0 (46.7%)
Missing		2.0 (13.3%)
Visit 4	14
I		6.0 (40.0%)
II		8.0 (53.3%)
Missing		1.0 (6.7%)
Visit 5	15
I		7.0 (46.7%)
II		8.0 (53.3%)
Visit 6	14
I		7.0 (46.7%)
II		7.0 (46.7%)
Missing		1.0 (6.7%)
Visit 7	13
I		6.0 (40.0%)
II		7.0 (46.7%)
Missing		2.0 (13.3%)
Visit 8	15
I		8.0 (53.3%)
II		7.0 (46.7%)
Visit 9	15
I		8.0 (53.3%)
II		7.0 (46.7%)
Visit 10	13
I		7.0 (46.7%)
II		6.0 (40.0%)
Missing		2.0 (13.3%)
1N without missing values
2Median (Min, Max); n (%)

Add visit group

We can the visit group to group visits accordingly. Visitgroup needs to be an ordered factor.

summaryByVisit(data,
         vars = c("LDH", "Lymphocytes", "stage"),
         visitgroup = "visitgroup",
         visit = "visit")

Visit	N = 151
LDH
Baseline
Visit 1	0.3 (-1.0, 2.3)
Visit 2	0.2 (-1.6, 2.6)
Treatment
Visit 3	0.4 (-1.9, 1.9)
Visit 4	-0.3 (-1.4, 1.5)
Follow-up
Visit 5	0.0 (-2.0, 0.9)
Visit 6	-0.2 (-1.5, 1.8)
Visit 7	-0.1 (-1.8, 0.8)
Visit 8	-0.3 (-2.8, 1.5)
Visit 9	-0.2 (-2.6, 2.7)
Visit 10	-0.1 (-2.3, 1.6)
Lymphocytes
Baseline
Visit 1	0.3 (-2.2, 1.7)
Visit 2	0.1 (-1.8, 1.2)
Treatment
Visit 3	0.5 (-2.2, 2.2)
Visit 4	-0.4 (-2.8, 2.5)
Follow-up
Visit 5	-0.5 (-1.4, 0.6)
Visit 6	-0.2 (-1.7, 1.9)
Visit 7	0.1 (-0.8, 1.8)
Visit 8	0.1 (-2.4, 1.2)
Visit 9	-0.2 (-2.8, 1.3)
Visit 10	-0.2 (-1.6, 2.2)
stage
Baseline
Visit 1
I	6.0 (40.0%)
II	8.0 (53.3%)
Missing	1.0 (6.7%)
Visit 2
I	7.0 (46.7%)
II	8.0 (53.3%)
Treatment
Visit 3
I	6.0 (40.0%)
II	7.0 (46.7%)
Missing	2.0 (13.3%)
Visit 4
I	6.0 (40.0%)
II	8.0 (53.3%)
Missing	1.0 (6.7%)
Follow-up
Visit 5
I	7.0 (46.7%)
II	8.0 (53.3%)
Visit 6
I	7.0 (46.7%)
II	7.0 (46.7%)
Missing	1.0 (6.7%)
Visit 7
I	6.0 (40.0%)
II	7.0 (46.7%)
Missing	2.0 (13.3%)
Visit 8
I	8.0 (53.3%)
II	7.0 (46.7%)
Visit 9
I	8.0 (53.3%)
II	7.0 (46.7%)
Visit 10
I	7.0 (46.7%)
II	6.0 (40.0%)
Missing	2.0 (13.3%)
1Median (Min, Max); n (%)

By group

We can stratify the table by groups via the group argument. The overall column can still be shown if desired, using the overall = TRUE argument. A maximum of 3 groups are supported.

summaryByVisit(data = data,
         vars = c("LDH", "Lymphocytes", "stage"),
         group = "arm",
         visitgroup = "visitgroup",
         visit = "visit")

Visit	Arm A N = 71	Arm B N = 81
LDH
Baseline
Visit 1	0.3 (0.0, 1.8)	0.1 (-1.0, 2.3)
Visit 2	-0.5 (-1.6, 0.4)	0.4 (-0.5, 2.6)
Treatment
Visit 3	0.1 (-1.4, 0.8)	0.4 (-1.9, 1.9)
Visit 4	-0.3 (-1.0, 0.3)	0.1 (-1.4, 1.5)
Follow-up
Visit 5	0.2 (-1.2, 0.9)	-0.1 (-2.0, 0.5)
Visit 6	0.4 (-0.9, 1.8)	-0.5 (-1.5, 0.9)
Visit 7	-0.1 (-0.9, 0.3)	-0.3 (-1.8, 0.8)
Visit 8	-0.4 (-1.9, 1.2)	0.0 (-2.8, 1.5)
Visit 9	-1.0 (-2.6, 2.0)	0.3 (-1.2, 2.7)
Visit 10	-0.4 (-2.3, 1.6)	0.1 (-1.9, 1.5)
Lymphocytes
Baseline
Visit 1	0.4 (-0.6, 1.4)	0.3 (-2.2, 1.7)
Visit 2	-0.5 (-1.1, 1.2)	0.7 (-1.8, 1.2)
Treatment
Visit 3	0.7 (-0.1, 2.2)	-0.1 (-2.2, 2.0)
Visit 4	-0.4 (-2.8, 1.5)	-0.3 (-1.6, 2.5)
Follow-up
Visit 5	-0.5 (-1.4, 0.0)	-0.3 (-1.2, 0.6)
Visit 6	0.0 (-1.7, 0.7)	-0.2 (-0.8, 1.9)
Visit 7	0.1 (-0.8, 1.0)	0.3 (-0.7, 1.8)
Visit 8	0.2 (-1.0, 0.8)	-0.2 (-2.4, 1.2)
Visit 9	-0.3 (-2.8, 1.3)	0.3 (-1.0, 1.2)
Visit 10	-0.2 (-1.6, 2.2)	0.0 (-1.5, 1.8)
stage
Baseline
Visit 1
I	0.0 (0.0%)	6.0 (75.0%)
II	7.0 (100.0%)	1.0 (12.5%)
Missing	0.0 (0.0%)	1.0 (12.5%)
Visit 2
I	0.0 (0.0%)	7.0 (87.5%)
II	7.0 (100.0%)	1.0 (12.5%)
Treatment
Visit 3
I	0.0 (0.0%)	6.0 (75.0%)
II	6.0 (85.7%)	1.0 (12.5%)
Missing	1.0 (14.3%)	1.0 (12.5%)
Visit 4
I	0.0 (0.0%)	6.0 (75.0%)
II	7.0 (100.0%)	1.0 (12.5%)
Missing	0.0 (0.0%)	1.0 (12.5%)
Follow-up
Visit 5
I	0.0 (0.0%)	7.0 (87.5%)
II	7.0 (100.0%)	1.0 (12.5%)
Visit 6
I	0.0 (0.0%)	7.0 (87.5%)
II	7.0 (100.0%)	0.0 (0.0%)
Missing	0.0 (0.0%)	1.0 (12.5%)
Visit 7
I	0.0 (0.0%)	6.0 (75.0%)
II	7.0 (100.0%)	0.0 (0.0%)
Missing	0.0 (0.0%)	2.0 (25.0%)
Visit 8
I	0.0 (0.0%)	8.0 (100.0%)
II	7.0 (100.0%)	0.0 (0.0%)
Visit 9
I	0.0 (0.0%)	8.0 (100.0%)
II	7.0 (100.0%)	0.0 (0.0%)
Visit 10
I	0.0 (0.0%)	7.0 (87.5%)
II	6.0 (85.7%)	0.0 (0.0%)
Missing	1.0 (14.3%)	1.0 (12.5%)
1Median (Min, Max); n (%)

summaryByVisit(data,
         vars = c("LDH", "Lymphocytes", "stage"),
         group = "arm",
         overall = TRUE,
         visitgroup = "visitgroup",
         visit = "visit")

Visit	Arm A N = 71	Arm B N = 81	Overall N = 151
LDH
Baseline
Visit 1	0.3 (0.0, 1.8)	0.1 (-1.0, 2.3)	0.3 (-1.0, 2.3)
Visit 2	-0.5 (-1.6, 0.4)	0.4 (-0.5, 2.6)	0.2 (-1.6, 2.6)
Treatment
Visit 3	0.1 (-1.4, 0.8)	0.4 (-1.9, 1.9)	0.4 (-1.9, 1.9)
Visit 4	-0.3 (-1.0, 0.3)	0.1 (-1.4, 1.5)	-0.3 (-1.4, 1.5)
Follow-up
Visit 5	0.2 (-1.2, 0.9)	-0.1 (-2.0, 0.5)	0.0 (-2.0, 0.9)
Visit 6	0.4 (-0.9, 1.8)	-0.5 (-1.5, 0.9)	-0.2 (-1.5, 1.8)
Visit 7	-0.1 (-0.9, 0.3)	-0.3 (-1.8, 0.8)	-0.1 (-1.8, 0.8)
Visit 8	-0.4 (-1.9, 1.2)	0.0 (-2.8, 1.5)	-0.3 (-2.8, 1.5)
Visit 9	-1.0 (-2.6, 2.0)	0.3 (-1.2, 2.7)	-0.2 (-2.6, 2.7)
Visit 10	-0.4 (-2.3, 1.6)	0.1 (-1.9, 1.5)	-0.1 (-2.3, 1.6)
Lymphocytes
Baseline
Visit 1	0.4 (-0.6, 1.4)	0.3 (-2.2, 1.7)	0.3 (-2.2, 1.7)
Visit 2	-0.5 (-1.1, 1.2)	0.7 (-1.8, 1.2)	0.1 (-1.8, 1.2)
Treatment
Visit 3	0.7 (-0.1, 2.2)	-0.1 (-2.2, 2.0)	0.5 (-2.2, 2.2)
Visit 4	-0.4 (-2.8, 1.5)	-0.3 (-1.6, 2.5)	-0.4 (-2.8, 2.5)
Follow-up
Visit 5	-0.5 (-1.4, 0.0)	-0.3 (-1.2, 0.6)	-0.5 (-1.4, 0.6)
Visit 6	0.0 (-1.7, 0.7)	-0.2 (-0.8, 1.9)	-0.2 (-1.7, 1.9)
Visit 7	0.1 (-0.8, 1.0)	0.3 (-0.7, 1.8)	0.1 (-0.8, 1.8)
Visit 8	0.2 (-1.0, 0.8)	-0.2 (-2.4, 1.2)	0.1 (-2.4, 1.2)
Visit 9	-0.3 (-2.8, 1.3)	0.3 (-1.0, 1.2)	-0.2 (-2.8, 1.3)
Visit 10	-0.2 (-1.6, 2.2)	0.0 (-1.5, 1.8)	-0.2 (-1.6, 2.2)
stage
Baseline
Visit 1
I	0.0 (0.0%)	6.0 (75.0%)	6.0 (40.0%)
II	7.0 (100.0%)	1.0 (12.5%)	8.0 (53.3%)
Missing	0.0 (0.0%)	1.0 (12.5%)	1.0 (6.7%)
Visit 2
I	0.0 (0.0%)	7.0 (87.5%)	7.0 (46.7%)
II	7.0 (100.0%)	1.0 (12.5%)	8.0 (53.3%)
Treatment
Visit 3
I	0.0 (0.0%)	6.0 (75.0%)	6.0 (40.0%)
II	6.0 (85.7%)	1.0 (12.5%)	7.0 (46.7%)
Missing	1.0 (14.3%)	1.0 (12.5%)	2.0 (13.3%)
Visit 4
I	0.0 (0.0%)	6.0 (75.0%)	6.0 (40.0%)
II	7.0 (100.0%)	1.0 (12.5%)	8.0 (53.3%)
Missing	0.0 (0.0%)	1.0 (12.5%)	1.0 (6.7%)
Follow-up
Visit 5
I	0.0 (0.0%)	7.0 (87.5%)	7.0 (46.7%)
II	7.0 (100.0%)	1.0 (12.5%)	8.0 (53.3%)
Visit 6
I	0.0 (0.0%)	7.0 (87.5%)	7.0 (46.7%)
II	7.0 (100.0%)	0.0 (0.0%)	7.0 (46.7%)
Missing	0.0 (0.0%)	1.0 (12.5%)	1.0 (6.7%)
Visit 7
I	0.0 (0.0%)	6.0 (75.0%)	6.0 (40.0%)
II	7.0 (100.0%)	0.0 (0.0%)	7.0 (46.7%)
Missing	0.0 (0.0%)	2.0 (25.0%)	2.0 (13.3%)
Visit 8
I	0.0 (0.0%)	8.0 (100.0%)	8.0 (53.3%)
II	7.0 (100.0%)	0.0 (0.0%)	7.0 (46.7%)
Visit 9
I	0.0 (0.0%)	8.0 (100.0%)	8.0 (53.3%)
II	7.0 (100.0%)	0.0 (0.0%)	7.0 (46.7%)
Visit 10
I	0.0 (0.0%)	7.0 (87.5%)	7.0 (46.7%)
II	6.0 (85.7%)	0.0 (0.0%)	6.0 (40.0%)
Missing	1.0 (14.3%)	1.0 (12.5%)	2.0 (13.3%)
1Median (Min, Max); n (%)

Add N

Sample size can be shown for each column, if the option add_n is set to TRUE.

summaryByVisit(data,
         vars = c("LDH", "Lymphocytes", "stage"),
         group = "arm",
         overall = TRUE,
         visitgroup = "visitgroup",
         visit = "visit",
         add_n = TRUE)

Visit	N1	Arm A N = 72	N1	Arm B N = 82	N1	Overall N = 152
LDH
Baseline
Visit 1	7	0.3 (0.0, 1.8)	8	0.1 (-1.0, 2.3)	15	0.3 (-1.0, 2.3)
Visit 2	7	-0.5 (-1.6, 0.4)	8	0.4 (-0.5, 2.6)	15	0.2 (-1.6, 2.6)
Treatment
Visit 3	6	0.1 (-1.4, 0.8)	8	0.4 (-1.9, 1.9)	14	0.4 (-1.9, 1.9)
Visit 4	7	-0.3 (-1.0, 0.3)	8	0.1 (-1.4, 1.5)	15	-0.3 (-1.4, 1.5)
Follow-up
Visit 5	7	0.2 (-1.2, 0.9)	8	-0.1 (-2.0, 0.5)	15	0.0 (-2.0, 0.9)
Visit 6	7	0.4 (-0.9, 1.8)	8	-0.5 (-1.5, 0.9)	15	-0.2 (-1.5, 1.8)
Visit 7	7	-0.1 (-0.9, 0.3)	8	-0.3 (-1.8, 0.8)	15	-0.1 (-1.8, 0.8)
Visit 8	7	-0.4 (-1.9, 1.2)	8	0.0 (-2.8, 1.5)	15	-0.3 (-2.8, 1.5)
Visit 9	7	-1.0 (-2.6, 2.0)	8	0.3 (-1.2, 2.7)	15	-0.2 (-2.6, 2.7)
Visit 10	7	-0.4 (-2.3, 1.6)	8	0.1 (-1.9, 1.5)	15	-0.1 (-2.3, 1.6)
Lymphocytes
Baseline
Visit 1	7	0.4 (-0.6, 1.4)	8	0.3 (-2.2, 1.7)	15	0.3 (-2.2, 1.7)
Visit 2	7	-0.5 (-1.1, 1.2)	8	0.7 (-1.8, 1.2)	15	0.1 (-1.8, 1.2)
Treatment
Visit 3	7	0.7 (-0.1, 2.2)	8	-0.1 (-2.2, 2.0)	15	0.5 (-2.2, 2.2)
Visit 4	7	-0.4 (-2.8, 1.5)	8	-0.3 (-1.6, 2.5)	15	-0.4 (-2.8, 2.5)
Follow-up
Visit 5	7	-0.5 (-1.4, 0.0)	8	-0.3 (-1.2, 0.6)	15	-0.5 (-1.4, 0.6)
Visit 6	7	0.0 (-1.7, 0.7)	8	-0.2 (-0.8, 1.9)	15	-0.2 (-1.7, 1.9)
Visit 7	7	0.1 (-0.8, 1.0)	8	0.3 (-0.7, 1.8)	15	0.1 (-0.8, 1.8)
Visit 8	7	0.2 (-1.0, 0.8)	8	-0.2 (-2.4, 1.2)	15	0.1 (-2.4, 1.2)
Visit 9	7	-0.3 (-2.8, 1.3)	8	0.3 (-1.0, 1.2)	15	-0.2 (-2.8, 1.3)
Visit 10	7	-0.2 (-1.6, 2.2)	8	0.0 (-1.5, 1.8)	15	-0.2 (-1.6, 2.2)
stage
Baseline
Visit 1	7		7		14
I		0.0 (0.0%)		6.0 (75.0%)		6.0 (40.0%)
II		7.0 (100.0%)		1.0 (12.5%)		8.0 (53.3%)
Missing		0.0 (0.0%)		1.0 (12.5%)		1.0 (6.7%)
Visit 2	7		8		15
I		0.0 (0.0%)		7.0 (87.5%)		7.0 (46.7%)
II		7.0 (100.0%)		1.0 (12.5%)		8.0 (53.3%)
Treatment
Visit 3	6		7		13
I		0.0 (0.0%)		6.0 (75.0%)		6.0 (40.0%)
II		6.0 (85.7%)		1.0 (12.5%)		7.0 (46.7%)
Missing		1.0 (14.3%)		1.0 (12.5%)		2.0 (13.3%)
Visit 4	7		7		14
I		0.0 (0.0%)		6.0 (75.0%)		6.0 (40.0%)
II		7.0 (100.0%)		1.0 (12.5%)		8.0 (53.3%)
Missing		0.0 (0.0%)		1.0 (12.5%)		1.0 (6.7%)
Follow-up
Visit 5	7		8		15
I		0.0 (0.0%)		7.0 (87.5%)		7.0 (46.7%)
II		7.0 (100.0%)		1.0 (12.5%)		8.0 (53.3%)
Visit 6	7		7		14
I		0.0 (0.0%)		7.0 (87.5%)		7.0 (46.7%)
II		7.0 (100.0%)		0.0 (0.0%)		7.0 (46.7%)
Missing		0.0 (0.0%)		1.0 (12.5%)		1.0 (6.7%)
Visit 7	7		6		13
I		0.0 (0.0%)		6.0 (75.0%)		6.0 (40.0%)
II		7.0 (100.0%)		0.0 (0.0%)		7.0 (46.7%)
Missing		0.0 (0.0%)		2.0 (25.0%)		2.0 (13.3%)
Visit 8	7		8		15
I		0.0 (0.0%)		8.0 (100.0%)		8.0 (53.3%)
II		7.0 (100.0%)		0.0 (0.0%)		7.0 (46.7%)
Visit 9	7		8		15
I		0.0 (0.0%)		8.0 (100.0%)		8.0 (53.3%)
II		7.0 (100.0%)		0.0 (0.0%)		7.0 (46.7%)
Visit 10	6		7		13
I		0.0 (0.0%)		7.0 (87.5%)		7.0 (46.7%)
II		6.0 (85.7%)		0.0 (0.0%)		6.0 (40.0%)
Missing		1.0 (14.3%)		1.0 (12.5%)		2.0 (13.3%)
1N without missing values
2Median (Min, Max); n (%)

Change labels

The labels can be changed by using the “labels” statement.

summaryByVisit(data,
         vars = c("LDH", "Lymphocytes", "stage"),
         labels=list(LDH="LDH (g/L)", Lymphocytes="Lymphocytes (g)", stage="Stage"),
         group = "arm",
         overall = TRUE,
         visitgroup = "visitgroup",
         visit = "visit",
         add_n = TRUE)

Visit	N1	Arm A N = 72	N1	Arm B N = 82	N1	Overall N = 152
LDH (g/L)
Baseline
Visit 1	7	0.3 (0.0, 1.8)	8	0.1 (-1.0, 2.3)	15	0.3 (-1.0, 2.3)
Visit 2	7	-0.5 (-1.6, 0.4)	8	0.4 (-0.5, 2.6)	15	0.2 (-1.6, 2.6)
Treatment
Visit 3	6	0.1 (-1.4, 0.8)	8	0.4 (-1.9, 1.9)	14	0.4 (-1.9, 1.9)
Visit 4	7	-0.3 (-1.0, 0.3)	8	0.1 (-1.4, 1.5)	15	-0.3 (-1.4, 1.5)
Follow-up
Visit 5	7	0.2 (-1.2, 0.9)	8	-0.1 (-2.0, 0.5)	15	0.0 (-2.0, 0.9)
Visit 6	7	0.4 (-0.9, 1.8)	8	-0.5 (-1.5, 0.9)	15	-0.2 (-1.5, 1.8)
Visit 7	7	-0.1 (-0.9, 0.3)	8	-0.3 (-1.8, 0.8)	15	-0.1 (-1.8, 0.8)
Visit 8	7	-0.4 (-1.9, 1.2)	8	0.0 (-2.8, 1.5)	15	-0.3 (-2.8, 1.5)
Visit 9	7	-1.0 (-2.6, 2.0)	8	0.3 (-1.2, 2.7)	15	-0.2 (-2.6, 2.7)
Visit 10	7	-0.4 (-2.3, 1.6)	8	0.1 (-1.9, 1.5)	15	-0.1 (-2.3, 1.6)
Lymphocytes (g)
Baseline
Visit 1	7	0.4 (-0.6, 1.4)	8	0.3 (-2.2, 1.7)	15	0.3 (-2.2, 1.7)
Visit 2	7	-0.5 (-1.1, 1.2)	8	0.7 (-1.8, 1.2)	15	0.1 (-1.8, 1.2)
Treatment
Visit 3	7	0.7 (-0.1, 2.2)	8	-0.1 (-2.2, 2.0)	15	0.5 (-2.2, 2.2)
Visit 4	7	-0.4 (-2.8, 1.5)	8	-0.3 (-1.6, 2.5)	15	-0.4 (-2.8, 2.5)
Follow-up
Visit 5	7	-0.5 (-1.4, 0.0)	8	-0.3 (-1.2, 0.6)	15	-0.5 (-1.4, 0.6)
Visit 6	7	0.0 (-1.7, 0.7)	8	-0.2 (-0.8, 1.9)	15	-0.2 (-1.7, 1.9)
Visit 7	7	0.1 (-0.8, 1.0)	8	0.3 (-0.7, 1.8)	15	0.1 (-0.8, 1.8)
Visit 8	7	0.2 (-1.0, 0.8)	8	-0.2 (-2.4, 1.2)	15	0.1 (-2.4, 1.2)
Visit 9	7	-0.3 (-2.8, 1.3)	8	0.3 (-1.0, 1.2)	15	-0.2 (-2.8, 1.3)
Visit 10	7	-0.2 (-1.6, 2.2)	8	0.0 (-1.5, 1.8)	15	-0.2 (-1.6, 2.2)
Stage
Baseline
Visit 1	7		7		14
I		0.0 (0.0%)		6.0 (75.0%)		6.0 (40.0%)
II		7.0 (100.0%)		1.0 (12.5%)		8.0 (53.3%)
Missing		0.0 (0.0%)		1.0 (12.5%)		1.0 (6.7%)
Visit 2	7		8		15
I		0.0 (0.0%)		7.0 (87.5%)		7.0 (46.7%)
II		7.0 (100.0%)		1.0 (12.5%)		8.0 (53.3%)
Treatment
Visit 3	6		7		13
I		0.0 (0.0%)		6.0 (75.0%)		6.0 (40.0%)
II		6.0 (85.7%)		1.0 (12.5%)		7.0 (46.7%)
Missing		1.0 (14.3%)		1.0 (12.5%)		2.0 (13.3%)
Visit 4	7		7		14
I		0.0 (0.0%)		6.0 (75.0%)		6.0 (40.0%)
II		7.0 (100.0%)		1.0 (12.5%)		8.0 (53.3%)
Missing		0.0 (0.0%)		1.0 (12.5%)		1.0 (6.7%)
Follow-up
Visit 5	7		8		15
I		0.0 (0.0%)		7.0 (87.5%)		7.0 (46.7%)
II		7.0 (100.0%)		1.0 (12.5%)		8.0 (53.3%)
Visit 6	7		7		14
I		0.0 (0.0%)		7.0 (87.5%)		7.0 (46.7%)
II		7.0 (100.0%)		0.0 (0.0%)		7.0 (46.7%)
Missing		0.0 (0.0%)		1.0 (12.5%)		1.0 (6.7%)
Visit 7	7		6		13
I		0.0 (0.0%)		6.0 (75.0%)		6.0 (40.0%)
II		7.0 (100.0%)		0.0 (0.0%)		7.0 (46.7%)
Missing		0.0 (0.0%)		2.0 (25.0%)		2.0 (13.3%)
Visit 8	7		8		15
I		0.0 (0.0%)		8.0 (100.0%)		8.0 (53.3%)
II		7.0 (100.0%)		0.0 (0.0%)		7.0 (46.7%)
Visit 9	7		8		15
I		0.0 (0.0%)		8.0 (100.0%)		8.0 (53.3%)
II		7.0 (100.0%)		0.0 (0.0%)		7.0 (46.7%)
Visit 10	6		7		13
I		0.0 (0.0%)		7.0 (87.5%)		7.0 (46.7%)
II		6.0 (85.7%)		0.0 (0.0%)		6.0 (40.0%)
Missing		1.0 (14.3%)		1.0 (12.5%)		2.0 (13.3%)
1N without missing values
2Median (Min, Max); n (%)

Further steps:

improve costumization