Subtract Signatures/Catalogues

Subtracts signature2 from signature1 and returns result.

Usage

sig_subtract(signature1, signature2)

signature1 %-% signature2

Arguments

signature1, signature2

sigverse signature/catalogue data.frames. See sigshared::example_signature() or sigshared::example_catalogue().

Value

a data.frame representing a sigverse signature.

Examples

library(sigstash)

# Load a signature collection
signatures <- sig_load("COSMIC_v3.3.1_SBS_GRCh38")

# Subtract signatures
signatures[['SBS3']] %-% signatures[['SBS4']]
#>     type channel      fraction
#> 193  C>A A[C>A]A -2.153099e-02
#> 194  C>A A[C>A]C -1.664639e-02
#> 195  C>A A[C>A]G -1.430841e-02
#> 196  C>A A[C>A]T -1.739770e-02
#> 197  C>G A[C>G]A  1.288248e-02
#> 198  C>G A[C>G]C  8.775283e-03
#> 199  C>G A[C>G]G -1.061015e-03
#> 200  C>G A[C>G]T  1.392321e-02
#> 201  C>T A[C>T]A  5.530060e-03
#> 202  C>T A[C>T]C  8.141122e-03
#> 203  C>T A[C>T]G  1.845617e-03
#> 204  C>T A[C>T]T  7.891082e-03
#> 205  T>A A[T>A]A -3.918470e-03
#> 206  T>A A[T>A]C  3.694537e-03
#> 207  T>A A[T>A]G -2.753820e-03
#> 208  T>A A[T>A]T  1.894499e-03
#> 209  T>C A[T>C]A  8.519717e-03
#> 210  T>C A[T>C]C  7.073942e-03
#> 211  T>C A[T>C]G  7.276744e-03
#> 212  T>C A[T>C]T  1.570053e-02
#> 213  T>G A[T>G]A  2.890030e-03
#> 214  T>G A[T>G]C  2.451044e-03
#> 215  T>G A[T>G]G  4.864805e-03
#> 216  T>G A[T>G]T  3.744614e-03
#> 217  C>A C[C>A]A -5.789319e-02
#> 218  C>A C[C>A]C -5.375123e-02
#> 219  C>A C[C>A]G -2.207407e-02
#> 220  C>A C[C>A]T -5.331806e-02
#> 221  C>G C[C>G]A  1.140520e-02
#> 222  C>G C[C>G]C  7.140612e-03
#> 223  C>G C[C>G]G -1.907500e-03
#> 224  C>G C[C>G]T  1.539370e-02
#> 225  C>T C[C>T]A  3.378900e-03
#> 226  C>T C[C>T]C  8.591852e-03
#> 227  C>T C[C>T]G -1.074480e-03
#> 228  C>T C[C>T]T  6.958618e-03
#> 229  T>A C[T>A]A -6.461466e-03
#> 230  T>A C[T>A]C  5.613680e-03
#> 231  T>A C[T>A]G -1.930279e-02
#> 232  T>A C[T>A]T  9.035457e-04
#> 233  T>C C[T>C]A  3.648386e-03
#> 234  T>C C[T>C]C  1.217409e-02
#> 235  T>C C[T>C]G  4.598236e-03
#> 236  T>C C[T>C]T  1.332373e-02
#> 237  T>G C[T>G]A  3.499204e-03
#> 238  T>G C[T>G]C  6.617684e-03
#> 239  T>G C[T>G]G  6.908170e-03
#> 240  T>G C[T>G]T  5.866854e-03
#> 241  C>A G[C>A]A -2.452313e-02
#> 242  C>A G[C>A]C -2.362341e-02
#> 243  C>A G[C>A]G -1.400985e-02
#> 244  C>A G[C>A]T -1.116893e-02
#> 245  C>G G[C>G]A  7.423360e-03
#> 246  C>G G[C>G]C  2.844868e-03
#> 247  C>G G[C>G]G -8.893042e-04
#> 248  C>G G[C>G]T  1.162307e-02
#> 249  C>T G[C>T]A  1.233306e-02
#> 250  C>T G[C>T]C  8.833372e-03
#> 251  C>T G[C>T]G  4.636514e-05
#> 252  C>T G[C>T]T  7.465936e-03
#> 253  T>A G[T>A]A -2.148785e-03
#> 254  T>A G[T>A]C  3.864925e-03
#> 255  T>A G[T>A]G -1.509206e-03
#> 256  T>A G[T>A]T  6.876089e-03
#> 257  T>C G[T>C]A  7.037202e-03
#> 258  T>C G[T>C]C  4.443212e-03
#> 259  T>C G[T>C]G  8.568192e-03
#> 260  T>C G[T>C]T  9.580148e-03
#> 261  T>G G[T>G]A  3.799073e-03
#> 262  T>G G[T>G]C  2.079199e-03
#> 263  T>G G[T>G]G  8.507909e-03
#> 264  T>G G[T>G]T  5.657852e-03
#> 265  C>A T[C>A]A -3.158094e-02
#> 266  C>A T[C>A]C -3.022953e-02
#> 267  C>A T[C>A]G -9.694986e-03
#> 268  C>A T[C>A]T -3.165580e-02
#> 269  C>G T[C>G]A  9.776192e-03
#> 270  C>G T[C>G]C  1.440654e-02
#> 271  C>G T[C>G]G -1.338357e-03
#> 272  C>G T[C>G]T  2.032369e-02
#> 273  C>T T[C>T]A  2.031125e-03
#> 274  C>T T[C>T]C  6.334056e-03
#> 275  C>T T[C>T]G  7.230880e-04
#> 276  C>T T[C>T]T  9.583094e-04
#> 277  T>A T[T>A]A -2.342574e-03
#> 278  T>A T[T>A]C  6.308687e-03
#> 279  T>A T[T>A]G -2.853877e-03
#> 280  T>A T[T>A]T  7.409249e-03
#> 281  T>C T[T>C]A  1.263288e-02
#> 282  T>C T[T>C]C  8.650741e-03
#> 283  T>C T[T>C]G  5.230234e-03
#> 284  T>C T[T>C]T  1.331831e-02
#> 285  T>G T[T>G]A  6.788732e-03
#> 286  T>G T[T>G]C  6.223839e-03
#> 287  T>G T[T>G]G  5.774346e-03
#> 288  T>G T[T>G]T  9.972652e-03

# Identical approach using full function name
sig_subtract(signatures[['SBS3']], signatures[['SBS4']])
#>     type channel      fraction
#> 193  C>A A[C>A]A -2.153099e-02
#> 194  C>A A[C>A]C -1.664639e-02
#> 195  C>A A[C>A]G -1.430841e-02
#> 196  C>A A[C>A]T -1.739770e-02
#> 197  C>G A[C>G]A  1.288248e-02
#> 198  C>G A[C>G]C  8.775283e-03
#> 199  C>G A[C>G]G -1.061015e-03
#> 200  C>G A[C>G]T  1.392321e-02
#> 201  C>T A[C>T]A  5.530060e-03
#> 202  C>T A[C>T]C  8.141122e-03
#> 203  C>T A[C>T]G  1.845617e-03
#> 204  C>T A[C>T]T  7.891082e-03
#> 205  T>A A[T>A]A -3.918470e-03
#> 206  T>A A[T>A]C  3.694537e-03
#> 207  T>A A[T>A]G -2.753820e-03
#> 208  T>A A[T>A]T  1.894499e-03
#> 209  T>C A[T>C]A  8.519717e-03
#> 210  T>C A[T>C]C  7.073942e-03
#> 211  T>C A[T>C]G  7.276744e-03
#> 212  T>C A[T>C]T  1.570053e-02
#> 213  T>G A[T>G]A  2.890030e-03
#> 214  T>G A[T>G]C  2.451044e-03
#> 215  T>G A[T>G]G  4.864805e-03
#> 216  T>G A[T>G]T  3.744614e-03
#> 217  C>A C[C>A]A -5.789319e-02
#> 218  C>A C[C>A]C -5.375123e-02
#> 219  C>A C[C>A]G -2.207407e-02
#> 220  C>A C[C>A]T -5.331806e-02
#> 221  C>G C[C>G]A  1.140520e-02
#> 222  C>G C[C>G]C  7.140612e-03
#> 223  C>G C[C>G]G -1.907500e-03
#> 224  C>G C[C>G]T  1.539370e-02
#> 225  C>T C[C>T]A  3.378900e-03
#> 226  C>T C[C>T]C  8.591852e-03
#> 227  C>T C[C>T]G -1.074480e-03
#> 228  C>T C[C>T]T  6.958618e-03
#> 229  T>A C[T>A]A -6.461466e-03
#> 230  T>A C[T>A]C  5.613680e-03
#> 231  T>A C[T>A]G -1.930279e-02
#> 232  T>A C[T>A]T  9.035457e-04
#> 233  T>C C[T>C]A  3.648386e-03
#> 234  T>C C[T>C]C  1.217409e-02
#> 235  T>C C[T>C]G  4.598236e-03
#> 236  T>C C[T>C]T  1.332373e-02
#> 237  T>G C[T>G]A  3.499204e-03
#> 238  T>G C[T>G]C  6.617684e-03
#> 239  T>G C[T>G]G  6.908170e-03
#> 240  T>G C[T>G]T  5.866854e-03
#> 241  C>A G[C>A]A -2.452313e-02
#> 242  C>A G[C>A]C -2.362341e-02
#> 243  C>A G[C>A]G -1.400985e-02
#> 244  C>A G[C>A]T -1.116893e-02
#> 245  C>G G[C>G]A  7.423360e-03
#> 246  C>G G[C>G]C  2.844868e-03
#> 247  C>G G[C>G]G -8.893042e-04
#> 248  C>G G[C>G]T  1.162307e-02
#> 249  C>T G[C>T]A  1.233306e-02
#> 250  C>T G[C>T]C  8.833372e-03
#> 251  C>T G[C>T]G  4.636514e-05
#> 252  C>T G[C>T]T  7.465936e-03
#> 253  T>A G[T>A]A -2.148785e-03
#> 254  T>A G[T>A]C  3.864925e-03
#> 255  T>A G[T>A]G -1.509206e-03
#> 256  T>A G[T>A]T  6.876089e-03
#> 257  T>C G[T>C]A  7.037202e-03
#> 258  T>C G[T>C]C  4.443212e-03
#> 259  T>C G[T>C]G  8.568192e-03
#> 260  T>C G[T>C]T  9.580148e-03
#> 261  T>G G[T>G]A  3.799073e-03
#> 262  T>G G[T>G]C  2.079199e-03
#> 263  T>G G[T>G]G  8.507909e-03
#> 264  T>G G[T>G]T  5.657852e-03
#> 265  C>A T[C>A]A -3.158094e-02
#> 266  C>A T[C>A]C -3.022953e-02
#> 267  C>A T[C>A]G -9.694986e-03
#> 268  C>A T[C>A]T -3.165580e-02
#> 269  C>G T[C>G]A  9.776192e-03
#> 270  C>G T[C>G]C  1.440654e-02
#> 271  C>G T[C>G]G -1.338357e-03
#> 272  C>G T[C>G]T  2.032369e-02
#> 273  C>T T[C>T]A  2.031125e-03
#> 274  C>T T[C>T]C  6.334056e-03
#> 275  C>T T[C>T]G  7.230880e-04
#> 276  C>T T[C>T]T  9.583094e-04
#> 277  T>A T[T>A]A -2.342574e-03
#> 278  T>A T[T>A]C  6.308687e-03
#> 279  T>A T[T>A]G -2.853877e-03
#> 280  T>A T[T>A]T  7.409249e-03
#> 281  T>C T[T>C]A  1.263288e-02
#> 282  T>C T[T>C]C  8.650741e-03
#> 283  T>C T[T>C]G  5.230234e-03
#> 284  T>C T[T>C]T  1.331831e-02
#> 285  T>G T[T>G]A  6.788732e-03
#> 286  T>G T[T>G]C  6.223839e-03
#> 287  T>G T[T>G]G  5.774346e-03
#> 288  T>G T[T>G]T  9.972652e-03