小编jay*_*.sf的帖子

只是想知道是否有人遇到过他/她需要从非常高维多元正态分布(比如维度= 10,000)中随机抽取的问题,因为包的rmvnorm功能mvtnorm是不切实际的.

我知道这篇文章有一个包Rcpp的dmvnorm功能实现mvtnorm,所以我想知道是否存在等效的东西rmvnorm？

我一直在Objective-C中使用这个宏:

#define     RGBA(r, g, b, a) [UIColor colorWithRed:(r)/255.0 green:(g)/255.0 blue:(b)/255.0 alpha:(a)]

我试图找出如何在swift中获得最接近的东西.有任何想法吗？

我通常使用mfx包和logitmfx函数生成logit模型边际效应.然而,我目前使用的调查具有权重(由于某些人群中的过采样,其对样本中DV的比例具有很大影响)并且logitmfx似乎没有任何方式来包括权重.

我用svyglm为模型拟合如下:

library(survey)

survey.design <- svydesign(ids = combined.survey$id,
                                        weights = combined.survey$weight,
                                            data = combined.survey)

vote.pred.1 <- svyglm(formula = turnout ~ gender + age.group + 
                                    education + income, 
                                 design = survey.design)
summary(vote.pred.1)

如何从这些结果中产生边际效应？

出于发布的目的，我经常需要我的作品的 PDF 和 HTML 版本，包括回归表，并且我想使用 R Markdown。对于 PDF，stargazer和texreg软件包可生成精美的表格。现在，在尝试生成同样有吸引力的 HTML 输出时，我面临着不同的问题。

1. 这两种 HTML 输出方法在注释中都缺少重要的星号。由于它们是自动生成的，我不知道如何逃脱它们。（我认为这可能是一个小问题，因此我不想将其分成单独的问题。） 注意：子问题已在此处得到解答。

2. 在创建明确的输出之前，我经常必须更改数据或进行一些格式化。type='html'我发现总是在手动之间切换选项非常烦人type='pdf'。我想知道是否有更可行的方法来组合 html/pdf 输出texreg，例如/中的个案切换stargazer与整洁的输出？

我尝试了有希望的pander解决方案，但自 2014 年以来它似乎不再起作用了。而且它pixiedust也不是很令人满意，它最终变得有点手动，而不完全是我想要的。另一个示例似乎仅涉及普通表。

非常感谢任何帮助，谢谢！

knitr以下是我在 HTML 和 PDF 中尝试的总结：

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```
```{r table, results = "asis"}
library(car)
lm1 <- lm(prestige ~ income + education, data=Duncan)

## html
# stargazer
library(stargazer)
stargazer(lm1, type="html", …

我想使用 R Markdown 在 HTML 文档中显示回归输出。我尝试了texreg和stargazer包。我现在的问题是，在笔记中我无法将重要的星星带入生活。由于自动生成，我似乎无法摆脱它们。我一直对这个和这个感到困惑，但没有成功。我错过了什么？非常感谢！！

这是一些代码：

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```
```{r data}
library(car)
lm1 <- lm(prestige ~ income + education, data=Duncan)
```
## with STARGAZER
```{r table1, results = "asis", message=FALSE}
library(stargazer)
stargazer(lm1, type="html", notes="stargazer html 1") # nothing
stargazer(lm1, type="html", notes="stargazer html 2", star.char = "\\*") # nothing, even gone in table
```
## with TEXREG
```{r table2, results = "asis", message=FALSE}
library(texreg)
htmlreg(lm1, custom.note="%stars. htmlreg") # …

我无法在R中复制Stata margins命令的特定用例： margins var1, over(var2) 我一直在尝试使用marginsR中的包进行复制。

为了提供可重现的示例，我使用了mtcars数据集并将其从R导出到Stata中，因此我们在两个程序中都使用了相同的数据集：

R代码：

library(foreign)
library(margins)
write.dta(mtcars, “mtcars.dta")

Stata代码：

use "mtcars.dta", clear

在两个程序中创建示例线性回归模型

Stata代码：

quietly regress mpg cyl i.am c.wt##c.hp

R代码：

x <- lm(mpg ~ cyl + factor(am) + hp * wt, data = mtcars)

两个程序之间的模型输出（未显示）相同

比较模型中每个变量的平均边际效应表

Stata代码和输出：

margins, dydx(*)

Average marginal effects                          Number of obs   =         32
Model VCE: OLS

Expression   : Linear prediction, predict() dy/dx w.r.t. : cyl 1.am wt hp

------------------------------------------------------------------------------
             |            Delta-method
             |      dy/dx   Std. Err.      t    P>|t|     [95% …

我很想在Rmarkdown中以Tufte风格进行“光束演示” 。在那里，我可以把在YAML头两种output: beamer_presentation或

output:
  tufte::tufte_html: default    
  tufte::tufte_handout:
    citation_package: natbib
    latex_engine: xelatex

有什么办法可以结合起来吗？还是任何其他简单的解决方案？

注意：我之前在tex.stackexchange上问过这个问题，没有解决方案。

我在他的书的网站上读过Hadley Wickham的Advanced R. 我在R中找到了关于替换函数的问题.根据他的书给出了以下结果.

library(pryr)
x <- 1:10
address(x)
#> [1] "0x103945110"

x[2] <- 7L
address(x)
#> [1] "0x103945110"

如果我们只是替换x的第二个元素,他认为x的地址不会改变.但是,当我这样做时,x的物理地址实际上已经改变了.那么,有谁告诉我为什么？

我有一个列表L（来自矢量分割）。

L <- strsplit(c("1 5 9", "", "3 7 11", ""), " ")

# [[1]]
# [1] "1" "5" "9"
# 
# [[2]]
# character(0)
# 
# [[3]]
# [1] "3"  "7"  "11"
# 
# [[4]]
# character(0)

当我rbind按如下所示进行普通操作时，我将丢失所有character(0)行。

do.call(rbind, L)
#      [,1] [,2] [,3]
# [1,] "1"  "5"  "9" 
# [2,] "3"  "7"  "11"

我是否总是必须执行lapply以下操作或错过了什么？

do.call(rbind, lapply(L, function(x) 
    if (length(x) == 0)  rep("", 3) else x))
#      [,1] [,2] [,3]
# …

我有一个看似很简单的问题，与之相比，我需要更快的R实现。

我为此示例初始化随机种子和尺寸：

set.seed(1)
d1<-400
d2<-20000
d3<-50

我有一个矩阵X，尺寸为d1 x d2：

X<-as.data.frame(matrix(rnorm(d1*d2),nrow=d1,ncol=d2))
rownames(X)<-paste0("row",1:nrow(X))
colnames(X)<-paste0("col",1:ncol(X))

以及具有d1行索引的向量u ：

u<-sample(rownames(X),nrow(X),replace=TRUE)

我也有一个矩阵C，其命名的行和尺寸为d3 x d2：

C<-matrix(rnorm(d3*d2),nrow=d3,ncol=d2)
rownames(C)<-sample(rownames(X),nrow(C),replace=FALSE)

现在，通过以下非常慢的循环，我用匹配的X行的总和填充矩阵C：

system.time(
    for(i in 1:nrow(C)){
        indexes<-which(u==rownames(C)[i])
        C[i,] <- colSums(X[indexes,])
    }
)

在我的PC上，此操作大约需要11.5秒，但是我确信可以通过避免for循环来加快速度。有任何想法吗？非常感谢！