## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  fig.width = 7,
  fig.height = 6
)
set.seed(42)
library(balnet)

## -----------------------------------------------------------------------------
n <- 100
p <- 25
X <- matrix(rnorm(n * p), n, p)
W <- rbinom(n, 1, 1 / (1 + exp(1 - X[, 1])))

## -----------------------------------------------------------------------------
fit <- balnet(X, W)

## -----------------------------------------------------------------------------
print(fit)

## -----------------------------------------------------------------------------
smd.baseline <- (colMeans(X[W == 1, ]) - colMeans(X)) /
  (apply(X, 2, sd) * sqrt((n - 1) / n))
max(abs(smd.baseline))

## -----------------------------------------------------------------------------
plot(fit)

## -----------------------------------------------------------------------------
plot(fit, lambda = 0)

## -----------------------------------------------------------------------------
fit.att <- balnet(X, W, target = "ATT")
plot(fit.att, lambda = 0)