statistical_power.R
In effectsize: Indices of Effect Size

## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  warning = FALSE,
  message = FALSE,
  eval = requireNamespace("pwr", quietly = TRUE)
)

options(digits = 4L, knitr.kable.NA = "")

set.seed(123)

## -----------------------------------------------------------------------------
# library(pwr)
# library(effectsize)

## -----------------------------------------------------------------------------
# t <- t.test(mpg ~ am, data = mtcars)

## ----eval = FALSE-------------------------------------------------------------
# pwr.t2n.test(
#   n1 = ..., n2 = ...,
#   d = ...,
#   sig.level = ...,
#   power = ...,
#   alternative = ...
# )

## ----eval = FALSE-------------------------------------------------------------
# t_alt <- t.test(mtcars$mpg[mtcars$am == 0], mtcars$mpg[mtcars$am == 1])
# 
# effectsize(t_alt, type = "cohens_d")

## ----eval = FALSE-------------------------------------------------------------
# cohens_d(mpg ~ am, data = mtcars)

## -----------------------------------------------------------------------------
# t_to_d(
#   t = t$statistic,
#   df_error = t$parameter
# )

## -----------------------------------------------------------------------------
# (result <- cohens_d(mpg ~ am, data = mtcars))
# (Ns <- table(mtcars$am))
# 
# pwr.t2n.test(
#   n1 = Ns[1], n2 = Ns[2],
#   d = result[["Cohens_d"]],
#   sig.level = 0.05,
#   alternative = "two.sided"
# )