analisis DRISs

analisis DRIS,realizado en base al curso de Nutricon vegetal Universidad de Caldas 2016.

datos iniciales de analisis foliares del cultivo de papa en fincas del Departamento de Caldas.

primero,se sube el archivo de analisis foliares que llamaremos Relaciones_binarias,el cual se ha subido en formato “data.frame” con Rstudio

library(readxl)
Relaciones_binarias <- read_excel("~/Descargas/1.xlsx")
View(Relaciones_binarias)

segundo,se calculan las variables para el analisis DRIS,por el metodo en R de agregar columnas al “data.frame” relaciones con Nitrogeno

#relacion nitrogeno fosforo
Relaciones_binarias$`N/P`<-Relaciones_binarias$`% N`/Relaciones_binarias$`% P`

Relaciones_binarias$`P/N`<-Relaciones_binarias$`% P`/Relaciones_binarias$`% N`

#relacion nitrogeno potasio

Relaciones_binarias$`N/K`<-Relaciones_binarias$`% N`/Relaciones_binarias$`% K`
Relaciones_binarias$`K/N`<-Relaciones_binarias$`% K`/Relaciones_binarias$`% N`

#relacion nitrogeno Calcio

Relaciones_binarias$`N/Ca`<-Relaciones_binarias$`% N`/Relaciones_binarias$`% Ca`
Relaciones_binarias$`Ca/N`<-Relaciones_binarias$`% Ca`/Relaciones_binarias$`% N`

#relacion nitrogeno Magnesio
Relaciones_binarias$`N/Mg`<-Relaciones_binarias$`% N`/Relaciones_binarias$`% Mg`
Relaciones_binarias$`Mg/N`<-Relaciones_binarias$`% Mg`/Relaciones_binarias$`% N`

#relacion nitrogeno Azufre
Relaciones_binarias$`N/S`<-Relaciones_binarias$`% N`/Relaciones_binarias$`% S`
Relaciones_binarias$`S/N`<-Relaciones_binarias$`% S`/Relaciones_binarias$`% N`

#para aquellos elementos que se presentan como partes por millon se debe multiplicar el elemento presentado como %, por 100.

#relacion nitrogeno Hierro
Relaciones_binarias$`N/Fe`<-((Relaciones_binarias$`% N`)*100)/Relaciones_binarias$`Fe ppm`
Relaciones_binarias$`Fe/N`<-Relaciones_binarias$`Fe ppm`/((Relaciones_binarias$`% N`)*100)

#relacion nitrogeno Boro

Relaciones_binarias$`N/B`<-((Relaciones_binarias$`% N`)*100)/Relaciones_binarias$`B ppm`
Relaciones_binarias$`B/N`<-Relaciones_binarias$`B ppm`/((Relaciones_binarias$`% N`)*100)

#relacion nitrogeno Zinc

Relaciones_binarias$`N/Zn`<-((Relaciones_binarias$`% N`)*100)/Relaciones_binarias$`Zn ppm`
Relaciones_binarias$`Zn/N`<-Relaciones_binarias$`Zn ppm`/((Relaciones_binarias$`% N`)*100)

#relacion nitrogeno Manganeso

Relaciones_binarias$`N/Mn`<-((Relaciones_binarias$`% N`)*100)/Relaciones_binarias$`Mn ppm`
Relaciones_binarias$`Mn/N`<-Relaciones_binarias$`Mn ppm`/((Relaciones_binarias$`% N`)*100)

#relacion nitrogeno Cobre

Relaciones_binarias$`N/Cu`<-((Relaciones_binarias$`% N`)*100)/Relaciones_binarias$`Cu ppm`
Relaciones_binarias$`Cu/N`<-Relaciones_binarias$`Cu ppm`/((Relaciones_binarias$`% N`)*100)

relaciones con Fosforo. como ya se han calculado las relaciones con Nitrogeno estas se han omitido.

#relacion fosforo potasio

Relaciones_binarias$`P/K`<-Relaciones_binarias$`% P`/Relaciones_binarias$`% K`
Relaciones_binarias$`K/P`<-Relaciones_binarias$`% K`/Relaciones_binarias$`% P`

#relacion fosforo Calcio

Relaciones_binarias$`P/Ca`<-Relaciones_binarias$`% P`/Relaciones_binarias$`% Ca`
Relaciones_binarias$`Ca/P`<-Relaciones_binarias$`% Ca`/Relaciones_binarias$`% P`

#relacion fosforo Magnesio
Relaciones_binarias$`P/Mg`<-Relaciones_binarias$`% P`/Relaciones_binarias$`% Mg`
Relaciones_binarias$`Mg/P`<-Relaciones_binarias$`% Mg`/Relaciones_binarias$`% P`

#relacion fosforo Azufre
Relaciones_binarias$`P/S`<-Relaciones_binarias$`% P`/Relaciones_binarias$`% S`
Relaciones_binarias$`S/P`<-Relaciones_binarias$`% S`/Relaciones_binarias$`% P`

#para aquellos elementos que se presentan como partes por millon se debe multiplicar el elemento presentado como %, por 100.

#relacion fosforo Hierro
Relaciones_binarias$`P/Fe`<-((Relaciones_binarias$`% P`)*100)/Relaciones_binarias$`Fe ppm`
Relaciones_binarias$`Fe/P`<-Relaciones_binarias$`Fe ppm`/((Relaciones_binarias$`% P`)*100)

#relacion fosforo Boro

Relaciones_binarias$`P/B`<-((Relaciones_binarias$`% P`)*100)/Relaciones_binarias$`B ppm`
Relaciones_binarias$`B/P`<-Relaciones_binarias$`B ppm`/((Relaciones_binarias$`% P`)*100)

#relacion fosforo Zinc

Relaciones_binarias$`P/Zn`<-((Relaciones_binarias$`% P`)*100)/Relaciones_binarias$`Zn ppm`
Relaciones_binarias$`Zn/P`<-Relaciones_binarias$`Zn ppm`/((Relaciones_binarias$`% P`)*100)

#relacion fosforo Manganeso

Relaciones_binarias$`P/Mn`<-((Relaciones_binarias$`% P`)*100)/Relaciones_binarias$`Mn ppm`
Relaciones_binarias$`Mn/P`<-Relaciones_binarias$`Mn ppm`/((Relaciones_binarias$`% P`)*100)

#relacion fosforo Cobre

Relaciones_binarias$`P/Cu`<-((Relaciones_binarias$`% P`)*100)/Relaciones_binarias$`Cu ppm`
Relaciones_binarias$`Cu/P`<-Relaciones_binarias$`Cu ppm`/((Relaciones_binarias$`% P`)*100)

relaciones con Potasio. como ya se han calculado las relaciones con Nitrogeno y fosforo estas se han omitido.

#relacion potasio Calcio

Relaciones_binarias$`K/Ca`<-Relaciones_binarias$`% K`/Relaciones_binarias$`% Ca`
Relaciones_binarias$`Ca/K`<-Relaciones_binarias$`% Ca`/Relaciones_binarias$`% K`

#relacion potasio Magnesio
Relaciones_binarias$`K/Mg`<-Relaciones_binarias$`% K`/Relaciones_binarias$`% Mg`
Relaciones_binarias$`Mg/K`<-Relaciones_binarias$`% Mg`/Relaciones_binarias$`% K`

#relacion potasio Azufre
Relaciones_binarias$`K/S`<-Relaciones_binarias$`% K`/Relaciones_binarias$`% S`
Relaciones_binarias$`S/K`<-Relaciones_binarias$`% S`/Relaciones_binarias$`% K`

#para aquellos elementos que se presentan como partes por millon se debe multiplicar el elemento presentado como %, por 100.

#relacion potasio Hierro
Relaciones_binarias$`K/Fe`<-((Relaciones_binarias$`% K`)*100)/Relaciones_binarias$`Fe ppm`
Relaciones_binarias$`Fe/K`<-Relaciones_binarias$`Fe ppm`/((Relaciones_binarias$`% K`)*100)

#relacion potasio Boro

Relaciones_binarias$`K/B`<-((Relaciones_binarias$`% K`)*100)/Relaciones_binarias$`B ppm`
Relaciones_binarias$`B/K`<-Relaciones_binarias$`B ppm`/((Relaciones_binarias$`% K`)*100)

#relacion potasio Zinc

Relaciones_binarias$`K/Zn`<-((Relaciones_binarias$`% K`)*100)/Relaciones_binarias$`Zn ppm`
Relaciones_binarias$`Zn/K`<-Relaciones_binarias$`Zn ppm`/((Relaciones_binarias$`% K`)*100)

#relacion potasio Manganeso

Relaciones_binarias$`K/Mn`<-((Relaciones_binarias$`% K`)*100)/Relaciones_binarias$`Mn ppm`
Relaciones_binarias$`Mn/K`<-Relaciones_binarias$`Mn ppm`/((Relaciones_binarias$`% K`)*100)

#relacion potasio Cobre

Relaciones_binarias$`K/Cu`<-((Relaciones_binarias$`% K`)*100)/Relaciones_binarias$`Cu ppm`
Relaciones_binarias$`Cu/K`<-Relaciones_binarias$`Cu ppm`/((Relaciones_binarias$`% K`)*100)

relaciones con Calcio. como ya se han calculado las relaciones con Nitrogeno, Fosforo y Potasio estas se han omitido.

#relacion Calcio Magnesio
Relaciones_binarias$`Ca/Mg`<-Relaciones_binarias$`% Ca`/Relaciones_binarias$`% Mg`
Relaciones_binarias$`Mg/Ca`<-Relaciones_binarias$`% Mg`/Relaciones_binarias$`% Ca`

#relacion Calcio Azufre
Relaciones_binarias$`Ca/S`<-Relaciones_binarias$`% Ca`/Relaciones_binarias$`% S`
Relaciones_binarias$`S/Ca`<-Relaciones_binarias$`% S`/Relaciones_binarias$`% Ca`

#para aquellos elementos que se presentan como partes por millon se debe multiplicar el elemento presentado como %, por 100.

#relacion Calcio Hierro
Relaciones_binarias$`Ca/Fe`<-((Relaciones_binarias$`% Ca`)*100)/Relaciones_binarias$`Fe ppm`
Relaciones_binarias$`Fe/Ca`<-Relaciones_binarias$`Fe ppm`/((Relaciones_binarias$`% Ca`)*100)

#relacion Calcio Boro

Relaciones_binarias$`Ca/B`<-((Relaciones_binarias$`% Ca`)*100)/Relaciones_binarias$`B ppm`
Relaciones_binarias$`B/Ca`<-Relaciones_binarias$`B ppm`/((Relaciones_binarias$`% Ca`)*100)

#relacion Calcio Zinc

Relaciones_binarias$`Ca/Zn`<-((Relaciones_binarias$`% Ca`)*100)/Relaciones_binarias$`Zn ppm`
Relaciones_binarias$`Zn/Ca`<-Relaciones_binarias$`Zn ppm`/((Relaciones_binarias$`% Ca`)*100)

#relacion Calcio Manganeso

Relaciones_binarias$`Ca/Mn`<-((Relaciones_binarias$`% Ca`)*100)/Relaciones_binarias$`Mn ppm`
Relaciones_binarias$`Mn/Ca`<-Relaciones_binarias$`Mn ppm`/((Relaciones_binarias$`% Ca`)*100)

#relacion Calcio Cobre

Relaciones_binarias$`Ca/Cu`<-((Relaciones_binarias$`% Ca`)*100)/Relaciones_binarias$`Cu ppm`
Relaciones_binarias$`Cu/Ca`<-Relaciones_binarias$`Cu ppm`/((Relaciones_binarias$`% Ca`)*100)

relaciones con Magnesio. como ya se han calculado las relaciones con Nitrogeno, Fosforo,Potasio y Calcio estas se han omitido.

#relacion Magnesio Azufre
Relaciones_binarias$`Mg/S`<-Relaciones_binarias$`% Mg`/Relaciones_binarias$`% S`
Relaciones_binarias$`S/Mg`<-Relaciones_binarias$`% S`/Relaciones_binarias$`% Mg`

#para aquellos elementos que se presentan como partes por millon se debe multiplicar el elemento presentado como %, por 100.

#relacion Magnesio Hierro
Relaciones_binarias$`Mg/Fe`<-((Relaciones_binarias$`% Mg`)*100)/Relaciones_binarias$`Fe ppm`
Relaciones_binarias$`Fe/Mg`<-Relaciones_binarias$`Fe ppm`/((Relaciones_binarias$`% Mg`)*100)

#relacion Magnesio Boro

Relaciones_binarias$`Mg/B`<-((Relaciones_binarias$`% Mg`)*100)/Relaciones_binarias$`B ppm`
Relaciones_binarias$`B/Mg`<-Relaciones_binarias$`B ppm`/((Relaciones_binarias$`% Mg`)*100)

#relacion Magnesio Zinc

Relaciones_binarias$`Mg/Zn`<-((Relaciones_binarias$`% Mg`)*100)/Relaciones_binarias$`Zn ppm`
Relaciones_binarias$`Zn/Mg`<-Relaciones_binarias$`Zn ppm`/((Relaciones_binarias$`% Mg`)*100)

#relacion Magnesio Manganeso

Relaciones_binarias$`Mg/Mn`<-((Relaciones_binarias$`% Mg`)*100)/Relaciones_binarias$`Mn ppm`
Relaciones_binarias$`Mn/Mg`<-Relaciones_binarias$`Mn ppm`/((Relaciones_binarias$`% Mg`)*100)

#relacion Magnesio Cobre

Relaciones_binarias$`Mg/Cu`<-((Relaciones_binarias$`% Mg`)*100)/Relaciones_binarias$`Cu ppm`
Relaciones_binarias$`Cu/Mg`<-Relaciones_binarias$`Cu ppm`/((Relaciones_binarias$`% Mg`)*100)

relaciones con Azufre. como ya se han calculado las relaciones con Nitrogeno, Fosforo,Potasio,Calcio y Magnesio estas se han omitido.

#para aquellos elementos que se presentan como partes por millon se debe multiplicar el elemento presentado como %, por 100.

#relacion Azufre Hierro
Relaciones_binarias$`S/Fe`<-((Relaciones_binarias$`% S`)*100)/Relaciones_binarias$`Fe ppm`
Relaciones_binarias$`Fe/S`<-Relaciones_binarias$`Fe ppm`/((Relaciones_binarias$`% S`)*100)

#relacion Azufre Boro

Relaciones_binarias$`S/B`<-((Relaciones_binarias$`% S`)*100)/Relaciones_binarias$`B ppm`
Relaciones_binarias$`B/S`<-Relaciones_binarias$`B ppm`/((Relaciones_binarias$`% S`)*100)

#relacion Azufre Zinc

Relaciones_binarias$`S/Zn`<-((Relaciones_binarias$`% S`)*100)/Relaciones_binarias$`Zn ppm`
Relaciones_binarias$`Zn/S`<-Relaciones_binarias$`Zn ppm`/((Relaciones_binarias$`% S`)*100)

#relacion Azufre Manganeso

Relaciones_binarias$`S/Mn`<-((Relaciones_binarias$`% S`)*100)/Relaciones_binarias$`Mn ppm`
Relaciones_binarias$`Mn/S`<-Relaciones_binarias$`Mn ppm`/((Relaciones_binarias$`% S`)*100)

#relacion Azufre Cobre

Relaciones_binarias$`S/Cu`<-((Relaciones_binarias$`% S`)*100)/Relaciones_binarias$`Cu ppm`
Relaciones_binarias$`Cu/S`<-Relaciones_binarias$`Cu ppm`/((Relaciones_binarias$`% S`)*100)

relaciones con Hierro. como Hierro es presentado en el analisi de suelos como ppm,no es necesario multiplicar por 100 para las variables que se presentan como ppm. como ya se han calculado las relaciones con Nitrogeno, Fosforo,Potasio,Calcio,Magnesio y Azufre estas se han omitido.

#relacion Hierro Boro

Relaciones_binarias$`Fe ppm/B`<-(Relaciones_binarias$`Fe ppm`)/Relaciones_binarias$`B ppm`
Relaciones_binarias$`B/Fe ppm`<-Relaciones_binarias$`B ppm`/(Relaciones_binarias$`Fe ppm`)

#relacion Hierro Zinc

Relaciones_binarias$`Fe ppm/Zn`<-(Relaciones_binarias$`Fe ppm`)/Relaciones_binarias$`Zn ppm`
Relaciones_binarias$`Zn/Fe ppm`<-Relaciones_binarias$`Zn ppm`/(Relaciones_binarias$`Fe ppm`)

#relacion Hierro Manganeso

Relaciones_binarias$`Fe ppm/Mn`<-(Relaciones_binarias$`Fe ppm`)/Relaciones_binarias$`Mn ppm`
Relaciones_binarias$`Mn/Fe ppm`<-Relaciones_binarias$`Mn ppm`/(Relaciones_binarias$`Fe ppm`)

#relacion Hierro Cobre

Relaciones_binarias$`Fe ppm/Cu`<-(Relaciones_binarias$`Fe ppm`)/Relaciones_binarias$`Cu ppm`
Relaciones_binarias$`Cu/Fe ppm`<-Relaciones_binarias$`Cu ppm`/(Relaciones_binarias$`Fe ppm`)

relaciones con Boro. como Boro es presentado en el analisi de suelos como ppm,no es necesario multiplicar por 100 para las variables que se presentan como ppm. como ya se han calculado las relaciones con Nitrogeno, Fosforo,Potasio,Calcio,Magnesio,Azufre y Hierro estas se han omitido.

#relacion Boro Zinc

Relaciones_binarias$`B ppm/Zn`<-(Relaciones_binarias$`B ppm`)/Relaciones_binarias$`Zn ppm`
Relaciones_binarias$`Zn/B ppm`<-Relaciones_binarias$`Zn ppm`/(Relaciones_binarias$`B ppm`)

#relacion Boro Manganeso

Relaciones_binarias$`B ppm/Mn`<-(Relaciones_binarias$`B ppm`)/Relaciones_binarias$`Mn ppm`
Relaciones_binarias$`Mn/B ppm`<-Relaciones_binarias$`Mn ppm`/(Relaciones_binarias$`B ppm`)

#relacion Boro Cobre

Relaciones_binarias$`B ppm/Cu`<-(Relaciones_binarias$`B ppm`)/Relaciones_binarias$`Cu ppm`
Relaciones_binarias$`Cu/B ppm`<-Relaciones_binarias$`Cu ppm`/(Relaciones_binarias$`B ppm`)

relaciones con Zinc. como Zinc es presentado en el analisi de suelos como ppm,no es necesario multiplicar por 100 para las variables que se presentan como ppm. como ya se han calculado las relaciones con Nitrogeno, Fosforo,Potasio,Calcio,Magnesio,Azufre,Hierro y Boro estas se han omitido.

#relacion Zinc  Manganeso

Relaciones_binarias$`Zn ppm/Mn`<-(Relaciones_binarias$`Zn ppm`)/Relaciones_binarias$`Mn ppm`
Relaciones_binarias$`Mn/Zn ppm`<-Relaciones_binarias$`Mn ppm`/(Relaciones_binarias$`Zn ppm`)

#relacion Zinc Cobre

Relaciones_binarias$`Zn ppm/Cu`<-(Relaciones_binarias$`Zn ppm`)/Relaciones_binarias$`Cu ppm`
Relaciones_binarias$`Cu/Zn ppm`<-Relaciones_binarias$`Cu ppm`/(Relaciones_binarias$`Zn ppm`)

relaciones con Cobre y Manganeso. como Cobre y Manganeso es presentado en el analisi de suelos como ppm,no es necesario multiplicar por 100 para las variables que se presentan como ppm. como ya se han calculado las relaciones con Nitrogeno, Fosforo,Potasio,Calcio,Magnesio,Azufre,Hierro,Boro y Zinc

#relacion Manganeso Cobre

Relaciones_binarias$`Mn ppm/Cu`<-(Relaciones_binarias$`Mn ppm`)/Relaciones_binarias$`Cu ppm`
Relaciones_binarias$`Cu/Mn ppm`<-Relaciones_binarias$`Cu ppm`/(Relaciones_binarias$`Mn ppm`)
head(Relaciones_binarias)

si desearamos corregir datos del data.frame

##  fix(Relaciones_binarias)

estadisticas de todas las columnas

seleccion<-Relaciones_binarias[,14:123]
media<-apply(seleccion,2,mean)
varianza<-apply(seleccion,2,var)
desviacion_estandar<-apply(seleccion,2,sd)

data.frame con la estasdisticas

estadisticas<-data.frame(media,varianza,desviacion_estandar)

agregamos la columna coeficiente de variacion que nos servira mas adelante para el analisis DRIS

estadisticas$coeficiente_de_variacion<-estadisticas$desviacion_estandar/estadisticas$media

para arreglar la tabla de estadisticas en una manera comoda de ver utilizamos la funcion t(),para transponer los datos

estadisticas_en_filas<-t(estadisticas)