Skip to content

Instantly share code, notes, and snippets.

View bschneidr's full-sized avatar

Ben Schneider bschneidr

26 followers · 28 following
View GitHub Profile
@bschneidr
bschneidr / sample-variance-as-quadratic-form.R
Last active October 28, 2022 12:33
Expresses sample variance as a quadratic form
n = 5
# Establish quadratic form
quad_form_matrix <- matrix(nrow = n, ncol = n)
for (i in seq(n)) {
for (j in seq(n)) {
if (i == j) {
quad_form_matrix[i,j] <- (1/n)
} else {
@bschneidr
bschneidr / alt-design-syntax.R
Created December 11, 2022 14:51
Alternative survey design syntax
describe_design(
srs_stage(id = "PSU_ID", method = "SRS",
stratum = "PSU_STRATUM_ID"),
nonresponse_stage(response_indicator = "PSU_RESPONDENT"),
pps_stage(id = "SSU_ID", method = "PPS",
stratum = "SSU_STRATUM_ID"),
nonresponse_stage(response_indicator = "SSU_RESPONDENT")
)
@bschneidr
bschneidr / chi-square-tests-for-variable-pairs.R
Created October 26, 2023 01:49
library(dplyr)
library(broom)
library(tidyr)
# Make example data
my_data <- data.frame(x = sample(1:5, size = 10, replace = TRUE),
y = sample(1:5, size = 10, replace = TRUE))
# Generate every pair of variables
var_list <- colnames(my_data)
@bschneidr
bschneidr / sample-groups-with-dplyr-and-tidyr.R
Created May 23, 2024 17:59
Sampling groups with 'dplyr' + 'tidyr'
library(dplyr)
library(tidyr)
library(survey)
# Load example data
# 'apipop' has one row per school
data('api', package = 'survey')
# Draw a sample of school districts
@bschneidr
bschneidr / gist:dbc2eb6a1e64d6408ef3281bb2c589d1
Last active September 30, 2024 04:08
Install cmdstanr with WSL and OpenBLAS 
sudo apt-get install libopenblas-dev
sudo apt-get install liblapacke-dev
sudo apt-get install liblapacke
sudo apt-get install libopenblas-serial-dev
sudo apt-get install libopenblas0-serial
sudo apt-get install libopenblas0



  


  




    


          
    

      

        
        

          

              @bschneidr
          

          

            
                bschneidr
                / aoc24-day-1.R
            
            

              Created
              December 2, 2024 23:44
            

          

        

      

        

  

    
    


        



  


  

        

          
          read.table("day-1/input", col.names = c("list_1", "list_2")) |>
        

        

          
            apply(MARGIN = 2, sort) |>
        

        

          
            apply(MARGIN = 1, FUN = \(x) abs(x[2] - x[1])) |>
        

        

          
            sum()
        

  





    


  




    


          
    

      

        
        

          

              @bschneidr
          

          

            
                bschneidr
                / aoc24-day-1.py
            
            

              Created
              December 2, 2024 23:45
            

          

        

      

        

  

    
    


        



  


  

        

          
          import polars as pl
        

        

          
          import pandas as pd
        

        

          
          

        

        

          
          # Read input data as a data frame with two columns
        

        

          
          input = pl.from_pandas(
        

        

          
            pd.read_table("day-1/input", sep = "   ", header = None, 
        

        

          
                          names = ["list_1", "list_2"])
        

        

          
          )
        

        

          
          

        

        

          
          # Sort each column separately,