#Se desea conocer de que manera repercute en una plantacion de frutales, riego adicional en los cultivos (20L,35L,50L) respecto a la cantidad habitual,donde se aplicaron 3 tratamientos en una parcel 4x10, donde su variable respuesta es Rendimiento(rto) y su covariable es porcentaje de humedad del suelo(hmdad)

####Datos trabajados de excel

library(readxl)
library(knitr)
library(lattice)
library(ape)
library(ggplot2)
library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
datosl <- read_excel("D:/Materias 2021/Dexperimentos/Tareas/tablaviernes.xlsx",
                     sheet = "LUZ")
datosl$tratamientos =factor(datosl$tratamientos)

###Grafico De crecimientos (Resultado)

g1=datosl %>% ggplot(aes(x = x, y=y, colour=rto))+
   geom_point(size = 10,shape=19)+
   scale_fill_continuous(type = "gradient")+
  labs(title = "RENDIMIENTO DEL CULTIVO")+
  xlab("")+ylab("")+theme (text = element_text(size=18));g1

###Grafico de cantidad de frutos(COvariable)

g2=datosl %>% ggplot(aes(x = x, y=y, colour=hmdad))+
   geom_point(size = 10,shape=19)+scale_color_continuous(type = 'viridis') + labs(title = "PORCENTAJE DE HUMEDAD EN EL SUELO (%)")+
  xlab("")+ylab("")+theme (text = element_text(size=15));g2

####Medias de tratamiento con el Rendimiento

medr=tapply(datosl$rto,datosl$tratamientos,mean);medr
##   T0   T1   T2   T3 
## 0.31 0.71 1.19 1.74

####Medias de tratamiento con los porcentajes de humedad

medh=tapply(datosl$hmdad,datosl$tratamientos,mean);medh
##   T0   T1   T2   T3 
## 55.7 64.2 71.9 84.0
library(lattice)
g3=xyplot(datosl$rto~datosl$hmdad,pch=18);g3

g4=xyplot(datosl$rto~datosl$hmdad|datosl$tratamientos,pch=19);g4

###Hipotesis de lo tratamientos,¿Son iguales los tratamientos?… \[[H_0: \mu_{T_0} = \mu_{T_1} =\mu_{T_2} =\mu_{T_3}]\] ###Hipotesis de la covariable,¿La covariable influye en los tratamientos?… \[[H_0:\delta =0 ]\]

ANCOVA_A <- aov(datosl$rto~datosl$hmdad+datosl$tratamientos)
summary(ANCOVA_A)
##                     Df Sum Sq Mean Sq F value   Pr(>F)    
## datosl$hmdad         1 10.659  10.659  561.87  < 2e-16 ***
## datosl$tratamientos  3  0.781   0.260   13.72 4.48e-06 ***
## Residuals           35  0.664   0.019                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Debido a que el valor p-valor de los tratamientos y la covaiable es menor al 5% se rechaza la Hipotesis nula por lo que los tratamientos son diferentes y la covariable si influye en los tratamientos

##La Prueba Shapiro-Wilk normality test mostro normalidad en los residuales

shapiro.test(ANCOVA_A$residuals)
## 
##  Shapiro-Wilk normality test
## 
## data:  ANCOVA_A$residuals
## W = 0.97162, p-value = 0.4046

#La prueba Bartlett test of homogeneity of variances mostro que las varianzas son similares

bartlett.test(datosl$rto,datosl$tratamientos)
## 
##  Bartlett test of homogeneity of variances
## 
## data:  datosl$rto and datosl$tratamientos
## Bartlett's K-squared = 0.28081, df = 3, p-value = 0.9636
ANCOVA.dists <- as.matrix(dist(cbind(datosl$x, datosl$y)))

ANCOVA.dists.inv <- 1/ANCOVA.dists
diag(ANCOVA.dists.inv) <- 0

###Sin Covariable (Moran)

inds=Moran.I(datosl$rto, ANCOVA.dists.inv);inds
## $observed
## [1] -0.06959209
## 
## $expected
## [1] -0.02564103
## 
## $sd
## [1] 0.02397495
## 
## $p.value
## [1] 0.06677161

###Con Covariable (Moran)

indc=Moran.I(datosl$hmdad, ANCOVA.dists.inv);indc
## $observed
## [1] -0.06500269
## 
## $expected
## [1] -0.02564103
## 
## $sd
## [1] 0.0239125
## 
## $p.value
## [1] 0.09974921

#Relacion calcio magnesio en el suelo

##1.El calcio participa en la elongación y división celular, permeabilidad de las membranas y activación de algunas enzimas críticas para el desarrollo.

##2.Es fundamental en al estructura del suelo debido a que es el cation (Ca2+)que predomina en el suelo bajo condiciones de suelos con poco Aluminio(Al3+)es decir suelos basicos(El calcio se usa para tratar suelos acidos mediante la tecnica del “encalado” ).

##3.Debido al comportamiento del calcio y ser poco movil en la la planta, en el momento en que hay una deficiencia se pueden observar en las zona meristematicas de las plantas donde si baja mucho el nivel de calcio los metales como Mg,Mn,Zn suben a niveles toxicos

##4.La cantidad Ca y Mg en el suelo esta directamente relacionada con la meteorización de los minerales y el grado de lixiviación

##5.El magnesio hace parte de la molecula de la clorofila, debido a esto influye en la fotosintesis de las plantas, ademas de suma importancia en la sintesis de aceites, proteinas y la actividad enzimatica.

##6.En el momento en que hay un deficiencia de Mg(magnesio), el sintoma mas caracteristico en la clorosis internerval en las hojas viejas de las plantas, esta deficiencia se da priincipalmente en suelos arenosos o suelos de baja CIC.

##7.La mejor relación Calcio/Magnesio (Ca/Mg) es de 1/1 a 15/1.

##8.En un sistema forestal el principal metodo de perdida de Ca se da mediante la deforestacion o tala de arboles

##9.Hay mas probabilidad de tener una deficiencia de Magnesio que de calcio , debido a que se encuentra en el suelo en menor medida

##10.Las enmiendas son utiles para balancear las pérdidas de Calcio y Magnesio,donde no sólo mantienen los niveles de Ca y Mg sino que también generan una estabilidad física y química en el suelo.

##https://www.edafologia.org/app/download/7956234176/Calcio+y+Magnesio+Xi.pdf?t=1563476239