Showing content from https://cran.r-project.org/web/packages/rms/../effectsize/vignettes/statistical_power.R below:
## ----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" )
RetroSearch is an open source project built by @garambo
| Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4