R plotly循环add_trace添加多轨迹仅显示最后一条问题咨询

在R语言环境使用plotly包批量绘制多轨迹图表时，通过for循环调用add_trace()添加轨迹，最终图表仅展示最后一次循环添加的轨迹；手动通过管道符逐次链式调用add_trace()添加相同轨迹时，可正常展示全部内容。

• 复现异常的循环代码：

Fig3abis<-plot_ly(data=dataresults3ab,x=dataresults3ab[[1]])
    
for(j in c(100,500,1000,1500))
{
    Fig3abis <- add_trace(Fig3abis,
                   y=~dataresults3ab[,j],
                   name=paste("N1",as.character(j),sep = "_"),
                   type="scatter",
                   mode="markers",
                   marker=list(size=4,color="black"))
}

Fig3abis <- Fig3abis%>% layout(title="Bifurcation diagram for five species competing for five resources. Local minima and maxima of species 1, 
                     from t=2,000 to t=4,000 days, as a function of the half-saturation constant K41",
                         showlegend=F,
                         xaxis=list(title="Half-saturation constant K41, of species 1",range=c(0.1,0.5)),
                         yaxis=list(title="Abundancie of species 1",range=c(0,100))) 

• 可正常运行的手动链式调用代码：

Fig3a <- plot_ly(data=dataresults3ab,x=dataresults3ab[[1]])

Fig3a <- add_trace(Fig3a,
                      y=~dataresults3ab[,100],
                      name="N1_2",
                      type="scatter",
                      mode="markers",
                      marker=list(size=4,color="black"))%>%
       add_trace(y=~dataresults3ab[,500],
                   name="N_3",
                   type="scatter",
                   mode="markers",
                   marker=list(size=4,color="black"))%>%
       add_trace(y=~dataresults3ab[,1000],
                   name="N1_4",
                   type="scatter",
                   mode="markers",
                   marker=list(size=4,color="black"))%>%
        add_trace(y=~dataresults3ab[,1500],
                   name="N1_5",
                   type="scatter",
                   mode="markers",
                   marker=list(size=4,color="black"))

Fig3a <- Fig3a%>% layout(title="Bifurcation diagram for five species competing for five resources. Local minima and maxima of species 1, from t=2,000 to t=4,000 days, as a function of the half-saturation constant K41",
       showlegend=F,
       xaxis=list(title="Half-saturation constant K41, of species 1",range=c(0.1,0.5)),
       yaxis=list(title="Abundancie of species 1",range=c(0,100)))

实际业务场景需批量绘制约8000条轨迹，数据源为标准data.frame结构：第一列为X轴对应取值，其余每一列分别对应一条待绘制轨迹的Y轴取值。

异常核心是plotly的公式接口（参数前加~的写法）采用延迟求值机制：

• 循环中写y=~dataresults3ab[,j]时，代码不会在当前循环步立即解析j的取值并提取对应列，只会保存这个表达式；等图表最终渲染时才会执行表达式取值，此时for循环已经运行结束，所有trace的表达式中引用的j都指向循环结束后的最终值（即测试用例中的1500），因此所有轨迹都绘制为最后一列的数据，视觉上仅展示最后一条轨迹。
• 手动链式调用时传入的是固定常量列索引（100/500/1000/1500），不存在循环变量动态变化的问题，因此渲染时可正常取到对应列的数据。

方案1（强烈推荐，适配8000条轨迹的大数量场景） #

无需循环逐次添加trace，将宽格式数据转换为长格式后，通过分组参数一次性生成所有轨迹，渲染效率远高于逐次添加trace，且完全规避延迟求值bug。

library(plotly)
library(tidyr)
library(dplyr)

# 宽表转长表：保留第一列为X轴字段，其余列转换为轨迹标识、Y值两个字段
long_data <- dataresults3ab %>%
  pivot_longer(
    cols = -1, # 排除第一列X轴字段
    names_to = "trace_id",
    values_to = "y_value"
  )

# 一次性绘图，按trace_id自动拆分独立轨迹
Fig3abis <- plot_ly(
  data = long_data,
  x = ~.[[1]], # 取第一列为X值
  y = ~y_value,
  split = ~trace_id, # 按轨迹id分组生成独立轨迹
  type = "scatter",
  mode = "markers",
  marker = list(size = 4, color = "black"),
  showlegend = F
) %>%
  layout(
    title="Bifurcation diagram for five species competing for five resources. Local minima and maxima of species 1, from t=2,000 to t=4,000 days, as a function of the half-saturation constant K41",
    xaxis=list(title="Half-saturation constant K41, of species 1",range=c(0.1,0.5)),
    yaxis=list(title="Abundancie of species 1",range=c(0,100))
  )

方案2（适配必须使用for循环的场景） #

循环中避免对循环变量使用~延迟引用，在当前循环步强制提取Y值向量，直接传入已计算完成的数值向量，阻断延迟求值：

Fig3abis<-plot_ly(data=dataresults3ab,x=dataresults3ab[[1]])
    
for(j in c(100,500,1000,1500))
{
    # 当前循环步立即提取对应列的Y值，完成求值
    current_y <- dataresults3ab[[j]]
    Fig3abis <- add_trace(Fig3abis,
                   y = current_y, # 直接传值，不使用~前缀
                   name = paste0("N1_",j),
                   type = "scatter",
                   mode = "markers",
                   marker = list(size=4,color="black"))
}

# layout配置与原代码一致

注意：该方案逐次添加8000条trace时渲染性能较差，大数量场景优先选择方案1。

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