## Script: estadMsist.R
## Sobre?: estadigrafos muestreo sistematico
##
## Profesor: Christian Salas Eljatib
## E-mail: cseljatib AT gmail DOT com
## Web: https://eljatib.com
#======================================

##ejercicio de apuntes de clases, sobre
# numero de plantulas infectadas por hongo/metro lineal
# en un vivero forestal.
y <- c(98, 85, 71, 97, 87, 94,
       90, 71, 89, 72, 81, 81, 80,
       84, 83, 88, 90, 51, 96, 78,
       100, 79, 51, 88, 54, 66, 86,
       62)
n <- length(y)
N <- 140

a<- 5

m.y<-mean(y)
m.y

m.y <- (a/N)*sum(y)
m.y


##
var.y <- var(y)
var.y

sd.y <- sd(y)
sd.y

cv.y <- sd.y*100/m.y
cv.y

f<- n/N
##Error estandar estimador-MAS
se.med.y <- sqrt((var.y/n) *
                   (1 - f))
se.med.y

se.med.y.p<-se.med.y*100/m.y
se.med.y.p


# Error estandar estimador-MSist
var.sd <- sum((diff(y))^2)/(2 *
                              (n - 1)) * ((1/n) - (1/N))
se.sd <- sqrt(var.sd)
se.sd

##error de muestreo
df <- n - 1
conf <- 90
alpha <- 1 - (conf/100)
alpha.2 <- alpha/2
t.value <- abs(qt(1 - alpha.2,
                  df))
t.value

EM<- se.sd*t.value
EM #en plantas/mlineal

EM2<- se.med.y*t.value
EM2 #en plantas/mlineal

#Intervalo confidencial
#intervalo de confianza
lim.inf<-m.y - EM
lim.sup<-m.y + EM
lim.inf2<-m.y - EM2
lim.sup2<-m.y + EM2

SE.MAS<-c(lim.inf2, lim.sup2)
SE.Msist<-c(lim.inf, lim.sup)

rbind(SE.MAS,SE.Msist)

#@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
#*´¨)
#¸.•´¸.•*´¨) ¸.•*¨)
#(¸.•´ (¸.•` ¤ Fin del script