This documents a sample view of `funkycells`

. Every
function will be examined throughout this vignette, feel free to see the
process from beginning to end or jumping to a given step in the
document. Although each step in the `funkycells`

process is
addressed, detailed information on each function is left to their
respective documentation.

This package is interested in investigating spatial data. Although potentially combined, in the most raw form the data of interest would be (1) agent, or cell, information (locations, patient, image, etc.) and (2) meta-data (patient demographics, etc.).

Rather than loading in data, we will generate data with such
properties. The cell data is constructed using
`simulatePP()`

, as given below.

```
set.seed(123)
cell_data <- simulatePP(
agentVarData =
data.frame(
"outcome" = c(0, 1),
"A" = c(0, 0),
"B" = c(1 / 100, 1 / 500),
"C" = c(1 / 500, 1 / 250),
"D" = c(1 / 100, 1 / 100)
),
agentKappaData =
data.frame(
"agent" = c("A", "B", "C", "D"),
"clusterAgent" = c(NA, "A", NA, NA),
"kappa" = c(20, 5, 15, 15)
),
unitsPerOutcome = 40,
replicatesPerUnit = 1,
silent = FALSE
)
#> Outcome: 0 (1/2)
#> Outcome: 1 (2/2)
```

We look at some pictures using `plotPP()`

.

```
plotPP(cell_data[cell_data$replicate == 1, c("x", "y", "type")],
dropAxes = TRUE,
colorGuide = "none",
xlim = c(0, 1),
ylim = c(0, 1)
)
```

In funkycells we often explain the cell interactions using \(K\) functions, although many functions are feasible. For example, K functions could be computed as the following function.

```
k_data <- getKFunction(
data = cell_data[cell_data$unit == "u1", -1],
agents = c("A", "A"),
unit = "unit",
replicate = "replicate",
rCheckVals = seq(0, 0.25, 0.01),
xRange = c(0, 1),
yRange = c(0, 1)
)
```

However, `funkyModel()`

does not use functions directly,
but the principle components of functions. Because of the commonness of
\(K\) functions, currently we
implemented the steps in the same function. Although in our data, the
types are in a single column, the functions accepts when the cell
properties are defined in many columns (see documentation on
`getKsPCAData()`

). We also note that we only consider
one-direction of the interactions to reduce variables and improve model
power.

```
cells <- unique(cell_data$type)
cells_interactions <- rbind(
data.frame(t(combn(cells, 2))),
data.frame("X1" = cells, "X2" = cells)
)
pca_data <- getKsPCAData(
data = cell_data,
outcome = "outcome",
unit = "unit",
replicate = "replicate",
rCheckVals = seq(0, 0.25, 0.01),
agents_df = cells_interactions,
xRange = c(0, 1), yRange = c(0, 1),
nPCs = 3
)
#> PCA Pairs (10): 1, 2, 3, 4, 5, 6, 7, 8, 9, 10
```

Upon this data, we attach the meta information. For example, to this data we will add gender (with no effect), and age (with some effect).

We now jump into the model. We use `funkyModel()`

as it
employs cross-validation and permutation to create accurate models.

```
set.seed(123)
model_fc <- funkyModel(
data = pcaMeta,
outcome = "outcome",
unit = "unit",
metaNames = c("gender", "age")
)
```

The variable importance plot is returned (plot not shown as this was not run for computational reasons).

The results of the model from `funkyModel()`

also includes
the model from `funkyForest()`

on the original data, which is
necessary when trying to predict new data. Suppose we create a few new
people (i.e. get more data) and wish to predict. This would be done as
follows.

```
set.seed(12345)
cell_data_pred <- simulatePP(
agentVarData =
data.frame(
"outcome" = c(0, 1),
"A" = c(0, 0),
"B" = c(1 / 100, 1 / 500),
"C" = c(1 / 500, 1 / 250),
"D" = c(1 / 100, 1 / 100)
),
agentKappaData =
data.frame(
"agent" = c("A", "B", "C", "D"),
"clusterAgent" = c(NA, "A", NA, NA),
"kappa" = c(20, 5, 15, 15)
),
unitsPerOutcome = 2,
replicatesPerUnit = 2,
silent = FALSE
)
pca_data_pred <- getKsPCAData(
data = cell_data_pred,
outcome = "outcome",
unit = "unit",
replicate = "replicate",
rCheckVals = seq(0, 0.25, 0.01),
agents_df = cells_interactions,
xRange = c(0, 1), yRange = c(0, 1),
nPCs = 3
)
set.seed(12345)
pcaMeta_pred <- simulateMeta(pca_data_pred,
outcome = "outcome",
metaInfo = data.frame(
"var" = c("gender", "age"),
"rdist" = c("rbinom", "rnorm"),
"outcome_0" = c("0.5", "30"),
"outcome_1" = c("0.5", "31")
)
)
predictions <- predict_funkyForest(model = model_fc$model, data_pred = pcaMeta_pred[-1])
```

We also can compare K functions by group.

```
ab_outcome0 <- getKFunction(cell_data[cell_data$outcome == 0, -1],
agents = c("A", "B"), unit = "unit",
replicate = "replicate",
rCheckVals = seq(0, 0.25, 0.01)
)
ab_outcome1 <- getKFunction(cell_data[cell_data$outcome == 1, -1],
agents = c("A", "B"), unit = "unit",
replicate = "replicate",
rCheckVals = seq(0, 0.25, 0.01)
)
ab_outcome0_long <- tidyr::pivot_longer(data = ab_outcome0, cols = K1:K15)
ab_outcome1_long <- tidyr::pivot_longer(data = ab_outcome1, cols = K1:K15)
data_k_plot <- rbind(
data.frame(
"r" = ab_outcome0_long$r,
"K" = ab_outcome0_long$value,
"unit" = ab_outcome0_long$name,
"outcome" = "0"
),
data.frame(
"r" = ab_outcome1_long$r,
"K" = ab_outcome1_long$value,
"unit" = paste0(ab_outcome1_long$name, "_1"),
"outcome" = "1"
)
)
plot_K_functions(data_k_plot)
```

We can also consider ROC curves (Figure may not be shown due to computational reasons).

There are many uses for funkycells, but we hope this gives a taste of the current functionality and aids in understanding of the documentation provided.