小编Prr*_*dep的帖子

让我来说,这可能更适合作为html/css问题,但是(a)我对这些技术并不太熟悉,(b)如果可能的话,我想把它全部留在家庭(R家族)中.

我想用Rmarkdown knitr编写学术风格报告(APA 6类型指南).我已经解决了大部分方面但没有分页.我可以手动设置分页符,例如:

# report

```{r setup, include=FALSE}
# set global chunk options
opts_chunk$set(cache=TRUE)
```
------
## Page 1

```{r plot1}
plot(1:10, 1:10)
```
------
## Page 2

在下面.Rmd我想在n行/ n距离之后以编程方式设置那些中断.所以,让我们说8英寸或140行之后.

# report

```{r setup, include=FALSE}
# set global chunk options
opts_chunk$set(cache=TRUE)
```
Initial Text. Yay!

```{r plot1}
plot(1:10, 1:10)
```

More Text.  Outstanding.  What Hadley's not calling it plyr2?

```{r plot2, fig.width=4, fig.height=4}
plot(1:10, 1:10)
```

`r paste(rep("So much text so little time!", 10000))`

如何在n距离之后以编程方式设置分页符.这与LaTeX …

我正在绘制ROC并测量部分AUC作为生态位模型质量的度量.当我在R工作时,我正在使用ROCR和pROC包.我会选择一个使用,但就目前而言,我只是想看看它们是如何表现的,以及是否能更好地满足我的需求.

令我困惑的一件事是,在绘制ROC时,轴如下:

ROCR

x axis: 'true positive rate' 0 -> 1
y axis: 'false positive rate', 0 -> 1

PROC

x axis: 'sensitivity' 0 -> 1
y axis: 'specificity' 1 -> 0. 

但如果我使用两种方法绘制ROC,它们看起来完全相同.所以我只想确认一下:

true positive rate = sensitivity

false positive rate = 1 - specificity.

这是一个可重复的例子:

obs<-rep(0:1, each=50)
pred<-c(runif(50,min=0,max=0.8),runif(50,min=0.3,max=0.6))
plot(roc(obs,pred))

ROCRpred<-prediction(pred,obs)
plot(performance(ROCRpred,'tpr','fpr'))

我正在学习使用ggplot2,我正在寻找ggplot2能够重现以下base::plot结果的最小代码.我已经尝试了一些东西,它们最终都变得非常长,所以我正在寻找最小的表达式,理想情况下我想在x轴上设置日期(plot下面没有这些日期).

df = data.frame(date = c(20121201, 20121220, 20130101, 20130115, 20130201),
                val  = c(10, 5, 8, 20, 4))
plot(cumsum(rowsum(df$val, df$date)), type = "l")

我已将git版本控件集成到Rstudio(Windows)中的现有项目中.在此之后,我在Rmd文件中进行的每个更改都会弹出一个弹出窗口.我可以看到它是在我在文件中添加的每个命令上执行的git bash.我没有更多细节.

如何将Z - statistics的值作为来自glm对象的向量？例如,我有

fit <- glm(y ~ 0 + x,binomial)

如何Pr(>|z|)以与系数估计相同的方式访问列fit$coef？

是否可以直接使用knitr函数knit2pdf()与R Markdown(Rmd)文件？我已经看过各种教程/课堂笔记似乎表明它可以在这里和这里(Ctrl + F"knit2pdf"中的任何一个).

但是当我拿一个简单的rmd文件(保存为"test.rmd")时

---
title: "knit2pdf test"
author: "A Aaronson"
date: "Thursday, February 19, 2015"
output: pdf_document
---

This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see <http://rmarkdown.rstudio.com>.

When you click the **Knit** button a document will be generated that includes both content as well as the output of any embedded R code chunks within …

我想识别dataframe1中存在的行,这些行基于特定列在dataframe2中不存在.我使用下面的代码来获取所需的信息.

diffId <- anti_join(dat$ID,datwe$ID)

不幸的是,我遇到了一个错误:

UseMethod中的错误("anti_join"):
没有适用于"anti_join"的方法应用于类"factor"的对象

我检查了两个数据帧中所需列的类,结果证明是这样的factor.还试图将列分成单独的变量,假设它可以解决问题,但没有运气!

fac1 <- datwe$ID
fac2 <- dat$ID
diffId <- anti_join(fac2,fac1)

你能分享一下你的想法吗？

谢谢

我有一个闪亮的应用程序(使用navbarPage)有很多选项卡,并希望添加一个侧栏菜单,无论选择哪个选项卡都可以看到.侧栏中的输入值会影响所有选项卡的内容.此外,应该可以隐藏sidebarMenu,因为它在shinydashboard中.

我看到两种可能的方式:

(A)使用shinydashboard并以某种方式添加顶部导航栏或

(B)使用navbarPage并以某种方式添加可隐藏的侧边栏菜单.

(A)使用shinydashboard,最接近我想要的是(简化的MWE):

library("shiny")
library("shinydashboard")

cases <- list(A=seq(50,500, length.out=10), B=seq(1000,10000, length.out=10))

ui <- dashboardPage(
  dashboardHeader(title = "dash w/ navbarMenu"),
  dashboardSidebar(selectizeInput('case', 'Pick a case', selected="A", choices = c("A", "B"), multiple = FALSE), numericInput('num', 'Number', min = 1, max = 10, value = 1, step = 1)),
  dashboardBody(
    tabsetPanel(
      tabPanel(h4("Perspective 1"),
               tabsetPanel(
                 tabPanel("Subtab 1.1", plotOutput("plot11")),
                 tabPanel("Subtab 1.2")
               )),
      tabPanel(h4("Perspective 2"),
               tabsetPanel(
                 tabPanel("Subtab 2.1"),
                 tabPanel("Subtab 2.2")
               ))
    )
  )
)

server <- function(input, output) {
  output$plot11 <- renderPlot({ …

我希望将闪亮的图形输出高度和宽度调整为当前窗口大小.我试过使用下面但没有用.

ShinyUi <- fluidPage(

  # Application title
  titlePanel("title"),

  sidebarLayout(
    sidebarPanel(
      ... inputs ...
    ),

    mainPanel(
          plotlyOutput("distPlot", height = 'auto', width = 'auto')
      )
  ))

ShinyServer <- function(input, output, session) {

   output$distPlot <- renderPlotly({

    p <- ggplot(dataShow, aes(x=dataShow$X, y=dataShow$Y))  + 
geom_point(shape=1, alpha = 0.5, color = "grey50")

    ggplotly(p)

  })

}


# Run the application 
shinyApp(ui = ShinyUi, server = ShinyServer)

您是否知道可能在服务器功能中使用的任何其他选项而不是上述UI功能使用？

较小的窗口:

扩展窗口:

我有一个datwe37列的数据框.我有兴趣将第23到35列中的整数值(1,2,99)转换为字符值('是','否','NA').

datwe$COL23 <- sqldf("SELECT CASE COL23 WHEN 1 THEN 'Yes'
                                        WHEN 2 THEN 'No'
                                        WHEN 99 THEN 'NA'
                                   ELSE 'Name ittt' 
                              END as newCol
                              FROM datwe")$newCol

我一直在使用上面的sqldf语句分别转换每一列.我想知道是否有其他智能方法可以做到这一点,也许应用功能？

如果您需要任何可重现的数据来构建数据帧datwe,我将在此处添加它.谢谢.

编辑:示例 datwe

set.seed(12)
data.frame(replicate(37,sample(c(1,2,99),10,rep=TRUE)))