Objetivo general

Estimar el área pélvica a través de medidas externas y su relación con facilidad de parto y peso al nacer en un rebaño del trópico Alto.

Figura 1. abertura de la pelvis craneal se estrecha ligeramente ventralmente, y tiene la forma de un trapecio con ángulos redondeados.

Objetivos especificos

• Estimar el área pélvica usando las metodologías de Murray 2002
• Establecer la relación entre área pélvica de vientres (novillas de vientre y vacas), facilidad de parto y peso al nacer de las crías.
• Evaluar el Puntaje de condición corporal preparto de vacas y novillas en un hato de trópico Alto.
• Determinar el peso al nacer de crías y su escala de facilidad de parto en novillas y vacas Holstein en un hato de trópico Alto.

• Animales El estudio se efectuo con vacas lecheras de la raza Holstein con sus cruces (1/2Hol1/2Nor y 3/4Nor1/4Hol), se hizo la selección del rebaño, siendo escogidas por medio del software ganadero un grupo de vientres con fecha probable de parto en marzo a agosto de 2019.

• La investigación se realizo en el CIC Santa María del Puyón de la Universidad de la Salle ubicado en Sopo (Cundinamarca)
• Altitud media de 2587 m s. n. m, clima promedio de 14°C
• Dos periodos de lluvia al año: marzo a junio y septiembre a noviembre, promedio de lluvia anual entre 500 a 1000mm.

• Según metodóloga sugerida por Murray 2002 Área pélvica = - 122.2 + 23.2 x (ancho de entre coxales) + 24.3 (largo de coxales) - 0.3 x Ancho entre coxales x Longitud de coxales

Evaluar la diferencia entre área pélvica segun Murray (Basada en una ecucion realizada anteriormente para este valor) y la variabilidad que se encuentra en la razas holstein y en los cruces, para esto se hara uso de las metodologia de Anova, Kruskall Wallis.

Ho: las medidas pelvianas de la poblacion son iguales, por lo tanto, La raza no influyen en el área pélvica de las vacas.

Ha: Hay diferencia en al menos una raza (cruce), por consiguiente la raza influye en el tamaño del área pélvica de las vacas lecheras - Variable Cuantitativa: Area pelvica Murray 2002 - Variable cualitativa: Raza

library(readxl)
pelvimetria<-read_excel("datos pelvimetria.xlsx")
attach(pelvimetria)
names(pelvimetria)

##  [1] "ID"                       "Categoria"               
##  [3] "Raza"                     "Peso"                    
##  [5] "Altura"                   "Ancho anca"              
##  [7] "Largo anca"               "Area entre coxales"      
##  [9] "Condicion corporal"       "Area pelvica Murray 2002"

class(Raza)

## [1] "character"

class(`Area pelvica Murray 2002`)

## [1] "numeric"

boxplot(`Area pelvica Murray 2002`~Raza)

aov(`Area pelvica Murray 2002`~Raza)

## Call:
##    aov(formula = `Area pelvica Murray 2002` ~ Raza)
## 
## Terms:
##                     Raza Residuals
## Sum of Squares   7517.60  96567.59
## Deg. of Freedom        2        20
## 
## Residual standard error: 69.48654
## Estimated effects may be unbalanced
## 1 observation deleted due to missingness

ANOVA=aov(`Area pelvica Murray 2002`~Raza)
summary(ANOVA)

##             Df Sum Sq Mean Sq F value Pr(>F)
## Raza         2   7518    3759   0.778  0.473
## Residuals   20  96568    4828               
## 1 observation deleted due to missingness

• De esta manera se acepta la hipotesis nula ya que el P-valor es >0.5 por lo tanto no hay ninguna relación entre la raza y la variavilidad de area pelvia de las vacas. no hay necesidad de realizar la prueba de comparaciones multiples.

library(lmtest)

## Loading required package: zoo

## 
## Attaching package: 'zoo'

## The following objects are masked from 'package:base':
## 
##     as.Date, as.Date.numeric

Prueba de normalidad

errores=residuals(ANOVA)
shapiro.test(errores)

## 
##  Shapiro-Wilk normality test
## 
## data:  errores
## W = 0.95769, p-value = 0.4183

• Pueba de homocedasticidad

bptest(`Area pelvica Murray 2002`~Raza)

## 
##  studentized Breusch-Pagan test
## 
## data:  `Area pelvica Murray 2002` ~ Raza
## BP = 0.52832, df = 2, p-value = 0.7678

• prueba de independencia

dwtest(ANOVA,alternative = "two.sided")

## 
##  Durbin-Watson test
## 
## data:  ANOVA
## DW = 2.1509, p-value = 0.6605
## alternative hypothesis: true autocorrelation is not 0

Se cumplen los supuestos

• Con un nivel de significancia de 0.05 no se rechaza la Ho nula, ya que no hay diferencia en ningún método usado.

Esta prueba sera usada como alternativa al Anova

• Ho: No existe diferencia significativa entre la raza con el area pelvica de las vacas

• Ha: Existen diferencias significativas entre la raza con el area pelvica de las vacas

kruskal.test(`Area pelvica Murray 2002`~Raza)

## 
##  Kruskal-Wallis rank sum test
## 
## data:  Area pelvica Murray 2002 by Raza
## Kruskal-Wallis chi-squared = 2.1442, df = 2, p-value = 0.3423

Por lo tanto, se indica que el p-value = 0.3423 es mayor a 0.05, no se rechaza la Ho nula, lo que indica que no hay diferencia significativa entre la raza con el area pelvica.

Para esta prueba se tomo como variables: la categoria (novilla, vaca) y el área pelvica, mediendo la relacion de la categoria y su influencia en el tamaño del área pelviana.

Test U de Mann-Whitney

• Ho: Mediana (1) = Mediana (2), El área pelvica es igual en las novillas y vacas
• Ha:Mediana (1) ≠ Mediana (2), El área pelvica es diferente en las novilla y vacas

library(BSDA)

## Loading required package: lattice

## 
## Attaching package: 'BSDA'

## The following object is masked from 'package:datasets':
## 
##     Orange

Novilla=c(1247,1315,1180,1290,1281,1238,1298,1306,1315,1338,1107,1293,1240)
vaca=c(1347,1402,1362,1372,1259,1415,1275,1351,1307,1315)
wilcox.test(Novilla,vaca,paired = F)

## Warning in wilcox.test.default(Novilla, vaca, paired = F): cannot compute
## exact p-value with ties

## 
##  Wilcoxon rank sum test with continuity correction
## 
## data:  Novilla and vaca
## W = 21, p-value = 0.006925
## alternative hypothesis: true location shift is not equal to 0

• Con un nivel de significacia de 0.05, se rechaza la Ho nula y se acepta la alternativa, ya que si hubo diferecia en la medida área pelvica de novillas y vacas.