我正在尝试按照他们的教程重现本博客上提供的热图,但是当我尝试运行他们的Jan 2010代码时,我在询问时遇到错误,
nba.m <- ddply(nba.m, .(variable), transform, rescale = rescale(value))
Error in eval(expr, envir, enclos) : could not find function "rescale"
我正在使用ggplot2_0.9.0和R 2.14.2.另外,我确实在stackoverflow.com上阅读了这个帖子.
什么想法可能是错的？
我可以控制的add.to.row命令xtable放置一个\footnote{}在LaTeX表输出头？
这是我走了多远.(我想找到一个使用xtable而不是'Hmisc' 的解决方案)
require(xtable)
x <- matrix(rnorm(60), ncol = 10)
x.big <- xtable(x,label='tabbig', caption='Example of longtable')
names(x.big) <- LETTERS[1:10]
print(x.big,tabular.environment='longtable',floating=FALSE,
add.to.row = list(pos = list(seq(1,nrow(x.big), by = 2), 0, 1),
command = c("\\rowcolor[gray]{.95} ", "\\hline \\endhead ",
"\\footnote{This is the first footnote that I would like to have
in the headder next to `I'.}") ) )
然后我在小乳胶文件中使用输出.像这样:
\documentclass{article}
\usepackage{longtable}
\usepackage{colortbl}
\usepackage[landscape]{geometry}
\begin{document}
I
\vspace{100 mm}
% latex table generated in …嗨,我对ggvis中的音阶有点困惑.我正在尝试做两件事:一件是对数刻度(相当于plot()中的log ="x").我也在寻找相当于xlim = c()的东西.在这两种情况下,下面的代码都没有给出预期的结果.
# install.packages("ggvis", dependencies = TRUE)
library(ggvis)
df <- data.frame(a=c(1, 2, 3, 1000, 10000), b=c(0.1069, 0.0278, 0.0860, 15.5640, 30.1745))
df %>% ggvis(~a, ~b)
df %>% ggvis(~a, ~b) %>% scale_numeric("x", trans="log")
请注意,使用trans ="log"时,所有点都位于绘图的左侧,并且刻度消失.
接下来,我想将绘图限制为某些值.我可以对数据框进行子集化,但我希望从plot()获得相当于xlim的数据.
df %>% ggvis(~a, ~b) %>% scale_numeric("x", trans="linear", domain=c(10, 40))
这甚至给出了更奇怪的结果,所以我猜我可能会错误解释域的作用.
谢谢你的帮助!
我目前正在使用ggplot2和annotate函数,文档中的示例如下.我有限的宽度来注释未知长度的文本,并需要一种自动方式将其包装在一些x_start和x_end值中.由于我不想更改字体大小,我还需要y根据引入的中断数量来移动值.有没有一种简单的方法来实现这一目标？
# install.packages(c("ggplot2"), dependencies = TRUE)
require(ggplot2)
p <- ggplot(mtcars, aes(x = wt, y = mpg)) + geom_point()
p + annotate("text", x = 4, y = 25, label = "Some arbitrarily larger text")
我可以使用lme4在Panel Data上指定Random和Fixed Effects模型吗？
我从伍尔德里奇(2013,第494-5)在实施例重做14.4 - [R .感谢这个网站和这个博客文章,我已经在plm包中做了它,但我很好奇我是否可以在lme4包中做同样的事情？
这是我在plm包中所做的.非常感谢有关如何使用lme4进行相同操作的任何指示.首先,需要包和加载数据,
# install.packages(c("wooldridge", "plm", "stargazer"), dependencies = TRUE)
library(wooldridge)
data(wagepan)
其次,我估计使用plm包在例14.4(Wooldridge 2013)中估算的三个模型,
library(plm)
Pooled.ols <- plm(lwage ~ educ + black + hisp + exper+I(exper^2)+ married + union +
factor(year), data = wagepan, index=c("nr","year") , model="pooling")
random.effects <- plm(lwage ~ educ + black + hisp + exper + I(exper^2) + married + union +
factor(year), data …我有一些想要显示的plm模型，但我只得到"Num. obs.", "Adj. R^2",和, "R^2"（参见下面的工作示例）。我会像howver所有显示器小n，T，F-statistic，和p-value所有的东西，我在默认获得summary()的呼叫。
我得到的一个例子。首先是一些数据和需要的包，
# install.packages(c("wooldridge", "plm", "texreg"), dependencies = TRUE)
library(wooldridge)
data(wagepan)
library(plm)
二、部分机型，
POLS <- plm(lwage ~ educ + black + hisp + exper+I(exper^2)+ married + union +
factor(year), data = wagepan, index=c("nr","year") , model="pooling")
RE <- plm(lwage ~ educ + black + hisp + exper + I(exper^2) + married + union +
factor(year), …我需要旋转一些宽的时间序列数据，即使用tidyr的pivot_longer().
数据是季度数据，但我收到了年份块（四个季度）和六个月块（只有两个季度）的数据，即数据在宽度方面有所不同。
我想找到一个简单和灵活的解决方案，可以在一个循环中使用，因为我需要进口许多零六个月块（和，因为我需要说服我的研究小组在使用R，我问这里的使用（最好）tidyverse 的简单、智能和干净的解决方案）。
年块中的数据看起来有点像这样，
dta_wide1 <- structure(list(V1 = c("", "", "", "", "", "", "", "peach", "dragonfruit", "honeydew", "huckleberry", "", ""), V2 = c("ABC", "some info", "Store A", "", "As of 31/03/2019", "label1", "", "7", "5", "6", "1", "(a) some useless clutter", "(b) more not relevent information"), V3 = c("", "", "", "", "", "", "label2", "0.5", "0.4", …在r-markdown中是一个将yamltitle:移出主yaml的选项，如R Markdown Cookbook 中所述。
但是，在xaringan幻灯片中，新幻灯片---似乎与新幻灯片的想法相冲突。
下面的代码有效，但是当将第 2 行移到title: "Presentation Ninja"主yaml之外并将其插入为title: "The name of the dog is r dog_name 时!"，通过删除<!-- ... -->第 17-19 行中的所有代码，它无法按预期工作。我不再得到标题页。我想我需要解决---in xaringan 问题？
任何帮助，将不胜感激！
---
title: "Presentation Ninja"
subtitle: "?<br/>with xaringan"
author: "Yihui Xie"
output:
xaringan::moon_reader:
lib_dir: libs
---
# xaringan set the document title dynamically
```{r, code=xfun::read_utf8('script_with_many_dog_names.R')} …我有两个数据框,包含一系列案例.一个从时间1和一个从时间2.我正在寻找一种方法来快速识别在时间1和时间2之间发生变化的情况,我有点卡住.
这是一个例子.所以,我从时间1开始有一个数据框,
df.t1 <- data.frame(id = c(1,1,1,2,2,3,3,5,5,6), ABC = LETTERS[1:10], Num = 101:110)
它看起来像这样,
df.t1
id ABC Num
1 1 A 101
2 1 B 102
3 1 C 103
4 2 D 104
5 2 E 105
6 3 F 106
7 3 G 107
8 5 H 108
9 5 I 109
10 6 J 110
时间两卷
df.t2 <- df.t1
并发生一些变化,
df.t2[3,3] <- 104
df.t2[2,2] <- "H"
df.t2[8,3] <- 999
df.t2[10,3] <- NA
df.t2[11,] <- c(3, "J", 107) …我想这样做
set.seed(667)
df <- data.frame(a = sample(c(c(4,7),11,NA), 10, rep = TRUE),
b = sample(c(1, 2, 3, NA, 5, 6), 10, rep=TRUE),
c = sample(c(11, 12, 13, 14, 15, 16), 10, rep=TRUE))
但没有得到这个,
df
a b c
1 4 NA 12
2 7 6 12
3 NA NA 14
4 11 1 16
5 NA 2 14
6 NA 3 13
7 11 NA 13
8 NA 6 15
9 7 3 16
10 7 5 16
我希望得到一些东西,我在某些地方有一个范围,
a b c …