小编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% …

我用"rworldmap"包制作世界地图.我正在使用一个函数来访问UN Comtrade的交易数据.

我编辑了我原来的问题所以我可以展示一个我正在做的事情的真实例子.这是我可以制作的地图:

library(rjson)
library(rworldmap)

get.Comtrade <- function(url="http://comtrade.un.org/api/get?", maxrec=50000,
                         type="C", freq="A", px="HS", ps="now", r, p, rg="all",
                         cc="TOTAL", fmt="json") {
  string <- paste(url
                  , "max=", maxrec,"&"  # maximum no. of records returned
                  , "type=", type, "&"  # type of trade (c=commodities)
                  , "freq=", freq, "&"  # frequency
                  , "px=", px, "&"      # classification
                  , "ps=", ps, "&"      # time period
                  , "r=", r, "&"        # reporting area
                  , "p=", p, "&"        # partner country
                  , "rg=", rg, "&"      # trade flow …

Bookdown 有许多交叉引用功能。但是，我无法找到解决问题的方法。

假设您footnote[^1]在文档正文中有一个。然后将脚注定义为：

[^1]: Here is my footnote.

现在假设您想要交叉引用该脚注。例如，参见脚注\@ref(label)。我尝试了一些不同的选项，但没有成功（例如{#label}在脚注定义之后使用[^1]或\@ref(^1)就像^1标签一样使用）。

也许有一个简单的解决方案，但我错过了。或者可能需要有特定的图形\@ref(eq:label)和表格结构\@ref(tab:label)。比如我们可以在脚注的定义中\@ref(foot:label)插入{#label}或。(#foot:label)喜欢

[^1]: Here is my footnote.{#lab}

或者

[^1]: Here is my footnote.(#foot:lab)

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

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

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

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

我对R非常陌生，因此如果这是一个非常基本的问题，我会提前道歉。

我正在尝试绘制一个显示出水和悬浮泥沙量（SSL）的图表。但是，我想弄清楚条形图表示放电，而线图表示SSL。我有两个想法：

1. 为放电标签和SSL分别上色以分别对​​应于条形图和折线图，因此读者可以直观地知道哪个标签属于哪个标签。但是ggplot2不允许我这样做，因为它将以相同的颜色为两个y轴着色。

2. 建立一个图例，清楚地表明红线属于SSL，蓝框图属于排放。有一个帖子在这里它做类似的事情，但我似乎无法实现它。如果有人可以帮助我，我将不胜感激。

这就是我现在的图。

这是我的脚本：

library(ggplot2)
library(gridExtra)
library(RColorBrewer)
library(tibble)

P_Discharge <- Pyay$Mean.monthly.discharge
P_MaxTemp <- Pyay$Mean.monthly.max.temperature
P_MinTemp <- Pyay$Mean.monthly.minimum.temperature
P_Rain <- Pyay$Max.monthly.rainfall
P_SSL <- Pyay$Mean.suspended.sediment.load

#reorderingthemonths
Pyay$Month <- factor(Pyay$Month, 
                     levels=c("Jan", "Feb", "Mar", "Apr", "May", "Jun", 
                              "Jul", "Aug", "Sep", "Oct", "Nov", "Dec"))
#PlottingdischargeandSSL
Pgraph1 <- ggplot(Pyay, aes(x=Month, group=2))
Pgraph1 <- Pgraph1 + geom_bar(aes(y=P_Discharge), stat="identity", fill="blue")
Pgraph1 <- Pgraph1 + geom_line(aes(y=P_SSL), colour="red", size=1) +
  labs(y=expression(Q/(m^{3}))) +
  labs(x="")

#addingsecondaxis
Pgraph1 <- Pgraph1 + scale_y_continuous(sec.axis=sec_axis(~., name=expression(
  Suspended~sediment~load~(10^{6}~t))))

#colouringaxistitles
Pgraph1 <- Pgraph1 + …