R 如何根据用户在屏幕中的输入显示不同的绘图?
我有这个数据框:R 如何根据用户在屏幕中的输入显示不同的绘图?,r,plot,shiny,user-input,selectinput,R,Plot,Shiny,User Input,Selectinput,我有这个数据框: > df genes enst x y 1 Gene1 ENST1 25 14 2 Gene1 ENST2 60 25 3 Gene1 ENST3 12 5 4 Gene2 ENST1 9 34 5 Gene2 ENST2 14 12 6 Gene3 ENST1 10 1 我正在尝试创建一个闪亮的应用程序,它允许我选择基因和转录本。如果您选择一个基因,例如Gene1,您将有一个选项来选择在本例中您想要哪个转录本,ENST1,ENST2,ENST3。 问
> df
genes enst x y
1 Gene1 ENST1 25 14
2 Gene1 ENST2 60 25
3 Gene1 ENST3 12 5
4 Gene2 ENST1 9 34
5 Gene2 ENST2 14 12
6 Gene3 ENST1 10 1
我正在尝试创建一个闪亮的应用程序,它允许我选择基因和转录本。如果您选择一个基因,例如Gene1,您将有一个选项来选择在本例中您想要哪个转录本,ENST1,ENST2,ENST3。
问题是我想画两幅图。如果单击“基因级别:基因”,它将对该基因的所有值求和。
例如,对于具有3个转录本的第一个基因,x的总值将为20+60+12=92,y的总值将为14+25+5=44。因此,绘制基因1的值为:x=92,y=44
此外,我想画出每个成绩单。例如,如果选择Gene1和转录本1,绘图将使用x=25和y=14。但是,如果用户决定选择两个转录本,则用户将看到两个绘图。或者,如果用户选择3个转录本,用户将看到3个不同的绘图
现在,用我的代码:
如果你选择了这个基因,你会得到这个基因的图谱。
但是,它将显示同一情节中的所有转录本。如果用户需要,我只想显示一个或多个成绩单
我不知道如何继续
另一方面,有两件事我不知道如何实现
一个动作按钮可以看到绘图。-如果你点击它,它会显示绘图。
单选按钮级别如果你点击基因,它只会显示基因图。然而,如果你点击成绩单,它会显示成绩单图。
有人能帮我吗?提前谢谢
我的代码:
library(shiny)
################ DATA #############################
genes<- c("Gene1", "Gene1", "Gene1", "Gene2", "Gene2", "Gene3")
enst <- c("ENST1", "ENST2", "ENST3", "ENST1", "ENST2", "ENST1")
x <- c(25, 60, 12, 9, 14, 10)
y <- c(14, 25, 5, 34, 12, 1)
df<- data.frame(genes, enst, x, y)
###################################################
ui <- fluidPage(
# Application title
titlePanel("Barplot"),
sidebarLayout(
sidebarPanel(
uiOutput("selected_gene"),
uiOutput("selected_transcript"),
radioButtons("level", "Level:",
c("Gene" = "Gene",
"Transcript" = "Transcript")),
h5(strong("If you want to see the plot, you have to click the button")),
actionButton("add_plot", "See the plot"),
),
mainPanel(
plotOutput("plot"),
plotOutput("plot2"),
tableOutput("table1"),
tableOutput("table2")
)
)
)
server <- function(session, input, output) {
# This function gives us the list of genes.
genes_list <- reactive({
df$genes
})
transcripts_list <- reactive({
transcripts <- subset(df, df$genes==input$gene)
transcripts <- transcripts[,2]
return(transcripts)
})
# This function give us a select list input, in order to be able to select the gene that we want to see
output$selected_gene <- renderUI({
selectizeInput(inputId = "gene", "Select one gene", choices=genes_list(), options=list(maxOptions = length(genes_list())))
})
output$selected_transcript <- renderUI({
selectizeInput(inputId = "transcript", "Select one transcript", choices=transcripts_list(), options=list(maxOptions = length(transcripts_list())), multiple=T)
})
gene_values <- reactive({
values <- subset(df, df[1]==input$gene)
values$enst <- NULL
if(nrow(values)>1){ #for those genes who have more than 1 transcript
values_new <- values[2:length(values)]
values_new <- as.data.frame(t(colSums(values_new))) # sum the columns, transpose and transform into a dataframe
gene <- values[1,] #we take the first row, only one gene but all the info.
values <- cbind(values_new, gene[1]) # we bind both dataframes, however, we only want the gene name
values <- values[,c("genes",setdiff(names(values),"genes"))] # we move the last column at the beginning
}
return(values)
})
transc_values <- reactive({
values <- subset(df, df[1]==input$gene)
values$genes <- NULL
return(values)
})
plot_genes <- reactive({
gene_values <- gene_values()
barplot(c(gene_values$x, gene_values$y))
})
plot_transc <- reactive({
transc_values <- transc_values()
barplot(c(transc_values$x, transc_values$y))
})
v <- reactiveValues(plot = NULL)
observeEvent(input$add_plot, {
if(input$level == "Gene"){
v$plot <- plot_genes()
}
if(input$level == "Transcript"){
v$plot <- plot_transc()
}
})
# This function will draw the plot
# output$plot <- renderPlot({
# if (is.null(v$plot)){
# return()
# }
# v$plot
# })
output$table1 <- renderTable(gene_values())
output$table2 <- renderTable(transc_values())
output$plot <- renderPlot(plot_genes())
output$plot2 <- renderPlot(plot_transc())
}
shinyApp(ui, server)
也许您可以从这个开始,并根据您的需要进行修改
library(shiny)
library(ggplot2)
library(DT)
################ DATA #############################
genes<- c("Gene1", "Gene1", "Gene1", "Gene2", "Gene2", "Gene3")
enst <- c("ENST1", "ENST2", "ENST3", "ENST1", "ENST2", "ENST1")
x <- c(25, 60, 12, 9, 14, 10)
y <- c(14, 25, 5, 34, 12, 1)
df<- data.frame(genes, enst, x, y)
###################################################
ui <- fluidPage(
# Application title
titlePanel("Histogram"),
sidebarLayout(
sidebarPanel(
uiOutput("selected_gene"),
uiOutput("selected_transcript"),
radioButtons("level", "Level:",
c("Gene" = "Gene",
"Transcript" = "Transcript")),
h5(strong("If you want to see the plot, you have to click the button")),
div(actionButton("add_plot", "See the plot"),
actionButton("table", "See the table"),
actionButton("clear", "Clear All")
)
),
mainPanel(
plotOutput("plot"),
DTOutput("table")
)
)
)
server <- function(input, output, session) {
## This function gives us the list of genes.
genes_list <- reactive({
unique(df$genes)
})
transcripts_list <- reactive({
req(input$gene)
transcripts <- subset(df, df$genes==input$gene)
transcripts <- transcripts[,2]
return(unique(transcripts))
})
# This function give us a select list input, in order to be able to select the gene that we want to see
output$selected_gene <- renderUI({
selectizeInput(inputId = "gene", "Select one gene", choices=genes_list(), options=list(maxOptions = length(genes_list())))
})
output$selected_transcript <- renderUI({
selectizeInput(inputId = "transcript", "Select one transcript", choices=transcripts_list(), options=list(maxOptions = length(transcripts_list())), multiple=F)
})
gene_values <- reactive({
req(input$gene)
values <- subset(df, df[1]==input$gene)
values$enst <- NULL
if(nrow(values)>1){ #for those genes who have more than 1 transcript
values_new <- values[2:length(values)]
values_new <- as.data.frame(t(colSums(values_new))) # sum the columns, transpose and transform into a dataframe
gene <- values[1,] #we take the first row, only one gene but all the info.
values <- cbind(values_new, gene[1]) # we bind both dataframes, however, we only want the gene name
values <- values[,c("genes",setdiff(names(values),"genes"))] # we move the last column at the beginning
}
return(values)
})
transc_values <- reactive({
req(input$transcript)
values <- subset(df, df[2]==input$transcript)
values$genes <- NULL
return(values)
})
mydata <- reactive({
req(input$level)
if(input$level == "Gene"){
df <- req(gene_values())
}else if(input$level == "Transcript"){
df <- req(transc_values())
}else df <- NULL
df
})
# plot_genes <- reactive({
# gene_values <- req(gene_values())
# barplot(c(gene_values$x, gene_values$y))
#
# })
#
# plot_transc <- reactive({
# transc_values <- req(transc_values())
# barplot(c(transc_values$x, transc_values$y))
#
# })
v <- reactiveValues(plot = NULL, table=NULL)
observeEvent(input$add_plot, {
v$plot <- ggplot(mydata(), aes(x=x,y=y)) + geom_bar(stat = "identity")
v$table <- NULL ### display only plot
},ignoreInit = TRUE)
observeEvent(input$table, {
v$table <- req(mydata())
v$plot <- NULL ### display only table
},ignoreInit = TRUE)
observeEvent(input$clear, {
v$table <- NULL
v$plot <- NULL
},ignoreInit = TRUE)
## This function will draw the plot
output$plot <- renderPlot({ v$plot })
output$table <- renderDT({ v$table })
}
shinyApp(ui, server)
令人印象深刻,非常感谢!可以同时显示两个图吗?例如,如果一个基因有两个转录本,用户选择两个转录本,他能同时看到这两个图吗?或者,如果一个基因有3个转录本,用户只想看到2个,那么他只能看到两个图吗?是的。您可以显示您想要的任何内容-在observeEventinput$add_plot,{…}中定义它,然后创建一个新的输出$plot2…但是。。。如果您不知道用户需要多少个绘图?因为根据您所说的,您需要知道用户想要多少个绘图,以便在代码中写入输出$plotX