小编jay*_*.sf的帖子

我在 R 中运行下面的代码，它不返回原始十进制数。

as.hexmode(-8192)
"ffffe000"

strtoi(c("ffffe000"))
NA

我想将具有不同值的多个变量的函数应用于列表.我知道如何使用一个变化的变量来做到这一点

sapply(c(1:10), function(x) x * 2)
# [1]  2  4  6  8 10 12 14 16 18 20

但不是两个.我首先手动向你展示我想要的东西(实际上我使用 lapply()但sapply()在SO中更具概要性):

# manual
a <- sapply(c(1:10), function(x, y=2) x * y)
b <- sapply(c(1:10), function(x, y=3) x * y)
c <- sapply(c(1:10), function(x, y=4) x * y)
c(a, b, c)
# [1]  2  4  6  8 10 12 14 16 18 20  3  6  9 12 15 18 21 24 27 30  4  8 12 
# [24]  16 …

我正在创建一个包含一些绘图的 PDF 文件，但我还想在底部包含一些文本消息。由于我无法控制的原因，我无法在该系统上安装乳胶发行版，因此我无法编写 markdown 文件，但必须使用pdf().

当我只是使用print或catpdf 中没有显示任何内容时。我尝试sink()根据此处的答案使用，但这也不起作用：

pdf("filename.pdf")
sink("filename.pdf")
print("message")
sink()
dev.off()

未收到错误消息，但创建的文件没有页面。

有任何想法吗？我正在考虑仅绘制仅文本图的解决方法，但我希望有一个更合理的解决方案。

在矩阵中，例如M1，行是国家，列是年。这些国家在同一年没有观测值。我想找到多年来为我提供最多国家/地区的“最佳”交集。最小年份和最小国家/地区的数量将被预先定义。结果中包括哪些国家并不重要，年份不必连续。

> M1
      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15]
 [1,]   NA   NA   NA 2004   NA 2006   NA 2008 2009    NA  2011  2012    NA    NA    NA
 [2,]   NA 2002   NA 2004   NA   NA 2007   NA   NA  2010  2011    NA  2013  2014    NA
 [3,]   NA   NA   NA 2004 2005 2006 2007 2008 2009    NA    NA  2012  2013    NA  2015
 [4,]   NA 2002   NA 2004 2005 2006 2007 2008   NA  2010  2011 …

例子：

---
output: pdf_document
mainfont: Times New Roman
fontsize: 12pt
spacing: 1.1
geometry: "left=1in,right=1in,top=0.5in,bottom=1in"
vspace: 4pt
---

This line is normal text.

This line uses Garamond font and font size 14.

This line uses arial font and font size 16.

This line is normal text.

我尝试了几种选择，但没有一个有效。请告诉我是否有办法做到这一点。

我想知道为什么as.character(5.0)返回5但as.character(5.1)返回5.1R。我试图通过阅读文档获得答案，但没有运气。

我通常使用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)

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

我无法在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% …

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

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

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

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

我想要那些在刻面网格中带有标签的刻度线。所以我完成了这次尝试，并尝试将其调整为像这样的分面网格图：

定义中断和标签，次要和主要：

range.f <- range(unique(df1$weeks))
minor.f <- 1  # every 1 week, NOTE: range.f[2] should be divisible by minor.f!
major.f <- 5  # every 5 weeks

breaks.f <- seq(range.f[1], range.f[2], minor.f)

every_nth.lt <- function (x, nth) {x[1:nth != 1] <- ""; x}
# (lite version of /sf/answers/2417343141/
# works better for me than `insert_minor()`)

labels.f <- every_nth.lt(sequence(range.f[2]), major.f)

n_minor.f <- major.f / minor.f - 1

正常情节：

library(ggplot2)
p.f <- ggplot(df1, aes(weeks, births)) +
  geom_bar(stat="identity", fill="#F48024") + theme_bw() + …