I do not see anything obviously wrong but I think we need to see all of your code and some sample data. See
FAQ Asking Questions
A handy way to supply some sample data is the dput() function. In the case of a large dataset something like dput(head(mydata, 100)) should supply the data we need. Just do dput(mydata) where mydata is your data. Copy the output and paste it here between
```
What does drg2018_s50_fusion look like?