Thomas Lumley tslumley
😐
tslumley
/ multiresponse.Rmd
Created
July 23, 2019 04:00
Some ideas towards a multiple-response class
This file contains hidden or 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
| --- | |
| title: "Multiple response" | |
| author: "Thomas Lumley" | |
| date: "7/23/2019" | |
| output: html_document | |
| --- | |
| ```{r setup, include=FALSE} | |
| knitr::opts_chunk$set(echo = TRUE) | |
| ``` |
This file contains hidden or 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
| d<-strsplit("Northland | |
| No Motor Vehicle 3792 4077 | |
| One Motor Vehicle 20229 22161 | |
| Two Motor Vehicles 19578 20652 | |
| Three or More Motor Vehicles 7398 7407 | |
| Waitemata | |
| No Motor Vehicle 8871 8658 | |
| One Motor Vehicle 54753 57492 | |
| Two Motor Vehicles 66411 72768 | |
| Three or More Motor Vehicles 29757 33219 |
This file contains hidden or 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
| setClass("ggthing", representation(content="call")) | |
| setClass("ggthings",representation(contents="list")) | |
| bracket<-function(e1) NULL | |
| setGeneric("bracket") | |
| setMethod("bracket",c("ggthing"), | |
| function(e1){ | |
| new("ggthings",contents=list(e1@content)) |
This file contains hidden or 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
| nnot<-100 | |
| ndis<-50 | |
| tpr<-0.9 | |
| fpr<-0.5 | |
| rr<-replicate(100,{ | |
| tpos <-rbinom(1,ndis,tpr) | |
| fpos<-rbinom(1,nnot,fpr) |
This file contains hidden or 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
| use zip-merged.dta | |
| svyset SDMVPSU [pw=WTINT2YR], stra(SDMVSTRA) | |
| svy: zip malepartners RIDAGEYR i.RIDRETH1 DMDEDUC, infl(RIDAGEYR i.RIDRETH1 DMDEDUC) |
This file contains hidden or 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
| --- | |
| title: "Rogue polls" | |
| author: "Thomas Lumley" | |
| date: "27/07/2020" | |
| output: word_document | |
| --- | |
| ```{r setup, include=FALSE} | |
| knitr::opts_chunk$set(echo = TRUE) | |
| ``` |
tslumley
/ proflike.R
Created
August 5, 2020 06:02
Conditional profile likelihood for balanced accuracy
This file contains hidden or 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
| llprofci<-function(tpos,fpos,tneg,fneg){ | |
| ndis<-tpos+fneg | |
| nnot<-tneg+fpos | |
| lljoint<-function(a,d){ | |
| k<-length(a) |
This file contains hidden or 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
| svygofchisq<-function(formula, p, design,...){ | |
| means<-svytotal(formula, design,...) | |
| rval<-chisq.test(means,p=p) | |
| nm<-names(coef(means)) | |
| ncat<-length(coef(means)) | |
| means<-svycontrast(means, list(N=rep(1,ncat)), add=TRUE) | |
| pN<-split(cbind(p,0*diag(p)), paste0("p_",nm)) | |
| names(pN)<-paste0("p_",nm) | |
| means<-svycontrast(means, pN,add=TRUE) | |
| for(i in 1:length(nm)){ |
This file contains hidden or 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
| library(dplyr) | |
| library(rlang) | |
| ma<-read.csv("~/Downloads/weekly-covid-19-municipality-vaccination-report-5-20-2021 (1)/Age - municipality-Table 1.csv",skip=1) | |
| ma %>% | |
| mutate(pop = as.numeric(gsub(",","",Population)), vax=as.numeric(gsub(",","",Fully.vaccinated.individuals))) -> ma_all | |
| ma_all %>% | |
| filter(Age.Group %in% c("65-74 Years","75+ Years")) %>% | |
| select(County,Town,Age.Group,pop,vax) %>% | |
| group_by(County,Town) %>% |
This file contains hidden or 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
| library(rvest) | |
| library(dplyr) | |
| library(DHBins) | |
| a<-read_html("https://www.health.govt.nz/our-work/diseases-and-conditions/covid-19-novel-coronavirus/covid-19-data-and-statistics/covid-19-vaccine-data") | |
| a %>% | |
| html_node("body") %>% | |
| xml2::xml_find_all("//table") %>% nth(n=3) %>%html_text() ->b | |
| gsub("\t"""strsplit(b"\n")[[1]]fixed=TRUE)->d |