---
title: "_ChemmineOB_: Interface to a Subset of OpenBabel Functionalities"
author: "Authors: Kevin Horan, [Thomas Girke](mailto:thomas.girke@ucr.edu)"
date: "Last update: `r format(Sys.time(), '%d %B, %Y')`"
package: "`r pkg_ver('ChemmineOB')`"
output:
BiocStyle::html_document:
toc: true
toc_depth: 3
fig_caption: yes
fontsize: 14pt
bibliography: references.bib
---
```{r style, echo = FALSE, results = 'asis'}
BiocStyle::markdown()
options(width=100, max.print=1000)
knitr::opts_chunk$set(
eval=as.logical(Sys.getenv("KNITR_EVAL", "TRUE")),
cache=as.logical(Sys.getenv("KNITR_CACHE", "TRUE")))
```
```{r setup, echo=FALSE, messages=FALSE, warnings=FALSE}
suppressPackageStartupMessages({
library(ChemmineOB)
})
```
Note: the most recent version of this tutorial can be found here and a short overview slide show [here](http://faculty.ucr.edu/~tgirke/HTML_Presentations/Manuals/Workshop_Dec_5_8_2014/Rcheminfo/Cheminfo.pdf).
# Introduction
`ChemmineOB` provides an R interface to a subset of
cheminformatics functionalities implemented by the OpelBabel C++ project
[@18328109; @21982300]. OpenBabel is an open source
cheminformatics toolbox that includes utilities for structure format
interconversions, descriptor calculations, compound similarity searching
and more. `ChemineOB` aims to make a subset of these
utilities available from within R. For non-developers,
`ChemineOB` is primarily intended to be used from
`ChemmineR` [@Cao2008c; @Backman2011a; @Wang2013a] as an add-on package
rather than used directly.
# Installation
To use the `ChemmineOB` package on Linux or Mac, OpenBabel
2.3.0 or greater needs to be installed on a system. On Linux systems,
the OpenBabel header files are also required in order to compile `ChemmineOB`. The windows distribution
will include its own version of OpenBabel. The OpenBabel site
() provides excellent
instructions for installing the OpenBabel software on Mac or Linux
systems. The `ChemmineR` and `ChemmineOB`
packages can be installed from within R with:
```{r eval=FALSE, tidy=FALSE}
if (!requireNamespace("BiocManager", quietly=TRUE))
install.packages("BiocManager")
BiocManager::install(c("ChemmineR", "ChemmineOB"))
library("ChemmineR")
library("ChemmineOB")
```
If the installation fails on Linux, you may need to manually set the
locations of the open babel libraries and header files. This is best
done through configure flags. For example, at the command prompt do:
```Bash
$ R CMD INSTALL --configure-args='--with-openbabel-include=... --with-openbabel-lib=...'
```
where the '...' are replaced by the relevant paths. See the README file
for more details.
# User Manual in ChemmineR Vignette
Detailed instructions for using `ChemmineOB` are provided
in the vignette of the `ChemmineR` package instead of this
document. The main reason for consolidating the documentation in one
central document rather than distributing it across several vignettes is
that it helps minimizing duplications and inconsistencies. It also is
the more suitable format for providing a task-oriented description of
functionalities for users. To obtain an overview of the OpenBabel
utilities supported by `ChemmineOB`, we recommend
consulting the *OpenBabel Functions* section of the
`ChemmineR` vignette. To open the `ChemmineR`
vignette from R, one can use the following command.
```{r eval=FALSE, tidy=FALSE}
vignette("ChemmineR")
```
# SWIG Interface (For R developers)
`ChemmineOB` now includes wrapper functions for all of OpenBabel, as genereted by
[SWIG](http://www.swig.org). We still maintain our own set of functions to provide
better integration with R in general and `ChemmineR` specifically.
If you are familiar with the [Open Babel API](http://openbabel.org/api/2.3/), using the SWIG wrapper should be
similar, once you know a few conventions used. You can look at the R code in this
package to see examples of these.
- New objects are created with a function with the same name as the object you
want to construct. So instead of
```
OBConversion *x = new OBConversion(...)
```
in R you would have:
```
x = OBConversion(...)
```
- Methods on objects are called through a function whose name is the name of the
object type concatenated with the method name, separated by an underscore. Then the
first argument of this method is the instance of the object. For example,
instead of:
```
x->AddOption(...)
```
we have:
```
OBConversion_AddOption(x,...)
```
- If you need a char* you can create one in R with the `stringp` function.
The char* pointer can be accessed with the `cast` slot. The value can
be retrieved from the `value` slot. For example:
```
result = stringp()
OBDescriptor_GetStringValue(... , result$cast())
stringValue = result$value()
```
- STL vectors of various types are available via functions named as "vector{TYPE NAME}". For example,
a vector of unsigned int is created wit the function "vectorUnsignedInt". You can get the size of
this vector with "vectorUnsignedInt_size(vectorVariable)".
There are still many special cases however. The [SWIG documentation](http://www.swig.org/Doc2.0/SWIGDocumentation.html)
can help, as well as browsing the generated R code in R/ChemmineOB.R.
# Version Information
```{r sessionInfo, results='asis'}
sessionInfo()
```
# Funding
This software was developed with funding from the National Science
Foundation:
[ABI-0957099](http://www.nsf.gov/awardsearch/showAward.do?AwardNumber=0957099),
2010-0520325 and IGERT-0504249.
# References