acorg / Racmacs

Racmacs R package for performing antigenic cartography
https://acorg.github.io/Racmacs/
GNU Affero General Public License v3.0
20 stars 9 forks source link

Trying to remove all antigens or sera from a map causes a crash #116

Open drserajames opened 2 years ago

drserajames commented 2 years ago

If you use removeAntigens and try to remove all the antigens from the map, this causes RStudio to crash. I ran the code in the terminal - I don't know if this error is more informative? See below.

I don't know why anyone would want to remove all antigens/sera, so an informative error message would be good e.g. Error: not all antigens/sera can be removed from map

antigens

> library(Racmacs)
> packageVersion("Racmacs")
[1] ‘1.1.35’
> set.seed(850909)
> dat <- matrix(10*2^round(10*runif(100)), ncol=10)
> map <- make.acmap(dat)
Performing 100 optimizations
================================================================================
Optimization runs complete
Took 0.09 secs

> 
> removeAntigens(map, 1:10)
R(70834,0x10edd1e00) malloc: *** error for object 0x7ffee4918db0: pointer being freed was not allocated
R(70834,0x10edd1e00) malloc: *** set a breakpoint in malloc_error_break to debug
Abort trap: 6

sera

> packageVersion("Racmacs")
[1] ‘1.1.35’
> set.seed(850909)
> dat <- matrix(10*2^round(10*runif(100)), ncol=10)
> map <- make.acmap(dat)
Performing 100 optimizations
================================================================================
Optimization runs complete
Took 0.17 secs

> 
> removeSera(map, 1:10)
R(67423,0x10595ce00) malloc: *** error for object 0x7ffeeee9fdb0: pointer being freed was not allocated
R(67423,0x10595ce00) malloc: *** set a breakpoint in malloc_error_break to debug
Abort trap: 6