##############################
# concentracion
###################################
conc<-scan(text='
617.2 617.2 617.3 617.4 617.4 617.5 617.6 617.6 617.6 617.7 617.7 617.7 617.7 617.8 617.8 617.8
618.0 618.0 618.2 618.5 619.9 621.9 623.7 626.7 628.1 632.6 648.0 652.7')
 
#pdf("conc.pdf")
par(mfrow=c(2,2))  #divide la pantalla grafica en 4
boxplot(conc,main="Boxplot conc")
hist(conc,probability=T)
qqnorm(conc)
par(mfrow=c(1,1)) #vuelve la pantalla grafica a un solo grafico
#dev.off()
shapiro.test(conc)
(conc-620>0)
sum(conc-620>0)
#pvalor
pbinom(7,length(conc),prob=0.5)
library(BSDA)
SIGN.test(conc, md = 620, alternative = "less", conf.level = 0.95)

t.test(conc,mu = 620)   #podemos usar este test?

###################################
# suero
###################################
suero<-scan(text = '
35 36 37 37 39 44 48 56 60 76 81 83 83 84 85')

boxplot(suero)
hist(suero,probability=T)
qqnorm(suero)
shapiro.test(suero)
(suero-40)
abs(suero-40)
rangos<-rank(abs(suero-40))
(suero-40)>0
positivas<-((suero-40)>0)
positivas
estadistico<-sum(rangos[positivas])  #sumamos los rangos de las observaciones mayores a 40
estadistico
wilcox.test(suero ,alternative ="two.sided", mu=40,exact=FALSE)
SIGN.test(suero ,alternative ="two.sided", md=40)  #que conclusion obtiene?
t.test(suero ,alternative ="two.sided", mu=40)  #es correcto usarlo?


###################################
# perfumes
###################################

Moder<-scan(text = '
 20 25 40 44 43 13 32 34 35 11 12 46 13 17 47')

Trad<-scan(text = '
 5  7 10 11 12 17 21 28 33 35 40 40 41 44 45')

shapiro.test(Moder)
shapiro.test(Trad)

par(mfrow=c(2,2))
hist(Moder,probability=T)
hist(Trad,probability=T)
boxplot(Moder,Trad,names=c("moderno","tradicional"))
par(mfrow=c(1,1))
wilcox.test(Moder, Trad, paires=FALSE,alternative = "two.sided")
wilcox.test(Moder, Trad, paires=FALSE,alternative = "two.sided",exact=FALSE)
t.test(Moder,Trad, paired=FALSE,alternative = "two.sided") 
#se puede usar el test t? Coincide la conclusion con el anterior?