back_project: back_project

In ropensci/DoOR.functions: Integrating Heterogeneous Odorant Response Data into a Common Response Model: A DoOR to the Complete Olfactome

back_project

Description

project the model response values back into your scale of interest (spike, deltaF/F...)

Usage

back_project(template.data, responding.unit,
  response_matrix = door_default_values("door_response_matrix"))

Arguments

template.data	data frame, the data that provides the scale to transform to, containing InChIKeys in a column called "odorants" and responses in a column called "responses"
responding.unit	character, the name of the receptor/OSN/glomerulus which responses should be transformed
response_matrix	DoOR response matrix, the source data is picked from here

Details

The process of back projection is the following:

• 1. rescale both data sets to [0,1],

• 2. find the best fitting model between "bp.data" and "cons.data" (lowest MD value),

• 3. project the consensus data onto the fitted curv, these are now our normalized, back projected responses

• 4. rescale all responses to the scale of the original data via a linear fit.

Value

Output of back_project is a list containing a data frame with the back_projected data, the original data, the data used as a template and the data that resulted from fitting source and template (before rescaling to the template scale). additionaly the parameters of the linear fit between the source and template response scale is returned.

Author(s)

Daniel Münch <[e-mail]>

Shouwen Ma <[e-mail]>

Examples

# load some data sets
data(Or22a)
data(door_response_matrix)

# create example data we are going to back project
template.data <- data.frame(odorants = Or22a$InChIKey,
                            responses = Or22a$Hallem.2004.EN)

# run back_project and plot the results
bp <- back_project(template.data, "Or22a")

plot(bp$back_projected$original.data,
     bp$back_projected$back_projected.data,
     xlab = "DoOR consensus response",
     ylab = "back_projected data [spikes, Hallem.2004.EN]"
)

ropensci/DoOR.functions documentation built on Feb. 22, 2024, 9:44 a.m.