Plotly滑块初始位置与指定值不符问题（含R语言代码）

问题描述 #

在R语言Plotly可视化代码中，设置初始滑块值为600（对应t_values序列的中间位置），但滑块实际起始于t=650的位置（序列第4位），不过图表显示的是t=600对应的正确图形；标题显示为t=700的问题暂时忽略。

问题原因 #

Plotly滑块配置的active参数使用0-based索引（从0开始计数），而R语言中which()函数返回的是1-based索引（从1开始计数）。代码中直接将which(t_values == initial_slider_value)的结果赋值给active，导致索引偏移一位：

• t_values = c(500, 550, 600, 650, 700)，600的1-based索引是3，对应0-based索引应为2
• 代码中用3作为active值，会激活第4个滑块选项（650）

修复方案 #

将initial_slider_index转换为0-based索引，即在which()结果基础上减1：

initial_slider_index <- which(t_values == initial_slider_value) - 1

完整修复代码 #

前置参数与函数定义 #

alpha = 0.86 # Expenditure share on all other goods
beta = 0.6 # Cost share of housing capital in housing production function
g = 0.0005 # Scaling on housing production function
theta = 3 # Radians available for construction (benchmark is 3)
target_population <- 2000000
target_border <- 40
inc <- 70000
t <- 600
ra <- 45000
u <- 5353.65
x <- 0

get_price <- function(x,inc,t,u) {
  (((alpha^alpha) * ((1 - alpha)^(1 - alpha)) * (inc - t * x)) / u)^(1 / (1 - alpha))
}

get_quantity <- function(x,inc,t,u) {
  (((1 - alpha) * (inc - t * x)) / (((alpha^alpha) * ((1 - alpha)^(1 - alpha)) * (inc - t * x)) / u)^(1 / (1 - alpha)))
}

get_capital_intensitinc <- function(x,inc,t,u) {
  (1 / (get_price(x,inc,t,u) * (beta) * (g)))^(1 / (beta - 1))
}

get_rent <- function(x,inc,t,u) {
  (get_price(x,inc,t,u) * (g)) * ((1/(get_price(x,inc,t,u) * (beta) * (g)))^(beta / (beta - 1))) - 1 * ((1 / (get_price(x,inc,t,u) * (beta) * (g)))^(1 / (beta - 1)))
}

get_far <- function(x,inc,t,u) {
  (g) * get_capital_intensitinc(x,inc,t,u)^(beta)
}

get_density <- function(x,inc,t,u) {
  get_far(x,inc,t,u) / get_quantity(x,inc,t,u)
}

get_border <- function(target_border,inc,t,u) {
  uniroot(\(x) get_rent(x,inc,t,u) - ra, lower = 0, upper = 100)$root
}

get_population <- function(inc,t,u) {
  theta * integrate(\(x) {x * get_density(x,inc,t,u)},lower = 0, upper = get_border(target_border,inc,t,u))$value
}

get_u <- function(inc,t,target_population) {
  uniroot(\(u) get_population(inc,t,u) - target_population, c(4500, 6500))$root
}

u <- get_u(inc,t,target_population)

修复后的可视化代码 #

library(plotly)

x <- seq(0, 50, by = 0.1)
t_values <- seq(500, 700, by = 50)  # Adjusted t values

# Fixed y-axis range
y_range <- c(0, get_density(0, inc, 700, get_u(inc, 700, target_population)))  # Adjust this to your desired y-axis range

# create steps and plot all traces
steps <- list()
fig <- plot_ly()

aval <- list()
slider_steps <- list()

initial_slider_value <- 600  # Initial slider value

for (i in 1:length(t_values)) {
  t_value <- t_values[i]
  # Update u value based on the current t value
  u <- get_u(inc, t_value, target_population)

  # Update the border and population based on the new u value
  border_value <- get_border(target_border, inc, t_value, u)
  population_value <- get_population(inc, t_value, u)

  # Update the title to include the current t_value
  title_text <- paste("Transport cost:", t_value)

  aval[[i]] <- list(visible = t_value == initial_slider_value,  # Show line for the initial slider value
                    name = paste0('v = ', t_value),
                    x = x,
                    y = get_density(x, inc, t_value, u))

  fig <- add_lines(fig, x = aval[i][[1]]$x, y = aval[i][[1]]$y, visible = aval[i][[1]]$visible,
                   name = aval[i][[1]]$name, type = 'scatter', mode = 'lines', hoverinfo = 'name',
                   line = list(color = '00CED1'), showlegend = FALSE)

  step <- list(args = list('visible', rep(FALSE, length(aval))),
               method = 'restyle')
  step$args[[2]][i] = TRUE
  steps[[i]] = step

  # Define the slider step
  slider_steps[[i]] <- list(
    label = t_value,
    method = 'restyle',
    args = list('visible', rep(FALSE, length(t_values)))
  )
  slider_steps[[i]]$args[[2]][i] = TRUE
}

# 修复：转换为0-based索引
initial_slider_index <- which(t_values == initial_slider_value) - 1

# add slider control to plot
fig <- fig %>%
  layout(sliders = list(
    list(currentvalue = list(prefix = "Transport cost: "),
         steps = slider_steps,
         active = initial_slider_index  # Set the initial active slider
    )
  ),
  title = title_text,  # Set the initial title（此问题暂时忽略）
  yaxis = list(range = y_range),  # Set fixed y-axis range
  xaxis = list(range = c(0, 50))  # Set the desired x-axis range
  )

fig

内容的提问来源于stack exchange，提问作者Victor Nielsen