Tim Read plasmid02
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| # data from http://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/population-distribution-demography/geostat | |
| # Originally seen at http://spatial.ly/2014/08/population-lines/ | |
| # So, this blew up on both Reddit and Twitter. Two bugs fixed (southern Spain was a mess, | |
| # and some countries where missing -- measure twice, submit once, damnit), and two silly superflous lines removed after | |
| # @hadleywickham pointed that out. Also, switched from geom_segment to geom_line. | |
| # The result of the code below can be seen at http://imgur.com/ob8c8ph | |
| library(tidyverse) |
JimGrange
/ effect_sizes.R
Last active
November 25, 2021 15:48
Low power = inaccurate effect size estimates
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #------------------------------------------------------------------------------ | |
| rm(list = ls()) | |
| set.seed(50) | |
| # function for generating random draws from multivariate distribution | |
| # n = number of draws; p = number of variables | |
| # u = mean of each variable; s = SD of each variable | |
| # corMat = correlation matrix | |
| mvrnorm <- function(n, p, u, s, corMat) { |