Censo 70a-71b

library(readxl)
p70a <- read_excel("deficiencia clorofiliana y desprendimiento basal(Parcelas 70a 71b).xlsx", 
    sheet = "2070a")
View(p70a)

library(readxl)
p71b <- read_excel("deficiencia clorofiliana y desprendimiento basal(Parcelas 70a 71b).xlsx", 
    sheet = "2071b")
View(p71b)

#incidencia
desprendimiento_71b=table(p71b$`Des`)
desprendimiento_71b

## 
##  No  Si 
##  85 245

incidencia_71b_desprendimeinto=round((prop.table(desprendimiento_71b)*100),2)
incidencia_71b_desprendimeinto

## 
##    No    Si 
## 25.76 74.24

library(Hmisc)

## 
## Attaching package: 'Hmisc'

## The following objects are masked from 'package:base':
## 
##     format.pval, units

desprendimeinto_71b=describe(p71b$`Des`)
desprendimeinto_71b

## p71b$Des 
##        n  missing distinct 
##      330        0        2 
##                       
## Value         No    Si
## Frequency     85   245
## Proportion 0.258 0.742

daño=describe(p71b$`Daño`)
daño

## p71b$Daño 
##        n  missing distinct 
##      330        0        2 
##                       
## Value         No    Si
## Frequency    164   166
## Proportion 0.497 0.503

Amarillamiento=describe(p71b$`Amarillamiento`)
Amarillamiento

## p71b$Amarillamiento 
##        n  missing distinct 
##      330        0        2 
##                       
## Value         No    Si
## Frequency    126   204
## Proportion 0.382 0.618

tabcon <- array(data = c(180,10, 108,90,50, 2, 12, 20), 
                   dim = c(2,2,2), 
                   dimnames = list("Cig" = c("S","N"),
                                   "Mar" = c("S","N"),
                                   "Alc" = c("S","N")))
tabcon

## , , Alc = S
## 
##    Mar
## Cig   S   N
##   S 180 108
##   N  10  90
## 
## , , Alc = N
## 
##    Mar
## Cig  S  N
##   S 50 12
##   N  2 20

ftable(tabcon, row.vars = c("Alc","Cig","Mar"))

## Alc Cig Mar     
## S   S   S    180
##         N    108
##     N   S     10
##         N     90
## N   S   S     50
##         N     12
##     N   S      2
##         N     20