构建指数拟合方程初始值选择的响应式Shiny应用及图表修复

需要基于以下指数模型构建Shiny应用，通过滑块输入模型参数初始值，绘制拟合曲线与原始数据：

y = H0_1 * (1 - exp(- (Tmax / beta) ^ theta)) + c

当前代码运行后无图表显示，原始代码如下：

x <- Tmax <- 443, 454, 451, 438, 451, 452, 453, 454, 453, 445, 449, 449, 454
y <- HI <- 43, 62, 63, 95, 105, 117.51, 119.07, 122, 122, 125, 131.8, 137, 139

#save this script as app.R
library(shiny)
library(ggplot2)

ui <- fluidPage(

Application title
titlePanel("Exponential Equation App"),

Sidebar with a slider input for number of bins 
sidebarLayout(
sidebarPanel(
sliderInput(inputId = "H0_1",label = "H0_1:",
min = 100,max = 1000,value = 100, step = 0.5),
sliderInput(inputId = "beta",label = "beta",
min = 300,max = 600,value = 300, step = 0.5),
sliderInput(inputId = "theta",label = "theta",
min = -200,max = -1,value = -5, step = 1),
sliderInput(inputId = "c",label = "c",
min = 0,max = 100,value = 0, step = 1)

    ),


mainPanel(
plotOutput("lineplot")
    )
  )
)


server <- function(input, output) {

output$lineplot <- renderPlot({
x <- seq(from = 427, to = 458, by = 1)
y <-  H0_1 * (1 - exp(- (Tmax / beta) ^ theta)) + c 

plot(x,y, col="red", lwd = 3, type = "l")
lines(y~For_RStudio$Tmax, col="blue", lwd=3)
legend("topleft",c("real data","constructed"), 
col=c("blue","red"), lwd=3)


  })
}
shinyApp(ui = ui, server = server)

以下是导致图表不显示的核心问题及解决办法：

• 数据初始化错误：R中向量必须用c()包裹，原始代码直接写x <- Tmax <- 443, 454,...会引发语法错误，需改为Tmax <- c(443, 454,...)。
• UI注释未标记：UI中的Application title等注释未加#，会被当作普通文本渲染，破坏页面结构，需添加注释符号。
• 未引用滑块输入值：服务器端计算拟合曲线时，直接使用变量名H0_1、beta等，R无法识别，必须通过input$H0_1、input$beta获取滑块参数。
• 变量长度不匹配：拟合曲线的x是序列seq(427,458,1)，但计算y时用了原始数据Tmax，导致x和y长度不一致，需将Tmax替换为x。
• 内置函数名冲突：参数c是R内置函数，用作变量名会引发错误，改为c_val避免冲突。
• 无效数据引用：For_RStudio$Tmax不存在，直接使用原始的Tmax和HI数据即可；原始数据是散点，用points()代替lines()更合适。

# 初始化数据
Tmax <- c(443, 454, 451, 438, 451, 452, 453, 454, 453, 445, 449, 449, 454)
HI <- c(43, 62, 63, 95, 105, 117.51, 119.07, 122, 122, 125, 131.8, 137, 139)

# 保存为app.R
library(shiny)

ui <- fluidPage(
  # Application title
  titlePanel("Exponential Equation App"),
  
  # Sidebar with parameter sliders
  sidebarLayout(
    sidebarPanel(
      sliderInput(inputId = "H0_1", label = "H0_1:",
                  min = 100, max = 1000, value = 150, step = 0.5),
      sliderInput(inputId = "beta", label = "beta:",
                  min = 300, max = 600, value = 450, step = 0.5),
      sliderInput(inputId = "theta", label = "theta:",
                  min = -20, max = -1, value = -5, step = 1), # 缩小范围避免数值溢出
      sliderInput(inputId = "c_val", label = "c:",
                  min = 0, max = 100, value = 0, step = 1)
    ),
    
    mainPanel(
      plotOutput("lineplot")
    )
  )
)

server <- function(input, output) {
  output$lineplot <- renderPlot({
    # 生成拟合曲线的x序列
    x_seq <- seq(from = 427, to = 458, by = 1)
    # 根据滑块参数计算拟合y值
    y_fit <- input$H0_1 * (1 - exp(- (x_seq / input$beta) ^ input$theta)) + input$c_val
    
    # 绘制图表
    plot(x_seq, y_fit, col = "red", lwd = 3, type = "l", 
         xlab = "Tmax", ylab = "HI", main = "Exponential Fit vs Raw Data")
    # 绘制原始数据散点
    points(Tmax, HI, col = "blue", pch = 16, cex = 1.2)
    # 添加图例
    legend("topleft", c("拟合曲线", "原始数据"), 
           col = c("red", "blue"), lwd = 3, pch = c(NA, 16), cex = 1)
  })
}

shinyApp(ui = ui, server = server)

• 调整了theta的滑块范围（-20到-1），避免过大的负数导致数值计算溢出。
• 用points()展示原始数据散点，比lines()更直观反映数据分布。
• 添加了坐标轴标签和标题，提升图表可读性。

内容的提问来源于stack exchange，提问作者Khal Aboudi