使用不同的来源，我编写了一个小函数，该函数创建一个包含标准误差、t 统计量和标准误差的表，这些表根据线性回归模型后的组变量“集群”进行聚类。代码如下
cl1 <- function(modl,clust) {
# model is the regression model
# clust is the clustervariable
# id is a unique identifier in ids
library(plm)
library(lmtest)
# Get Formula
form <- formula(modl$call)
# Get Data frame
dat <- eval(modl$call$data)
dat$row <- rownames(dat)
dat$id <- ave(dat$row, dat[[deparse(substitute(clust))]], FUN =seq_along)
pdat <- pdata.frame(dat,
index=c("id", deparse(substitute(clust)))
, drop.index= F, row.names= T)
# # Regression
reg <- plm(form, data=pdat, model="pooling")
# # Adjustments
G <- length(unique(dat[, deparse(substitute(clust))]))
N <- length(dat[,deparse(substitute(clust))])
# # Resid …我有一台带Sierra的新笔记本电脑.我使用Time Machine从我的旧Mac上带来了我的应用程序,其中大多数工作正常.然而,Macvim消失了.我试着安装它.我尝试通过安装Homebrew和brew安装macvim来重新安装它.安装看起来很成功,但我仍然找不到也没有使用macvim.我想这又是一个问题,应用程序文件不再在/ usr /但在Library /中,但老实说,我有点迷失如何弄清楚这是否是问题以及如何解决它.
你能帮我解决这个问题吗？
我如何安装Homebrew:
ruby -e "$(curl -fsSL https://raw.githubusercontent.com/Homebrew/install/master/install)"
我是如何安装Macvim的:
brew install macvim
这是我面临的问题的可重现示例。我正在尝试在多个阶段使用 ggplot2 创建地图。这是我面临的问题。考虑border具有墨西哥边境的美国各州border.county的多边形以及这些州的县的多边形的数据。以下代码允许您获取数据：
library(maps)
library(ggmap)
library(ggplot2)
USA <- get_googlemap(center = 'usa', zoom = 4,
style = 'administrative|element:labels|visibility:off')
us.df <- map_data("state")
border <- subset(us.df,
region %in% c("california","arizona","new mexico","texas"))
counties <- map_data("county")
border.county <- subset(counties,
region %in% c("california","arizona","new mexico","texas"))
现在我想创建一张地图，以谷歌地图的地图为背景，带有州多边形和县边界。如果我执行以下操作，它会很好地工作：
Allmap <- ggmap(USA) +
geom_polygon(aes(x = long, y = lat, fill = region, group = group),
data=border, color = "white") +
geom_polygon(aes(x = long, y = lat, group = group),
data=border.county, fill=NA, color="red")
现在，如果我想分多个阶段创建此地图，就会遇到问题。我只想要背景信息的县边界（作为一种“重复主题”），我将创建多个地图，其中包含州级不断变化的信息。所以我用县创建了“背景图”，效果很好：
Countmap <- ggmap(USA) + …