小编luc*_*ano的帖子

我有一个如下所示的数据框:

df <- data.frame(Site=rep(paste0('site', 1:5), 50),
           Month=sample(1:12, 50, replace=T),
           Count=(sample(1:1000, 50, replace=T)))

我想删除所有网站上计数总是<每月最大月数的5%的网站.

所有网站的最高月度计数为:

library(plyr)
ddply(df, .(Month), summarise, Max.Count=max(Count))

如果将1的计数分配给site5,则其计数始终<所有站点的最大月计数的5%.因此我希望删除site5.

df$Count[df$Site=='site5'] <- 1

但是,在为site2分配新值后,其中一些计数<最大月度计数的5%,而其他计数> 5%.因此我不希望删除site2.

df$Count[df$Site=='site2'] <- ceiling(seq(1, 1000, length.out=20))

如何将数据框子集化以删除计数总是<每月最大计数的5%的任何网站？如果问题不清楚,请告诉我,我会修改.

我正在使用 knitr 在 Markdown (.Rmd) 文件中嵌入 R 代码。我正在使用 Pandoc 将文件转换为 PDF 文件。我的 .Rmd 文件的内容如下所示：

Report
========================================================

This is my report. Some analysis is included in LINK TO CHUNK NAMED MY.ANALYSIS HERE


```{r my.analysis, echo=FALSE}

summary(cars)

```

它在哪里说 LINK TO CHUNK NAMED MY.ANALYSIS HERE，我能否在输出的 PDF 中提供名为 my.analysis 的代码块的链接？

我相信 OP 在这里问了一个类似但略有不同的问题：图形标题、使用 knitr 和降价到 html 的引用

我在两个向量中包含了很长的单词列表

第一个向量看起来像这样:

x <- c("considerably", "much", "far")

第二个向量看起来像这样:

y <- c("higher", "lower")

我需要一个返回的向量,它列出了每个向量的单词的可能组合.使用x和y,我需要返回此向量

[1] "considerably higher" "considerably lower"  "much higher"         "much lower"         
[5] "far higher"          "far lower"

因此,向量x中的单词必须位于向量y中的单词之前.有这么快的方法吗？

我已经实现了这个混合效果模型:

mtcarsSub <- mtcars[,c("wt", "drat", "cyl")]
library(lme4)
mtcarsME <- lmer(drat ~ (1|cyl) + wt, data=mtcarsSub)

我现在想要从模型中添加预测mtcarsSub.我可以添加这些像这样的新变量:

mtcarsSub$fixed.effect <- predict(mtcarsME)
mtcarsSub$random.effect.cyl4 <- mtcarsMEFixed + ranef(mtcarsME)$cyl["4",]
mtcarsSub$random.effect.cyl6 <- mtcarsMEFixed + ranef(mtcarsME)$cyl["6",]
mtcarsSub$random.effect.cyl8 <- mtcarsMEFixed + ranef(mtcarsME)$cyl["8",]

请注意,当从随机效果添加预测时mtcarsSub,我重复自己3次.如何以编程方式添加随机效果的预测,可能使用函数并且可能在一行中？

scale似乎是一个有用的函数，但它从数字向量生成矩阵，如下所示：

class(scale(mtcars$drat))
# [1] "matrix"

是否有一个函数可以生成数字向量而不必执行类似的操作as.numeric(scale(x))？

我有一个.Rnw文件,如下所示:

\documentclass[a4paper,11pt]{article}
\begin{document}
\title{}
\author{me}
\date{\today}
\maketitle

\section{Header}

<<table_mtcars, results = "asis">>=

stargazer(rbind(mtcars, mtcars, mtcars), summary = F)

@
\FloatBarrier

\end{document}

我可以knitr::knit用来生成.tex文件,然后将其转换为PDF.结果PDF如下所示:

请注意表格如何从第2页的末尾消失,并且不再在第3页重新出现.如何使用观星者将表格继续在PDF的第3页上？

我正在使用ggmap包绘制地图.要通过互联网下载地图,我可以使用以下代码:

library(ggmap)
get_map(location = c(-1.81, 55.655), zoom = 12, maptype = "hybrid")

有没有办法避免通过互联网下载地图,而是从本地文件夹导入.png文件？或者换句话说,下载地图一次,缓存.png,然后从本地文件夹导入.png？我的连接速度很慢,我不断重新加载相同的基本地图浪费了宝贵的时间.

在LaTeX中,可以使用以下代码制作一个简短的标题(将出现在图表列表中):

  \caption[Short caption.]{Long caption.}

随着stargazer在RI试过这样:

stargazer(mtcars, summary = FALSE, title="[Short caption]Long caption")

这会为字幕生成此输出:

  \caption{[Short caption]Long caption} 

我怎样才能制作一个简短的标题stargazer？

像这样命名一个列表很好用：

list("mtcars" = mtcars)

但是当我尝试使用时paste，出现错误：

list(paste("mtcars") = mtcars)

Error: unexpected '=' in "list(paste("mtcars") ="

paste命名列表时如何使用？

考虑这两个t检验的P值

set.seed(1)
x <- c(rnorm(50,1), rnorm(50, 2))
y <- (c(rep("a", 50), rep("b", 50)))

t.test(x ~ y)$p.value

[1] 1.776808e-07

set.seed(2)
x <- c(rnorm(50), rnorm(50))
y <- (c(rep("a", 50), rep("b", 50)))

t.test(x ~ y)$p.value

[1] 0.3922354

第一个P值<0.001,第二个值> 0.001.如果P值恰好<0.001,我可以得到R将P值打印为<0.001.如果P值恰好大于0.001,我可以让R打印出实际的P值吗？因此,第一次t检验"<0.001"的结果应打印为,第二次t检验的结果应打印为0.3922354.

我正在使用knitr我的论文将R代码转换为乳胶.我的规则是只有P值> 0.001才应打印为实际P值.