context("NN_keras")
df <- na.omit(airquality)
responseVars <- "Ozone"
predInput <- df
varScale<- seq(-100, 100, length.out=ncol(df))
names(varScale)<- names(df)
dfCat <- iris[iris$Species %in% c("setosa", "versicolor"), ] # Only 2 categories supported
dfCat$Species <- as.character(dfCat$Species)
responseVarsCat <- "Species"
predInputCat <- dfCat
crossValStrategy<- c("Kfold", "bootstrap")
crossValRatio<- c(train=0.6, test=0.2, validate=0.2)
k<- 2
idVars<- character()
epochs<- 1
maskNA<- -999
replicates<- 2
repVi<- 2
summarizePred<- TRUE
shap<- TRUE
hidden_shape<- 2
batch_size<- "all"
scaleDataset<- FALSE
tempdirRaster<- tempdir()
dir.create(tempdirRaster, showWarnings=FALSE)
filenameRasterPred<- paste0(tempdirRaster, "/testMap.grd")
baseFilenameRasterPred<- paste0(tempdirRaster, "/testMap")
baseFilenameNN<- paste0(tempdir(), "/testNN")
nCoresRaster<- 2
variableResponse<- TRUE
DALEXexplainer<- TRUE
save_validateset<- TRUE
NNmodel<- TRUE
caseClass<- c(rep("A", 23), rep("B", 75), rep("C", 13)) ## TODO: use it on tests!
caseClassCat<- dfCat$Species ## TODO: use it on tests!
weight<- "class"
verbose<- 2
verbose<- 0
test_that("pipe_keras works", {
result<- list()
# future::plan(future::sequential, split=TRUE)
future::plan(future::multisession)
system.time(result$resp1summarizedPred<- pipe_keras(df=df, predInput=predInput, responseVars=responseVars,
epochs=epochs, repVi=repVi,
crossValStrategy=crossValStrategy[1], k=k, replicates=replicates,
batch_size=batch_size, hidden_shape=c(2, 2),
baseFilenameNN=paste0(baseFilenameNN, "-resp1summarizedPred"), DALEXexplainer=DALEXexplainer, variableResponse=variableResponse, save_validateset=save_validateset,
crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose))
system.time(result$resp2summarizedPred<- pipe_keras(df=df, predInput=predInput, responseVars=1:2,
epochs=epochs, maskNA=maskNA, repVi=repVi,
crossValStrategy=crossValStrategy[2], replicates=replicates,
batch_size=batch_size, hidden_shape=hidden_shape,
baseFilenameNN=paste0(baseFilenameNN, "-resp2summarizedPred"), DALEXexplainer=DALEXexplainer, variableResponse=variableResponse, save_validateset=save_validateset,
crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose))
system.time(result$resp1Cat<- pipe_keras(df=dfCat, predInput=rev(predInputCat), responseVars=responseVarsCat,
epochs=epochs, maskNA=maskNA, repVi=10, # check names with 2 digits
crossValStrategy=crossValStrategy[2], replicates=10, # check names with 2 digits
hidden_shape=hidden_shape, batch_size=batch_size, summarizePred=FALSE,
baseFilenameNN=paste0(baseFilenameNN, "-resp1"), DALEXexplainer=DALEXexplainer, variableResponse=variableResponse, save_validateset=save_validateset,
crossValRatio=crossValRatio[1], NNmodel=NNmodel, verbose=verbose))
system.time(result$resp2<- pipe_keras(df=df, predInput=rev(predInput), responseVars=1:2,
epochs=epochs, repVi=repVi,
crossValStrategy=crossValStrategy[1], k=k, replicates=replicates,
hidden_shape=hidden_shape, batch_size=batch_size, summarizePred=FALSE,
baseFilenameNN=paste0(baseFilenameNN, "-resp2"), DALEXexplainer=DALEXexplainer, variableResponse=variableResponse, save_validateset=save_validateset,
crossValRatio=c(train=0.8, test=0.2), NNmodel=NNmodel, verbose=verbose))
tmp<- lapply(result, function(x) expect_s3_class(x, class="pipe_result.keras"))
tmp<- lapply(result, function(x){
expect_s3_class(x$performance, class="data.frame")
reps<- nrow(x$performance)
if (x$params$crossValStrategy == "bootstrap"){
expect_equal(rownames(x$performance), expected=paste0("rep", formatC(1:reps, format="d", flag="0", width=nchar(reps))))
} else if (x$params$crossValStrategy == "Kfold") {
expect_equal(rownames(x$performance), expected=paste0("Fold", 2:k, ".Rep", rep(1:replicates, each=k-1))) # Fold2:k (Fold1 for validationset)
}
})
tmp<- lapply(result, function(x){
expect_type(x$scale, type="list")
expect_equal(unique(lapply(x$scale, names)), expected=list(c("mean", "sd")))
})
tmp<- expect_s3_class(result$resp1summarizedPred$shap, class = "shapviz")
tmp<- expect_s3_class(result$resp2summarizedPred$shap, class = "mshapviz")
tmp<- expect_s3_class(result$resp1$shap, class = "shapviz")
tmp<- expect_s3_class(result$resp2$shap, class = "mshapviz")
tmp<- lapply(result, function(x){
expect_type(x$vi, type="double")
reps<- nrow(x$performance)
repsVi<- x$params$repVi
expectedColnames<- paste0(rep(paste0("rep", formatC(1:reps, format="d", flag="0", width=nchar(reps))), each=repsVi), "_",
rep(paste0("perm", formatC(1:repsVi, format="d", flag="0", width=nchar(repsVi))), times=reps))
expect_equal(colnames(x$vi), expected=expectedColnames)
})
tmp<- lapply(result, function(x){
lapply(x$variableResponse, expect_s3_class, class="partial_dependence_explainer")
lapply(x$variableResponse, expect_s3_class, class="aggregated_profiles_explainer")
})
tmp<- lapply(result, function(x){
expect_type(x$variableCoef, type="list")
lapply(x$variableCoef, function(y){
if (ncol(y) > 4){
expectedColnames<- c("intercept", paste0("b", 1:(ncol(y) - 4)), "adj.r.squared", "r.squared", "degree")
} else {
expectedColnames<- c("intercept", "adj.r.squared", "r.squared", "degree")
}
expect_equal(colnames(y), expected=expectedColnames)
})
})
tmp<- lapply(result, function(x){
expect_type(x$predictions, type="list")
expect_type(x$predictions[[1]], type="double")
})
expectedColnames<- c(idVars, "Mean", "SD", "Naive SE", "2.5%", "25%", "50%", "75%", "97.5%")
tmp<- expect_equal(colnames(result$resp1summarizedPred$predictions[[1]]), expected=expectedColnames)
tmp<- expect_equal(unlist(unique(lapply(result$resp2summarizedPred$predictions, colnames))), expected=expectedColnames)
tmp<- expect_equal(colnames(result$resp1$predictions[[1]]), expected=c(idVars, paste0("rep", formatC(1:nrow(result$resp1$performance), format="d", flag="0", width=nchar(nrow(result$resp1$performance))))))
tmp<- expect_equal(unlist(unique(lapply(result$resp2$predictions, colnames))), expected=c(idVars, paste0("rep", formatC(1:nrow(result$resp2$performance), format="d", flag="0", width=nchar(nrow(result$resp2$performance))))))
tmp<- lapply(result, function(x){
expect_type(x$model, type="list")
lapply(x$model, function(y) expect_type(y, type="raw"))
lapply(x$model, function(y) expect_s3_class(keras::unserialize_model(y), class="keras.engine.training.Model"))
})
# dir(tempdir(), full.names=TRUE)
expect_true(any(grepl(baseFilenameNN, dir(tempdir(), full.names=TRUE))))
expect_equal(sum(grepl(baseFilenameNN, dir(tempdir(), full.names=TRUE))), sum(sapply(result, function(x) nrow(x$performance))))
tmp<- lapply(result, function(x){
expect_type(x$DALEXexplainer, type="list")
reps<- nrow(x$performance)
if (x$params$crossValStrategy == "bootstrap"){
expect_equal(names(x$DALEXexplainer), expected=paste0("rep", formatC(1:reps, format="d", flag="0", width=nchar(reps))))
} else if (x$params$crossValStrategy == "Kfold") {
expect_equal(names(x$DALEXexplainer), expected=paste0("Fold", 2:k, ".Rep", rep(1:replicates, each=k-1))) # Fold2:k (Fold1 for validationset)
}
lapply(x$DALEXexplainer, expect_s3_class, class="explainer")
})
tmp<- lapply(result, function(x){
expect_type(x$validateset, type="list")
reps<- nrow(x$performance)
if (x$params$crossValStrategy == "bootstrap"){
expect_equal(names(x$validateset), expected=paste0("rep", formatC(1:reps, format="d", flag="0", width=nchar(reps))))
} else if (x$params$crossValStrategy == "Kfold") {
expect_equal(names(x$validateset), expected=paste0("Fold", 2:k, ".Rep", rep(1:replicates, each=k-1))) # Fold2:k (Fold1 for validationset)
}
lapply(x$validateset, expect_s3_class, class="data.frame")
})
})
test_that("Predict with raster", {
predInputR<- raster::raster(nrows=3, ncols=5)
predInputR<- raster::stack(lapply(varScale, function(i){
raster::setValues(predInputR, runif(raster::ncell(predInputR)) * i)
}))
# Put some NAs to detect rotations
NAs<- expand.grid(col=1:ncol(predInputR), row=1:nrow(predInputR))
NAs<- NAs[NAs$row > NAs$col, ]
predInputR[NAs$row, NAs$col]<- NA
resultR<- list()
# predInput<- predInputR
## TODO: categorical variables for rasters
df<- df[, names(predInputR)]
suppressWarnings(future::plan(future::multicore))
filenameRasterPred<- paste0(tempdir(), "/testMap1.grd") # avoid overwrite
resultR$resp1summarizedPred<- pipe_keras(df, predInput=predInputR,
epochs=epochs, repVi=repVi,
crossValStrategy=crossValStrategy[1], k=k, replicates=replicates,
batch_size=batch_size, hidden_shape=hidden_shape, summarizePred=TRUE,
filenameRasterPred=filenameRasterPred, tempdirRaster=tempdirRaster, baseFilenameNN=baseFilenameNN,
DALEXexplainer=DALEXexplainer, crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose)
filenameRasterPred<- paste0(tempdir(), "/testMap2.grd") # avoid overwrite
resultR$resp1<- pipe_keras(df, predInput=predInputR[[rev(names(predInputR))]],
epochs=epochs, maskNA=maskNA, repVi=repVi,
crossValStrategy=crossValStrategy[2], replicates=replicates,
batch_size=batch_size, hidden_shape=hidden_shape, summarizePred=FALSE,
filenameRasterPred=filenameRasterPred, tempdirRaster=tempdirRaster, baseFilenameNN=baseFilenameNN,
DALEXexplainer=FALSE, crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose)
filenameRasterPred<- paste0(tempdir(), "/testMap3.grd") # avoid overwrite
resultR$resp2summarizedPred<- pipe_keras(df, predInput=predInputR, responseVars=1:2, epochs=epochs, maskNA=maskNA, repVi=repVi,
crossValStrategy=crossValStrategy[1], k=k, replicates=replicates, batch_size=batch_size, hidden_shape=hidden_shape,
summarizePred=TRUE, filenameRasterPred=filenameRasterPred, tempdirRaster=tempdirRaster, baseFilenameNN=baseFilenameNN,
DALEXexplainer=FALSE, crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose)
filenameRasterPred<- paste0(tempdir(), "/testMap4.grd") # avoid overwrite
resultR$resp2<- pipe_keras(df, predInput=predInputR, responseVars=1:2, epochs=epochs, repVi=repVi, crossValStrategy=crossValStrategy[2], replicates=replicates,
batch_size=batch_size, hidden_shape=hidden_shape,
summarizePred=FALSE, filenameRasterPred=filenameRasterPred, tempdirRaster=tempdirRaster, baseFilenameNN=baseFilenameNN,
DALEXexplainer=FALSE, crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose)
tmp<- lapply(resultR, function(x){
expect_s4_class(x$predictions, class="RasterBrick")
})
tmp<- expect_equal(names(resultR$resp1summarizedPred$predictions), expected=c("mean", "sd", "se"))
tmp<- expect_equal(names(resultR$resp1$predictions), expected=paste0("Ozone_rep", 1:replicates))
# lapply(resultR, function(x) names(x$predictions))
## Check NAs position
# plot(predInputR)
# plot(resultR$resp1summarizedPred$predictions)
# plot(resultR$resp1$predictions)
# plot(resultR$resp2summarizedPred$predictions)
# plot(resultR$resp2$predictions)
file.remove(dir(tempdir(), "testMap.+\\.gr(i|d)$", full.names=TRUE))
})
test_that("Future plans work", {
# options(future.globals.onReference = "error")
# Error in keras::reset_states(modelNN) : attempt to apply non-function
# Don't import/export python objects to/from code inside future for PSOCK and callR clusters
# https://cran.r-project.org/web/packages/future/vignettes/future-4-non-exportable-objects.html
future::plan(future::sequential, split=TRUE)
system.time(res<- pipe_keras(df=df, predInput=predInput, responseVars=responseVars, epochs=epochs, crossValStrategy=crossValStrategy[2], replicates=replicates, repVi=repVi, batch_size=batch_size,
hidden_shape=hidden_shape, DALEXexplainer=DALEXexplainer, crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose))
expect_s3_class(res, class="pipe_result.keras")
future::plan(future::multicore)
system.time(res<- pipe_keras(df=df, predInput=predInput, responseVars=responseVars, epochs=epochs, crossValStrategy=crossValStrategy[2], replicates=replicates, repVi=repVi, batch_size=batch_size,
hidden_shape=hidden_shape, DALEXexplainer=DALEXexplainer, crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose))
expect_s3_class(res, class="pipe_result.keras")
future::plan(future.callr::callr(workers=3))
system.time(res<- pipe_keras(df=df, predInput=predInput, responseVars=responseVars, epochs=epochs, crossValStrategy=crossValStrategy[2], replicates=replicates, repVi=repVi, batch_size=batch_size,
hidden_shape=hidden_shape, DALEXexplainer=DALEXexplainer, crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose))
expect_s3_class(res, class="pipe_result.keras")
future::plan(future::sequential)
system.time(res<- pipe_keras(df=df, predInput=predInput, responseVars=responseVars, epochs=epochs, crossValStrategy=crossValStrategy[2], replicates=replicates, repVi=repVi, batch_size=batch_size,
hidden_shape=hidden_shape, DALEXexplainer=DALEXexplainer, crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose))
expect_s3_class(res, class="pipe_result.keras")
})
test_that("scaleDataset", {
future::plan(future::multisession)
system.time(res<- pipe_keras(df=df, predInput=predInput, responseVars=responseVars, epochs=epochs, crossValStrategy=crossValStrategy[2], replicates=replicates, repVi=repVi,
batch_size=batch_size, scaleDataset=TRUE, hidden_shape=hidden_shape,
baseFilenameNN=baseFilenameNN, DALEXexplainer=DALEXexplainer, crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose))
expect_s3_class(res, class="pipe_result.keras")
predInputR<- raster::raster(nrows=5, ncols=5)
predInputR<- raster::stack(lapply(varScale, function(i){
raster::setValues(predInputR, runif(raster::ncell(predInputR)) * i)
}))
# predInput<- predInputR
## TODO: categorical variables for rasters
df<- df[, names(predInputR)]
filenameRasterPred<- paste0(tempdir(), "/testMapScaleDataset.grd") # avoid overwrite
res<- pipe_keras(df, predInput=predInputR, epochs=epochs, crossValStrategy=crossValStrategy[2], replicates=replicates, repVi=repVi, batch_size=batch_size,
scaleDataset=TRUE, hidden_shape=hidden_shape,
filenameRasterPred=filenameRasterPred, tempdirRaster=tempdirRaster, baseFilenameNN=baseFilenameNN,
DALEXexplainer=DALEXexplainer, crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose)
expect_s3_class(res, class="pipe_result.keras")
})
test_that("summary", {
future::plan(future::multisession)
system.time(res<- pipe_keras(df=df, predInput=predInput, responseVars=responseVars, epochs=epochs, crossValStrategy=crossValStrategy[2], replicates=replicates, repVi=repVi,
batch_size=batch_size, scaleDataset=TRUE, hidden_shape=hidden_shape,
baseFilenameNN=baseFilenameNN, DALEXexplainer=DALEXexplainer, crossValRatio=crossValRatio, NNmodel=NNmodel, verbose=verbose))
sres<- summary(res)
expect_s3_class(sres, class="summary.pipe_result.keras")
expect_type(sres, type="list")
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.