Nothing
library("compositions")
# ---- test name_invTrans -----------------------------------------------------
expect_equal(name_invTrans("ilr( CVAR )", "ilr"), "CVAR")
expect_equal(name_invTrans("ilr( CVAR, V)", "ilr"), "CVAR")
expect_equal(name_invTrans("ilr(V = ilrBase(D = 3), CVAR)", "ilr"), "CVAR")
expect_equal(name_invTrans("ilr(cbind(left, right, other))", "ilr"), "cbind(left, right, other)")
# ---- test whichTrans --------------------------------------------------------
expect_equal({
tt <- rice_yields[1:10,]
tt$"ilr(TEMPERATURES)" <- ilr(tt$TEMPERATURES)
whichTrans(tt["ilr(TEMPERATURES)"])[["name"]]
}, "ilr")
expect_equal({
tt <- rice_yields[1:10,]
tt$"ilr(TEMPERATURES)" <- ilr(tt$TEMPERATURES)
whichTrans(tt["ilr(TEMPERATURES)"])[["base_K"]]
}, {
V <- ilrBase(D = 3)
rownames(V) <- colnames(tt$TEMPERATURES)
colnames(V) <- paste0("ilr(TEMPERATURES)", 1:2)
V
})
expect_equal({
tt <- rice_yields[1:10,]
tt$"alr(TEMPERATURES)" <- alr(tt$TEMPERATURES)
whichTrans(tt["alr(TEMPERATURES)"])
}, {
F_alr <- cbind(diag(2),-1L)
dimnames(F_alr) <- list(c("alr(TEMPERATURES)1", "alr(TEMPERATURES)2"),c("LOW","MEDIUM","HIGH"))
K_alr <- clrBase(3)[,-3]
dimnames(K_alr) <- dimnames(t(F_alr))
list("name" = "alr", "base_F" = F_alr, "base_K" = K_alr)
})
expect_error({
whichTrans(rice_yields[1:10,]$TEMPERATURES)
})
# ---- alr_K2F ----------------------------------------------------------------
expect_equal(
clrBase(3)[,-3] |> alr_K2F(),
cbind(diag(2),-1),
info = "F from K with ref = 3")
expect_equal(
clrBase(3)[,-2] |> alr_K2F(),
cbind(c(1,0),-1,c(0,1)),
info = "F from K with ref = 2")
# ---- transformationSummary 1. (YX compo) ------------------------------------
expect_equal({
lr1 <- function(x) ilr(x, ilrBase(D = 3))
lr2 <- function(x) ilr(x, ilrBase(D = 2))
lr3 <- function(x) ilr(x, ilrBase(D = 3)[3:1,])
tt <- election[1:20,]
res <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~
lr2(cbind(Age_1839, Age_4064)) +
lr3(cbind(Educ_BeforeHighschool, Educ_Highschool, Educ_Higher)) +
log(unemp_rate),
data = tt))
res$COEF_CLR
},
{
V <- ilrBase(D = 3)
lr1 <- function(x) ilr(x, ilrBase(D = 3)[3:1,])
lr2 <- function(x) ilr(x, ilrBase(D = 2)[2:1,])
lr3 <- function(x) ilr(x, ilrBase(D = 3)[c(2,1,3),])
tt <- election[1:20,]
res <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~
lr2(cbind(Age_1839, Age_4064)) +
lr3(cbind(Educ_BeforeHighschool, Educ_Highschool, Educ_Higher)) +
log(unemp_rate), data = tt))
res$COEF_CLR
},info = "Three diffrent ilr bases in X and Y")
expect_equal({
lr1 <- function(x) alr(x,ivar = 1)
lr2 <- function(x) alr(x,ivar = 2)
lr3 <- function(x) alr(x,ivar = 3)
tt <- election[1:20,]
res1 <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~
lr2(cbind(Age_1839, Age_4064)) +
lr3(cbind(Educ_BeforeHighschool, Educ_Highschool, Educ_Higher)),
data = tt))
res1[,c("COEF_SIMPLEX","COEF_CLR")]
},
{
lr1 <- function(x) alr(x,ivar = 3)
lr2 <- function(x) alr(x,ivar = 1)
lr3 <- function(x) alr(x,ivar = 2)
tt <- election[1:20,]
res2 <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~
lr2(cbind(Age_1839, Age_4064)) +
lr3(cbind(Educ_BeforeHighschool, Educ_Highschool, Educ_Higher))
, data = tt))
res2[,c("COEF_SIMPLEX","COEF_CLR")]
},info = "Three diffrent alr bases in X and Y.")
expect_equal({
lr1 <- function(x) alr(x,ivar = 1)
lr2 <- function(x) alr(x,ivar = 2)
lr3 <- function(x) alr(x,ivar = 3)
tt <- election[1:20,]
res1 <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~
lr2(cbind(Age_1839, Age_4064, Age_65plus)) +
lr3(cbind(Educ_BeforeHighschool, Educ_Highschool, Educ_Higher)) +
unemp_rate, data = tt))
res1[,c("COEF_SIMPLEX","COEF_CLR")]
},
{
lr1 <- function(x) alr(x,ivar = 3)
lr2 <- function(x) alr(x,ivar = 1)
lr3 <- function(x) alr(x,ivar = 2)
tt <- election[1:20,]
res2 <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~
lr2(cbind(Age_1839, Age_4064, Age_65plus)) +
lr3(cbind(Educ_BeforeHighschool, Educ_Highschool, Educ_Higher)) +
unemp_rate, data = tt))
res2[,c("COEF_SIMPLEX","COEF_CLR")]
},info = "Two diffrent alr bases in X and a diffrent ilr Y.")
expect_equal({
lr1 <- function(x) alr(x,ivar = 1)
lr2 <- function(x) alr(x,ivar = 2)
lr3 <- function(x) alr(x,ivar = 3)
tt <- election[1:20,]
res1 <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~
lr2(cbind(Age_1839, Age_4064, Age_65plus)) +
lr3(cbind(Educ_BeforeHighschool, Educ_Highschool, Educ_Higher)) +
unemp_rate, data = tt))
res1[,c("COEF_SIMPLEX","COEF_CLR")]
},
{
lr1 <- function(x) alr(x,ivar = 3)
lr2 <- function(x) alr(x,ivar = 1)
lr3 <- function(x) alr(x,ivar = 2)
tt <- election[1:20,]
res2 <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~
lr2(cbind(Age_1839, Age_4064, Age_65plus)) +
lr3(cbind(Educ_BeforeHighschool, Educ_Highschool, Educ_Higher)) +
unemp_rate, data = tt))
res2[,c("COEF_SIMPLEX","COEF_CLR")]
},info = "Two diffrent ilr bases in X and a diffrent alr Y.")
# ---- transformationSummary 2. (Y compo) -------------------------------------
expect_equal({
lr1 <- function(x) ilr(x)
tt <- election[1:20,]
res1 <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~ unemp_rate, data = tt))
res1[,c("COEF_SIMPLEX","COEF_CLR")]
},
{
lr1 <- function(x) alr(x,ivar = 3)
tt <- election[1:20,]
res2 <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~ unemp_rate, data = tt))
res2[,c("COEF_SIMPLEX","COEF_CLR")]
},info = "Diffrent ilr and alr bases in Y.")
expect_equal({
lr1 <- function(x) alr(x,ivar = 1)
tt <- election[1:20,]
res1 <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~ unemp_rate, data = tt))
res1[,c("COEF_SIMPLEX","COEF_CLR")]
},
{
lr1 <- function(x) alr(x,ivar = 3)
tt <- election[1:20,]
res2 <- transformationSummary(lm(
lr1(cbind(left, right, extreme_right)) ~ unemp_rate, data = tt))
res2[,c("COEF_SIMPLEX","COEF_CLR")]
},info = "Diffrent alr bases in Y.")
# ---- transformationSummary 3. (X compo) -------------------------------------
expect_equal({
V <- ilrBase(D = 3)
lr <- function(x) ilr(x, V)
tt <- rice_yields[1:20,]
res <- transformationSummary(lm(YIELD ~ PRECIPITATION + lr(TEMPERATURES), data = tt))
# res["lr(TEMPERATURES)","COEF_COORD"][[1]] # differ
res["lr(TEMPERATURES)",c("COEF_SIMPLEX","COEF_CLR")]
},
{
V <- ilrBase(D = 3)[3:1,]
lr <- function(x) ilr(x, V)
tt <- rice_yields[1:20,]
res <- transformationSummary(lm(YIELD ~ PRECIPITATION + lr(TEMPERATURES), data = tt))
# res["lr(TEMPERATURES)","COEF_COORD"][[1]] # differ
res["lr(TEMPERATURES)",c("COEF_SIMPLEX","COEF_CLR")]
},info = "Two diffrent ilr bases.")
expect_equal({
lr <- function(x) alr(x, ivar = 1)
tt <- rice_yields[1:20,]
res <- transformationSummary(lm(YIELD ~ PRECIPITATION + lr(TEMPERATURES), data = tt))
res["lr(TEMPERATURES)",c("COEF_SIMPLEX","COEF_CLR")]
},
{
lr <- function(x) alr(x, ivar = 3)
tt <- rice_yields[1:20,]
res <- transformationSummary(lm(YIELD ~ PRECIPITATION + lr(TEMPERATURES), data = tt))
res["lr(TEMPERATURES)",c("COEF_SIMPLEX","COEF_CLR")]
},info = "Two diffrent alr bases.")
expect_equal({
lr <- compositions::ilr
tt <- rice_yields[1:20,]
res <- transformationSummary(lm(YIELD ~ PRECIPITATION + lr(TEMPERATURES), data = tt))
res["lr(TEMPERATURES)",c("COEF_CLR")]
},
{
lr <- compositions::alr
tt <- rice_yields[1:20,]
res <- transformationSummary(lm(YIELD ~ PRECIPITATION + lr(TEMPERATURES), data = tt))
res["lr(TEMPERATURES)",c("COEF_CLR")]
}, info = "Ilr and alr bases.")
expect_equivalent({
tt <- rice_yields[1:20,]
res <- transformationSummary(lm(YIELD ~ PRECIPITATION + ilr(TEMPERATURES), data = tt, weights = 1:20))
res["YIELD","VARCOV_CLR"]
},
{
tt <- rice_yields[1:20,]
res <- transformationSummary(lm(YIELD ~ PRECIPITATION + ilr(TEMPERATURES), data = tt))
res["YIELD","VARCOV_CLR"]
}, info = "Weighted and non-weighted regression.", tolerance = 10e10)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.