##########ejercicio 1 ###################################
conc<-scan("concentracion.txt")
boxplot(conc)
hist(conc)
shapiro.test(conc)
binom.test(sum(conc-620>0),length(conc) , p = 0.5,alternative ="less",conf.level = 0.95)
library(BSDA)
sign.test(conc, md = 620, alternative = "less", conf.level = 0.95)

########### ejercicio 2 ###################################
precios<-scan()
21.92  25.75  28.41  33.00  45.00  49.79  51.47  53.55  62.91  63.90  75.01  96.99 102.52 124.33 129.44 135.76 143.42 144.93 156.99 159.99

hist(precios,probability = T)
qqnorm(precios)
shapiro.test(precios)
wilcox.test(precios,mu=120,alternative ="two.sided")
#calculemos a mano el estadístico
precios-120
precios-120>0
abs(precios-120)
rank(abs(precios-120)) #rangos del valor absoluto de los D_i
sum(rank(abs(precios-120))[ precios-120>0]) #es el valor del estadístico del test de mann whitney wilcoxon para una muestra
sign.test(precios,md = 120)

###########    ejercicio 3   ####################
paracetamol<-read.table(file.choose(),header=TRUE) #paracetamol.txt
attach(paracetamol)
difpara<-UV-Infrarrojo
hist(difpara, probability = TRUE)
t.test(difpara,alternative = c("two.sided"),mu = 0)
t.test(UV,Infrarrojo,paired = TRUE, alternative = c("two.sided"),mu = 0)


############ ejercicio 4  ########################
energia<-read.table("obesas y delgadas.txt",header=TRUE)
attach(energia)
summary(energia)
boxplot(energia)
hist(Del)
hist(Ob)
stem(Del)
stem(Ob)
shapiro.test(Del)
shapiro.test(Ob)
var.test(Del,Ob, ratio = 1, alternative ="two.sided",conf.level = 0.95)
t.test(Del,Ob , alternative = c("two.sided"),mu = 0,paired = FALSE, var.equal=TRUE,conf.level = 0.95)

t.test(Del,Ob , alternative = c("two.sided"),mu = 0,paired = FALSE, var.equal=FALSE,conf.level = 0.95)


############## ejercicio 5 ###############################
perfumes<-read.table("perfumes.txt",header=TRUE)
attach(perfumes)
shapiro.test(Moder)
shapiro.test(Trad)
boxplot(Moder)
boxplot(Trad)
hist(Moder)
hist(Trad)
wilcox.test(Moder, Trad, paired=FALSE,alternative = "two.sided")
wilcox.test(Moder, Trad, paired=FALSE,alternative = "two.sided",exact=FALSE)

###########ejercicio 6#########
zapatos<-read.table("zapatos.txt",header=TRUE)
attach(zapatos)
boxplot(zapatos)
qqnorm(A)
qqnorm(B)
qqnorm(A-B)
boxplot(A-B)
shapiro.test(A-B)

###########ejercicio 7#########
claro<-scan()
22   29   13   16   15   18   14   6
 
denso<-scan()
17   20   170   315   22   190   64

par(mfrow=c(1,2))  #dividimos la pantalla gráfica en dos
hist(claro)
hist(denso)

qqnorm(claro)
qqnorm(denso)
par(mfrow=c(1,1)) #restauramos la pantalla gráfica a su forma habitual
wilcox.test(denso, claro)

###########ejercicio 8#########
aa<-scan()
 8 7 4 2 4 10 17 3 2 11

bb<-scan()
 10 6 4 5 7 11 15 6 3 14

hist(aa-bb)
wilcox.test(aa,bb,paired=T,alternative="two.sided")
wilcox.test(aa-bb,alternative="two.sided")
library(BSDA)
sign.test(aa-bb)
t.test(aa,bb,paired=TRUE)  #puede usarse esta instruccion?