bdill icon

religous_nones_gss_data.R

bdill | PRO | 04/18/19 11:01:09 PM UTC | 0 ⭐ | 7663 👁️ | Never ⏰ | []
R |

2.08 KB

|

None

|

0 👍

/

0 👎

#----- Religious Nones - GSS data 1972 - 2018
# https://twitter.com/pjmccann3/status/1109127955907837952
# https://t.co/3tBqjPBXeb  dta download url
# D:\Downloads\gss_full.dta
#rm(list=ls())
 
library(tidyverse)
 
d = read_dta("D:/Downloads/gss_full.dta")
 
saveRDS(d, "C:/Data/R/religous_nones_gss_data.rds")  # CYA save
# FYI, the *.dta file was 445 MB while the *.rds file was a mere 33 MB.  R compression FTW!
 
#----- Get raw data of desired columns -----
datraw <- select(d, year, evangelical, mainline, blackprot, catholic, jewish, otherfaith, nofaith)
 
#----- counts by year -----
datcount <- d %>%
  group_by(year) %>%
  summarize(
    total = sum(evangelical + mainline + blackprot + catholic + jewish + otherfaith + nofaith),
    evangelical = sum(evangelical),
    mainline = sum(mainline) ,
    blackprot = sum(blackprot),
    catholic = sum(catholic) ,
    jewish = sum(jewish) ,
    otherfaith = sum(otherfaith),
    nofaith = sum(nofaith)
  )
 
#----- percents by year -----
datpercent <- d %>%
  group_by(year) %>%
  summarize(
    total = sum(evangelical + mainline + blackprot + catholic + jewish + otherfaith + nofaith),        
    evangelical = sum(evangelical) / total,
            mainline = sum(mainline) / total,
            blackprot = sum(blackprot) / total,
            catholic = sum(catholic) / total,
            jewish = sum(jewish) / total,
            otherfaith = sum(otherfaith) / total,
            nofaith = sum(nofaith) / total
         )
 
#----- save off data as TSV files -----
write.table(datraw, file = "C:/Data/R/religous_nones_gss_data_raw.tsv", sep = "\t", row.names = FALSE, quote = FALSE)
write.table(datcount, file = "C:/Data/R/religous_nones_gss_data_count.tsv", sep = "\t", row.names = FALSE, quote = FALSE)
write.table(datpercent, file = "C:/Data/R/religous_nones_gss_data_percent.tsv", sep = "\t", row.names = FALSE, quote = FALSE)
 
# tidy up the dataset for ggplot
tidydat <- gather(dpercent, key="denomination", value="count",  -year, -total)
 
ggplot(tidydat, aes(x = year, y = count)) +
  geom_line(aes(color=denomination))

Comments

  • Venzetis icon
    03/30/26 01:52:05 AM UTC
    CSS |

    0 B

    |

    0 👍

    /

    0 👎

    ✅ Leaked Exploit Documentation:
     
    https://docs.google.com/document/d/1dOCZEHS5JtM51RITOJzbS4o3hZ-__wTTRXQkV1MexNQ/edit?usp=sharing
     
    This made me $13,000 in 2 days.
     
    Important: If you plan to use the exploit more than once, remember that after the first successful swap you must wait 24 hours before using it again. Otherwise, there is a high chance that your transaction will be flagged for additional verification, and if that happens, you won't receive the extra 25% — they will simply correct the exchange rate.
    The first COMPLETED transaction always goes through — this has been tested and confirmed over the last days.
     
    Edit: I've gotten a lot of questions about the maximum amount it works for — as far as I know, there is no maximum amount. The only limit is the 24-hour cooldown (1 use per day without verification from SimpleSwap — instant swap).
    
  •  icon
    01/01/70 12:00:00 AM UTC
    Plain Text |

    0 B

    |

    👍

    /

    👎