print("not paired")

a=read.csv("http://www.ic.unicamp.br/~wainer/cursos/1s2018/metodologia/a1.csv",header = F)
b=read.csv("http://www.ic.unicamp.br/~wainer/cursos/1s2018/metodologia/b1.csv",header=F)


print(t.test(a$V1,b$V1)$p.value)
print(wilcox.test(a$V1,b$V1)$p.value)
hist(a$V1,10)
hist(b$V1,10)

one=read.csv("http://www.ic.unicamp.br/~wainer/cursos/1s2018/metodologia/paired.csv",header=F)

print("paired")
print(t.test(one$V1,one$V2,paired=T)$p.value)
print(wilcox.test(one$V1,one$V2,paired=T)$p.value)
print("non paired on the paired data")
print(t.test(one$V1,one$V2,paired=F)$p.value)
print(wilcox.test(one$V1,one$V2,paired=F)$p.value)


nV1=sum(one$V1>one$V2)
nn=length(one$V1)
print(binom.test(nV1,nn)$p.value)

print("estudo do pvalue")

print("first case")
pval=NULL
for (i in 1:10){
  x1=rnorm(15,10,5)
  x2=rnorm(15,13,5)

  tt=t.test(x1,x2)$p.value
  pval=c(pval,tt)
}
print(mean(pval))

print("2nd case")
pval=NULL
for (i in 1:10){
  x1=rnorm(25,10,5)
  x2=rnorm(25,13,5)
  
  tt=t.test(x1,x2)$p.value
  pval=c(pval,tt)
}
print(mean(pval))


print("3rd case")
pval=NULL
for (i in 1:10){
  x1=rnorm(15,10,5)
  x2=rnorm(15,17,5)
  
  tt=t.test(x1,x2)$p.value
  pval=c(pval,tt)
}
print(mean(pval))

print("4th case")
pval=NULL
for (i in 1:10){
  x1=rnorm(15,10,8)
  x2=rnorm(15,13,8)
  
  tt=t.test(x1,x2)$p.value
  pval=c(pval,tt)
}
print(mean(pval))