p4 <- ggplot(upper.ws.df, aes(elapsed.sec, WindDir)) +
geom_point(colour = "tan", size = 0.25) +
stat_smooth(aes(color = '8'), size = 0.5, method = lm,
  formula = y ~ splines::bs(x, 8), se = FALSE ) +
  # more code
puwdir.smoothed <- ggplotly(p4)

axz <- list(title = 'Direction',
      ticketmode = 'array',
      ticktext = c('E', 'NE', 'N'),
      tickvals = c(0, 2*pi/8, 2*pi/4) )
# Create the plot
p3d.scatter <- plot_ly(upper.ws.df,
        x = ~elapsed.sec, y = ~WindSpeed, z = ~WindDir,
        type = "scatter3d", mode = 'markers',
        marker = list(size = 0.6, color= ~elapsed.sec,
        colorscale = list(c(0,'#BA52ED'), c(1,'#FCB040')))
        ) %>%,

  layout(
    scene = list(
      xaxis = list(title = 'Elapsed Time (s)'),
      yaxis = list(title = 'Speed (m/s'),
      zaxis = axz)
      )

p3d.scatter