-
Notifications
You must be signed in to change notification settings - Fork 0
/
transform_inverse.R
93 lines (82 loc) · 2.28 KB
/
transform_inverse.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
# DATA TRANSFORMATION: INVERSE
#' @include AllGenerics.R
NULL
# Backtransform ================================================================
## CLR -------------------------------------------------------------------------
#' @export
#' @rdname transform_inverse
#' @aliases transform_inverse,CLR,missing-method
setMethod(
f = "transform_inverse",
signature = c(object = "CLR", origin = "missing"),
definition = function(object) {
y <- methods::as(object, "matrix") # Drop slots
y <- exp(y)
y <- y / rowSums(y)
dimnames(y) <- list(rownames(object), object@parts)
.CompositionMatrix(
y,
totals = object@totals,
samples = object@samples,
groups = object@groups
)
}
)
## ALR -------------------------------------------------------------------------
#' @export
#' @rdname transform_inverse
#' @aliases transform_inverse,ALR,missing-method
setMethod(
f = "transform_inverse",
signature = c(object = "ALR", origin = "missing"),
definition = function(object) {
y <- exp(object)
y <- y / (1 + rowSums(y))
z <- 1 - rowSums(y)
y <- cbind(y, z)
dimnames(y) <- list(rownames(object), object@parts)
y <- y[, object@order]
.CompositionMatrix(
y,
totals = object@totals,
samples = object@samples,
groups = object@groups
)
}
)
## ILR -------------------------------------------------------------------------
#' @export
#' @rdname transform_inverse
#' @aliases transform_inverse,ILR,missing-method
setMethod(
f = "transform_inverse",
signature = c(object = "ILR", origin = "missing"),
definition = function(object) {
y <- tcrossprod(object, object@base)
y <- exp(y)
y <- y / rowSums(y)
dimnames(y) <- list(rownames(object), object@parts)
y <- y[, object@order]
.CompositionMatrix(
y,
totals = object@totals,
samples = object@samples,
groups = object@groups
)
}
)
#' @export
#' @rdname transform_inverse
#' @aliases transform_inverse,matrix,ILR-method
setMethod(
f = "transform_inverse",
signature = c(object = "matrix", origin = "ILR"),
definition = function(object, origin) {
y <- tcrossprod(object, origin@base)
y <- exp(y)
y <- y / rowSums(y)
dimnames(y) <- list(rownames(object), origin@parts)
y <- y[, origin@order]
y
}
)