View R
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
# Hi Glen. Here is comes: | |
countries <- unique(country) | |
everything <- sapply(countries, function(cy) { | |
parties <- unique(party[country==cy]) # all parties | |
this <- sapply(parties, function(z) sapply(1990:2013, function(y) mav(as.numeric(sapply(seq(y-3,y+3), function(x) experts.raw[country==cy & party==z & year==x]))))) | |
rownames(this) <- 1990:2013 | |
colnames(this) <- parties | |
return(this) | |
}) | |
print(everything ) |
View gist:91d0e6d2fc45444336da
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
<NotepadPlus> | |
<UserLang name="Markdown" ext="markdown mdown mkdn mdwn mkd md Rmd" udlVersion="2.1"> | |
<Settings> | |
<Global caseIgnored="yes" allowFoldOfComments="no" foldCompact="no" forcePureLC="2" decimalSeparator="0" /> | |
<Prefix Keywords1="yes" Keywords2="yes" Keywords3="no" Keywords4="yes" Keywords5="no" Keywords6="no" Keywords7="no" Keywords8="no" /> | |
</Settings> | |
<KeywordLists> | |
<Keywords name="Comments">00# 00## 00### 00#### 00##### 00###### 01 02 03[ 03![ 03<!-- 04] 04] 04--></Keywords> | |
<Keywords name="Numbers, prefix1"></Keywords> | |
<Keywords name="Numbers, prefix2"></Keywords> |
View fade-coefplot.R
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
library(denstrip) | |
fade <- function(x, labels=names(coef(x)), expo=FALSE, xlab="", ylab="", bty="n", ...) { | |
# argument: a regression, additional arguments passed on to plot() and text() | |
coe <- summary(x)$coefficients[,1] # extract coefficients | |
cse <- summary(x)$coefficients[,2] # standard errors | |
len <- length(coe) # how many coefficients (without intercept) | |
if(expo == TRUE) { # exponential form | |
coe <- exp(coe) | |
cse <- exp(cse) | |
} |
View splitworddocuments
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
Attribute VB_Name = "Module1" | |
Sub SplitNotes(delim As String, strFilename As String) | |
Dim doc As Document | |
Dim arrNotes | |
Dim I As Long | |
Dim X As Long | |
Dim Response As Integer | |
arrNotes = Split(ActiveDocument.Range, delim) | |
Response = MsgBox("This will split the document into " & UBound(arrNotes) + 1 & " sections. Do you wish to proceed?", 4) | |
If Response = 7 Then Exit Sub |
View plzcanton-tools.R
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
convid <- function(ID) { | |
# function to convert cantonal ID, v.1.0 (12 Nov 2013) didier.ruedin@wolfson.oxon.org | |
clabels <- c("ZH", "BE", "LU", "UR", "SZ", "OW", "NW", "GL", "ZG", "FR", "SO", "BS","BL", "SH", "AR", "AI", "SG", "GR", "AG", "TG", "TI", "VD", "VS", "NE", "GE", "JU") | |
id <- ifelse(is.numeric(ID), clabels[ID], match(ID, clabels)) | |
return(id) | |
} | |
plzvalid <- function(PLZ) { | |
# check if a number is a valid postcode; v.1.1 (1 Feb 2013) didier.ruedin@wolfson.oxon.org | |
return(!is.na(plzcanton(PLZ))) |
View plzcanton.R
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
plzcanton <- function(PLZ, format=1) { | |
# convert Swiss postcode to canton; v.0.4 (20 Nov 2020) didier.ruedin@unine.ch | |
# format 1 = numeric | |
# format 2 = string | |
# DATA: BFS canton numbers: | |
clabels <- c("ZH", "BE", "LU", "UR", "SZ", "OW", "NW", "GL", "ZG", "FR", "SO", "BS","BL", "SH", "AR", "AI", "SG", "GR", "AG", "TG", "TI", "VD", "VS", "NE", "GE", "JU") | |
# DATA: postcodes: | |
ZH <- c(8970, 8955, 8954, 8953, 8952, 8951, 8942, 8934, 8933, 8932, 8926, 8925, 8915, 8914, 8913, 8912, 8911, 8910, 8909, 8908, 8907, 8906, 8904, 8903, 8902, 8901, 8833, 8825, 8824, 8820, 8816, 8815, 8813, 8812, 8810, 8805, 8804, 8803, 8802, 8801, 8800, 8714, 8713, 8712, 8708, 8707, 8706, 8704, 8703, 8702, 8700, 8639, 8637, 8636, 8635, 8634, 8633, 8632, 8630, 8627, 8626, 8625, 8624, 8623, 8622, 8621, 8620, 8618, 8617, 8616, 8615, 8614, 8613, 8612, 8610, 8608, 8607, 8606, 8605, 8604, 8603, 8602, 8600, 8548, 8545, 8544, 8543, 8542, 8523, 8499, 8498, 8497, 8496, 8495, 8494, 8493, 8492, 8489, 8488, 8487, 8486, 8484, 8483, 8482, 8479 |
View gist:5842378
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
library(psych) | |
n <- 5000 | |
x <- rnorm(n) # random values | |
r <- sapply(1:n, function(i) alpha(data.frame(x, c(x[1:i],rep.int(0,n-i))))$total$std.alpha) | |
plot(r, type="l", ylim=c(0,1), xlim=c(0,n), axes=FALSE, ylab="Alpha", xlab="Zero-Inflation (%)") | |
axis(2) | |
axis(1, at=seq(from=0,to=n,by=n/5), labels=seq(from=0,to=100,by=20)) |
View namav
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
namav <- function(x,k=3){ | |
x <- c(rep(NA, k),x,rep(NA,k)) # add NA on both sides | |
n <- length(x) | |
return(sapply((k+1):(n-k), function(i) sum(x[(i-k):(i+k)],na.rm=TRUE)/(2*k+1-sum(is.na(x[(i-k):(i+k)]))))) | |
} |
View ethnic-and-racial-studies.csl
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
<?xml version="1.0" encoding="utf-8"?> | |
<style xmlns="http://purl.org/net/xbiblio/csl" class="in-text" version="1.0" page-range-format="chicago" demote-non-dropping-particle="sort-only"> | |
<info> | |
<title>Ethnic and Racial Studies</title> | |
<id>http://www.zotero.org/styles/ethnic-and-racial-studies</id> | |
<link href="http://www.zotero.org/styles/ethnic-and-racial-studies" rel="self"/> | |
<link href="www.tandfonline.com/action/authorSubmission?journalCode=rers20&page=instructions" rel="documentation"/> | |
<author> | |
<name>Didier Ruedin</name> | |
<email>didier.ruedin@wolfson.oxon.org</email> |
View as_pdf.R
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
as_pdf <- function(x){ | |
require(tools) | |
fname <- tempfile(pattern = "texview-", tmpdir = tempdir(), | |
fileext = ".tex") | |
header <- "\\documentclass{article} | |
\\usepackage[margin=10pt,font=small,labelformat=empty, | |
labelsep=none]{caption} | |
\\usepackage{dcolumn} |
NewerOlder