ODDs <- function(coeficientes, individuo, logit=1){
	z = sum(coeficientes*individuo)
	if (logit == 1){
		p = exp(z)/(1+exp(z))
	}
	else { 
		p = pnorm(z)
	}
	odd = p/(1-p)
	salida = data.frame(z, p, odd)
	return(salida)
}

#######################################################################

ODD_RATIO <- function(coeficientes, individuo1, individuo2, logit=1){
	odd1 = ODDs(coeficientes, individuo1, logit)[[3]]
	odd2 = ODDs(coeficientes, individuo2, logit)[[3]]
	odd_ratio = odd1/odd2
	return(odd_ratio)
}

#######################################################################

TASA_ACIERTOS <- function(umbral, valores_estimados, valores_observados)
{
  A = 0
  B = 0
  C = 0
  D = 0
  n = length(valores_observados)
  for (i in 1:n){
    if ((valores_estimados[i]>=umbral) & valores_observados[i]==1){A = A + 1}
    if ((valores_estimados[i]>=umbral) & valores_observados[i]==0){B = B + 1}
    if ((valores_estimados[i]<umbral) & valores_observados[i]==1){C = C + 1}
    if ((valores_estimados[i]<umbral) & valores_observados[i]==0){D = D + 1}
  }
  exito_aciertos_A = (A/n)*100
  exito_errores_B = (B/n)*100
  fracaso_errores_C = (C/n)*100
  fracaso_aciertos_D = (D/n)*100
  tasa_aciertos = ((A+D)/n)*100
  salida = data.frame(exito_aciertos_A, exito_errores_B, fracaso_errores_C, fracaso_aciertos_D, tasa_aciertos)
  return(salida)
}