Homework 4

library(tidyverse)

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.2.1     ✔ readr     2.2.0
## ✔ forcats   1.0.1     ✔ stringr   1.6.0
## ✔ ggplot2   4.0.3     ✔ tibble    3.3.1
## ✔ lubridate 1.9.5     ✔ tidyr     1.3.2
## ✔ purrr     1.2.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

library(readxl)
library(pastecs)

## 
## Attaching package: 'pastecs'
## 
## The following objects are masked from 'package:dplyr':
## 
##     first, last
## 
## The following object is masked from 'package:tidyr':
## 
##     extract

Religious_Landscape <- read_csv("2023-24 RLS Public Use File Feb 19.csv")

## Rows: 36908 Columns: 654
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## dbl (631): P_SUID, YEAR_FROM_CODATE, HAPPY, SATIS_A, SATIS_B, CHNG_A, CHNG_B...
## lgl  (23): REGION, DENOMREC2, FAMILY, PROTFAM, CURREL2, FRMREL2, SPREL, SPRE...
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.

pastecs::stat.desc(Religious_Landscape$RELTRAD)

##      nbr.val     nbr.null       nbr.na          min          max        range 
## 3.690800e+04 0.000000e+00 0.000000e+00 1.100000e+03 9.000000e+05 8.989000e+05 
##          sum       median         mean      SE.mean CI.mean.0.95          var 
## 1.745428e+09 1.000000e+04 4.729132e+04 5.253475e+02 1.029696e+03 1.018624e+10 
##      std.dev     coef.var 
## 1.009269e+05 2.134153e+00

This is the frequency rate that shows how often a religious tradition occurs within the avaible options that were listed on the survey that was sent out to all respondents. The minmal number a for each religious tradiion is one (theoretically zero). As a note, much of my data has “codes” which are assigned numerical values as a way to keep the data clean. I do not know if that affectred this data per say but I wanted to try with it specifically because I wanted charts of religious tradition frequency in my paper. I’ll be interested to see what you say as I may need direction on how to accomplish that properly.

hist(Religious_Landscape$RELTRAD)

Religious_Landscape_CLEAN <- Religious_Landscape %>% drop_na(RELTRAD)

hist(Religious_Landscape_CLEAN$RELTRAD)

Religious_Landscape_CLEAN <- Religious_Landscape_CLEAN %>% mutate(RELTRAD_SQRT=sqrt(RELTRAD))

hist(Religious_Landscape_CLEAN$RELTRAD_SQRT)