小编Isa*_*iah的帖子

使用Rstudio和Knitr的观星者

我正在努力使用RStudio在knitr中使用stargazer输出.例如,我将下面的代码粘贴到.Rmd文件中,然后单击"编织HTML".[和]之间的第一个块呈现为方程式.第二个街区来自观星者.它仍然是代码.当我将第二个块(少于[和]]粘贴到Sweave文件中然后单击"编译为PDF"时,代码将呈现为表格.我安装了MikTex和Stargazer的第3版.

在单击编译PDF时,将stargazer或xable表插入knitr文档的答案 适用于Sweave文件(Rnw).在Rmd文件中,单击Knit HTML时不会呈现tex.

如何将stargazer输出放入Rmd文件,以便Knit HTML将乳胶代码转换为表格?(我是Latex的新手,我不确定我可以删除哪些代码,所以请为长段道歉.)

\[
\begin{aligned}
\dot{x} & = \sigma(y-x) \\
\dot{y} & = \rho x - y - xz \\
\dot{z} & = -\beta z + xy
\end{aligned}
\]

\[
\documentclass{article}

\begin{document}


% Table created by StarGazer v.3.0.1 by Marek Hlavac, Harvard University. E-mail: hlavac at fas.harvard.edu
% Date and time: Sun, Feb 03, 2013 - 11:34:52 AM
\begin{table}[htb] \centering 
  \caption{} 
  \label{} 
\footnotesize 

\begin{tabular}{@{\extracolsep{5pt}}lc} 
\\[-1.8ex]\hline 
\hline \\[-1.8ex] 
 & \multicolumn{1}{c}{\textit{Dependent variable:}} \\ 
\cline{2-2} 
\\[-1.8ex] & Rate …
Run Code Online (Sandbox Code Playgroud)

r knitr stargazer

13
推荐指数
4
解决办法
3万
查看次数

setInternet2中的错误(TRUE):使用!= NA已不存在

自从我升级到R版本3.3.1后,当我尝试使用q()退出R时,我收到以下消息:"setInternet2中的错误(TRUE):使用!= NA已解除"

这发生在每个出口.

我没有使用过setInternet2,它不在.Rprofile或.Rprofile.site中.我没有.Renviron.site文件.

我检查了我使用正确的安装文件:

file.path(getwd(),".Rprofile")
file.path(getwd(),".Rprofile.Site")
file.path(getwd(),".Renviron.site")
Run Code Online (Sandbox Code Playgroud)

我正在使用R版本3.3.1(2016-06-21)和Windows 7 x64 Service Pack 1

r

7
推荐指数
1
解决办法
2986
查看次数

在模型列表上使用lapply

我已经生成了一个模型列表,并希望创建一个汇总表.

例如,这里有两个模型:

x <- seq(1:10)
y <- sin(x)^2
model1 <- lm(y ~ x)
model2 <- lm(y ~ x + I(x^2) + I(x^3))
Run Code Online (Sandbox Code Playgroud)

和两个公式,第一个从公式的组成部分生成方程

get.model.equation <- function(x) {
  x <- as.character((x$call)$formula)
  x <- paste(x[2],x[1],x[3])
}
Run Code Online (Sandbox Code Playgroud)

第二个生成模型名称作为字符串

get.model.name <- function(x) {
  x <- deparse(substitute(x))
}
Run Code Online (Sandbox Code Playgroud)

有了这些,我创建了一个汇总表

model.list <- list(model1, model2)
AIC.data <- lapply(X = model.list, FUN = AIC)
AIC.data <- as.numeric(AIC.data)
model.models <- lapply(X = model.list, FUN = get.model)
model.summary <- cbind(model.models, AIC.data)
model.summary <- as.data.frame(model.summary)
names(model.summary) <- c("Model", "AIC")
model.summary$AIC <- …
Run Code Online (Sandbox Code Playgroud)

r lapply

3
推荐指数
1
解决办法
2769
查看次数

使用tidytext和扫帚,但没有找到LDA_VEM整洁

tidytext书中有一些示例,主题模型更加整洁:

library(tidyverse)
library(tidytext)
library(topicmodels)
library(broom)

year_word_counts <- tibble(year = c("2007", "2008", "2009"),
+                            word = c("dog", "cat", "chicken"),
+                            n = c(1753L, 1157L, 1057L))

animal_dtm <- cast_dtm(data = year_word_counts, document = year, term = word, value = n)

animal_lda <- LDA(animal_dtm, k = 5, control = list( seed = 1234))

animal_lda <- tidy(animal_lda, matrix = "beta")

# Console output
Error in as.data.frame.default(x) : 
  cannot coerce class "structure("LDA_VEM", package = "topicmodels")" to a data.frame
In addition: Warning message:
In tidy.default(animal_lda, matrix …
Run Code Online (Sandbox Code Playgroud)

r broom tidytext

2
推荐指数
1
解决办法
1551
查看次数

标签 统计

r ×4

broom ×1

knitr ×1

lapply ×1

stargazer ×1

tidytext ×1