Displaying 1 result from an estimated 1 matches for "printing_and_protestants_data".
2017 Aug 03
0
Results of vcovCL (sandwich) and of cluster() in Stata
...ant robust cluster(territory)
I'm trying to translate this into R-code doing the following:
library(foreign)
library(dplyr)
library(lmtest)
library(sandwich)
# the data are here:
# https://drive.google.com/file/d/0B_QoCd-1jkVXRGdUMTlkYTNiNGc/edit?usp=sharing
cities <- read.dta("data/Printing_and_Protestants_Data-ReStat.dta")
# we filter the data
cities <- filter(cities, !is.na(pop1500))
cities <- filter(cities, city != "Z?rich" & city != "Mainz" &
city != "Wittenberg")
# the model
m1 <- lm(prot1530 ~ press - 1, data = cities)
# the clustered standard er...