Leer datos, cargar paquetes y grafico de datos crudos

Sin estandarizar

library(readr) # lex excel
library(dplyr) # data organization
library(GGally) # grafico coordenadas paralelas
library(ggplot2) # ggplot
library(plotly) # exporta a grafico interactivo


Datos_anali <- read_csv("C:/Users/diego.lizcano/Box Sync/CodigoR/ANLA/data/Analisis_integral_ARC11_Prueba_Graph.csv")
head(Datos_anali)
## # A tibble: 6 x 9
##   Narrativa Potencia_generacion_MW Longitud_rios  DORw Retencion_sediment~
##       <int>                  <dbl>         <dbl> <dbl>               <dbl>
## 1         1                 12823          8348.  5.52                38.5
## 2         1                  9737          8383.  3.02                37.7
## 3         1                  8187          8503.  3.76                37.5
## 4         1                 11511          8320.  5.21                40.0
## 5         1                  9074.         8317.  3.3                 39.2
## 6         1                 19841          7828. 13.0                 48.4
## # ... with 4 more variables: Bosque_seco <dbl>, Reservas_ley2 <dbl>,
## #   Riqueza_esp_amenazadas <dbl>, Riqueza_mamiferos_alto <dbl>
# st_dat <-  filter(Datos_anali, Narrativa=="1", Potencia_generacion_MW < r1)
# for(i)
# r1 <- (((max(st_dat[,2])-min(st_dat[,2]))/10) *1) + min(st_dat[,2])
# length(which(st_dat[,2] < r1))

#1er grafico intuiutivo?
ggparcoord(Datos_anali, columns = 2:9,  scale = 'globalminmax')

Armar graficos graficos ggplot

Datos_anali_scaled <- scale(Datos_anali[,2:9], center = TRUE)
Datos_anali$Narrativa <- as.factor(Datos_anali$Narrativa)

p1 <- ggparcoord(Datos_anali, columns = 2:9,  
           groupColumn = 1, 
           scale = "uniminmax",
           order = "anyClass",
           alphaLines = 0.2) + 
            theme(axis.text.x = element_text(angle = 270, hjust = 0)) +
            scale_colour_brewer(palette="Set1")

# p1



p2 <- ggparcoord(Datos_anali, columns = 2:9,  
           groupColumn = 1, 
           scale = "uniminmax",
           order = "anyClass",
           alphaLines = 0.1) + 
            facet_grid(Narrativa ~ .) +
            theme(axis.text.x = element_text(angle = 270, hjust = 0)) +
            scale_colour_brewer(palette="Set1")


# p2

Ver graficos ggplot (no interactivos)

p1

p2

Including plotly Plots (graficos interactivos)

first plot all in the same

Click on narrative leyend to on or off.

gl <- ggplotly(p1)
gl # see
# chart_link = api_create(pl, filename="ggparcoord_carlos") # too long

Wrap option by Narrativa

g2 <- ggplotly(p2)
g2 # see

Pure Plotly by Carlos R.

# MyData <- read.csv(file="C:/Users/carlos.rogeliz/Documents/2_TNC_Proyectos/5_Convenio TNC-ANLA_046_2016/10_R_Parallel_Axis/Analisis_integral_ARC11_Prueba_Graph.csv", header=TRUE, sep=",")

MyData <- Datos_anali

p <- MyData %>%
  plot_ly(width = 1000, height = 600) %>%
  add_trace(type = 'parcoords',
            line = list(color = ~Narrativa,
                        colorscale = 'Jet',
                        showscale = TRUE,
                        reversescale = TRUE,
                        cmin = 1,
                        cmax = 7),
            dimensions = list(
              list(range = c(~min(Potencia_generacion_MW),~max(Potencia_generacion_MW)),
                   label = 'Potencia Generacion (MW)', values = ~Potencia_generacion_MW),
              
              list(range = c(~min(Longitud_rios),9800),#~max(Longitud_rios)),
                   label = 'Longitud Rios Libres (Km)', values = ~Longitud_rios),
              
              list(range = c(~min(DORw),~max(DORw)),
                   label = 'DORw (%)', values = ~DORw),
              
              list(range = c(~min(Retencion_sedimentos),~max(Retencion_sedimentos)),
                   label = 'Alteracion Sedimentos (%)', values = ~Retencion_sedimentos),
              
              list(range = c(~min(Bosque_seco),~max(Bosque_seco)),
                   label = 'Bosque Seco (Ha)', values = ~Bosque_seco),
              
              list(range = c(~min(Reservas_ley2),~max(Reservas_ley2)),
                   label = 'Reservas Ley2 (Ha)', values = ~Reservas_ley2),
              
              list(range = c(~min(Riqueza_esp_amenazadas),~max(Riqueza_esp_amenazadas)),
                   label = 'Riqueza Spp Amenazadas', values = ~Riqueza_esp_amenazadas)
            )
  )


### ver grafica
p
# Link para publicacion en https://plot.ly/~dlizcano/32/#/
# chart_link = api_create(p, filename="parcoords-advanced")
# chart_link