-
Notifications
You must be signed in to change notification settings - Fork 46
Expand file tree
/
Copy pathllms.txt
More file actions
230 lines (183 loc) · 6.26 KB
/
Copy pathllms.txt
File metadata and controls
230 lines (183 loc) · 6.26 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
# ggsignif: Significance Brackets for ‘ggplot2’
[](https://cran.r-project.org/package=ggsignif)
[](https://github.com/const-ae/ggsignif/actions)
[](https://CRAN.R-project.org/package=ggsignif)
[](https://app.codecov.io/gh/const-ae/ggsignif?branch=main)
[](https://lifecycle.r-lib.org/articles/stages.html)
## Introduction
This package provides an easy way to indicate if two groups are
significantly different. Commonly this is shown by a bar on top
connecting the groups of interest which itself is annotated with the
level of significance (NS, \*, \*\*, \*\*\*). The package provides a
single layer (`geom_signif`) that takes the groups for comparison and
the test (t.test, wilcox etc.) and adds the annotation to the plot.
## Citation
If you wish to cite this package in a publication, you can run the
following command in your R console:
``` r
citation("ggsignif")
#> To cite 'ggsignif' in publications use:
#>
#> Ahlmann-Eltze, C., & Patil, I. (2021). ggsignif: R Package for
#> Displaying Significance Brackets for 'ggplot2'. PsyArxiv.
#> doi:10.31234/osf.io/7awm6
#>
#> A BibTeX entry for LaTeX users is
#>
#> @Article{,
#> title = {{ggsignif}: R Package for Displaying Significance Brackets for {'ggplot2'}},
#> author = {Ahlmann-Eltze Constantin and Indrajeet Patil},
#> year = {2021},
#> journal = {PsyArxiv},
#> url = {https://psyarxiv.com/7awm6},
#> doi = {10.31234/osf.io/7awm6},
#> }
```
## Example
You can first install this package from `CRAN`:
``` r
install.packages("ggsignif")
```
Or get the latest development version:
``` r
install.packages("remotes")
remotes::install_github("const-ae/ggsignif")
```
Plot significance
``` r
library(ggplot2)
library(ggsignif)
p1 <- ggplot(mpg, aes(class, hwy)) +
geom_boxplot() +
geom_signif(
comparisons = list(c("compact", "midsize"), c("minivan", "suv")),
map_signif_level = TRUE, textsize = 6
) +
ylim(NA, 48)
p1
```

Control the direction (either `x` or `y`) via `orientation`
``` r
p2 <- ggplot(
data = mpg,
mapping = aes(
x = hwy,
y = class
)
) +
geom_boxplot(
orientation = "y"
) +
geom_signif(
comparisons = list(
c("compact", "midsize"),
c("minivan", "suv")
),
map_signif_level = TRUE,
textsize = 6,
margin_top = 0.08,
step_increase = 0.05,
tip_length = 0.01,
orientation = "y"
)
p2
```

Compatible with coord_flip
``` r
p1 + coord_flip()
```

Setting the precise location
This is important if you use `position="dodge"`, because in that case I
cannot calculate the correct position of the bars automatically.
``` r
# Calculate annotation
anno <- t.test(
iris[iris$Petal.Width > 1 & iris$Species == "versicolor", "Sepal.Width"],
iris[iris$Species == "virginica", "Sepal.Width"]
)$p.value
# Make plot with custom x and y position of the bracket
ggplot(iris, aes(x = Species, y = Sepal.Width, fill = Petal.Width > 1)) +
geom_boxplot(position = "dodge") +
geom_signif(
annotation = formatC(anno, digits = 1),
y_position = 4.05, xmin = 2.2, xmax = 3,
tip_length = c(0.2, 0.04)
)
```

`ggsignif` is compatible with facetting (`facet_wrap` or `facet_grid`).
The significance label is calculated for each facet where the axis
labels listed in `comparisons` occur. Note that `ggsignif` fails to
calculate the significance if the data is grouped globally (e.g., by
setting `color`, `fill`, or `group` in `ggplot(aes(...))`). It is fine
to group the data per geom (e.g., set the fill within
`geom_boxplot(aes(fill = ...)))`.
``` r
ggplot(diamonds, aes(x = cut, y = carat)) +
geom_boxplot(aes(fill = color)) +
geom_signif(comparisons = list(
c("Fair", "Good"),
c("Very Good", "Ideal")
)) +
facet_wrap(~color) +
ylim(NA, 6.3)
```

## Advanced Example
Sometimes one needs to have a very fine tuned ability to set the
location of the the significance bars in combination with `facet_wrap`
or `facet_grid`. In those cases it you can set the flag `manual=TRUE`
and provide the annotations as a data.frame:
``` r
annotation_df <- data.frame(
color = c("E", "H"),
start = c("Good", "Fair"),
end = c("Very Good", "Good"),
y = c(3.6, 4.7),
label = c("Comp. 1", "Comp. 2")
)
annotation_df
#> color start end y label
#> 1 E Good Very Good 3.6 Comp. 1
#> 2 H Fair Good 4.7 Comp. 2
ggplot(diamonds, aes(x = cut, y = carat)) +
geom_boxplot() +
geom_signif(
data = annotation_df,
aes(xmin = start, xmax = end, annotations = label, y_position = y),
textsize = 3, vjust = -0.2,
manual = TRUE
) +
facet_wrap(~color) +
ylim(NA, 5.3)
```

You can ignore the warning about the missing aesthetics.
For further details, see:
<https://const-ae.github.io/ggsignif/articles/intro.html>
## Maintenance
This package is provided as is and we currently don’t have any plans and
the capacity to add any new features to it. If there is nonetheless a
feature which you would like to see in the package, you are always
welcome to submit pull request, which we will try to address as soon as
possible.
## Code of Conduct
Please note that the `ggsignif` project is released with a [Contributor
Code of
Conduct](https://const-ae.github.io/ggsignif/CODE_OF_CONDUCT.html). By
contributing to this project, you agree to abide by its terms.
# Package index
## All functions
- [`stat_signif()`](https://const-ae.github.io/ggsignif/reference/stat_signif.md)
[`geom_signif()`](https://const-ae.github.io/ggsignif/reference/stat_signif.md)
: Create significance layer
# Articles
### All vignettes
- [Introduction to 'ggsignif'
package](https://const-ae.github.io/ggsignif/articles/intro.md):