## Objetivo: criar função que crie uma tabela com os resultados (valores de p)

# A tabela deve ter 15 tempos (15 tempos diferentes possíveis),
# 7 linhas (7 fatores e combinações de fatores)
# e 10 colunas (2 espécies e 5 medidas)

# Fatores (1 a 7) - luz, ser, com, luz/ser, luz/com, ser/comp, luz/ser/com
# Medidas (6 a 15) - g.dsh, g.h, g.nl, g.hi, g.tm, r.dsh, ..., r.tm
# Tempos (1 a 15) - 01,02,04,06,07,12,14,16,17,24,26,27,46,47,67

res=function(d){

result=data.frame()
result[7,1]=0
result[7,2]=0
result[7,3]=0
result[7,4]=0
result[7,5]=0
result[7,6]=0
result[7,7]=0
result[7,8]=0
result[7,9]=0
result[7,10]=0
result[7,11]=0
result[,1]=c("luz","ser","com","luz/ser","luz/com","ser/com","luz/ser/com")
names(result)=c("fator","g.dsh","g.h","g.nl","g.hi","g.tm","r.dsh","r.h","r.nl","r.hi","r.tm")

## Tempo 

a6=aleat(999,6,d)
a7=aleat(999,7,d)
a8=aleat(999,8,d)
a9=aleat(999,9,d)
a10=aleat(999,10,d)
a11=aleat(999,11,d)
a12=aleat(999,12,d)
a13=aleat(999,13,d)
a14=aleat(999,14,d)
a15=aleat(999,15,d)

# Medida 6 (g.dsh); fixo: a6

result[1,2]=length(sort(a6[a6[,1]>=a6[1000,1],1]))/1000 # luz
result[2,2]=length(sort(a6[a6[,2]>=a6[1000,2],2]))/1000 # ser
result[3,2]=length(sort(a6[a6[,3]>=a6[1000,3],3]))/1000 # com
result[4,2]=length(sort(a6[a6[,4]>=a6[1000,4],4]))/1000 # luz/ser
result[5,2]=length(sort(a6[a6[,5]>=a6[1000,5],5]))/1000 # luz/com
result[6,2]=length(sort(a6[a6[,6]>=a6[1000,6],6]))/1000 # ser/com
result[7,2]=length(sort(a6[a6[,7]>=a6[1000,7],7]))/1000 # luz/ser/com

# Medida 7 (g.h); fixo: a7

result[1,3]=length(sort(a7[a7[,1]>=a7[1000,1],1]))/1000 # luz
result[2,3]=length(sort(a7[a7[,2]>=a7[1000,2],2]))/1000 # ser
result[3,3]=length(sort(a7[a7[,3]>=a7[1000,3],3]))/1000 # com
result[4,3]=length(sort(a7[a7[,4]>=a7[1000,4],4]))/1000 # luz/ser
result[5,3]=length(sort(a7[a7[,5]>=a7[1000,5],5]))/1000 # luz/com
result[6,3]=length(sort(a7[a7[,6]>=a7[1000,6],6]))/1000 # ser/com
result[7,3]=length(sort(a7[a7[,7]>=a7[1000,7],7]))/1000 # luz/ser/com

# Medida 8 (g.nl); fixo: a8

result[1,4]=length(sort(a8[a8[,1]>=a8[1000,1],1]))/1000 # luz
result[2,4]=length(sort(a8[a8[,2]>=a8[1000,2],2]))/1000 # ser
result[3,4]=length(sort(a8[a8[,3]>=a8[1000,3],3]))/1000 # com
result[4,4]=length(sort(a8[a8[,4]>=a8[1000,4],4]))/1000 # luz/ser
result[5,4]=length(sort(a8[a8[,5]>=a8[1000,5],5]))/1000 # luz/com
result[6,4]=length(sort(a8[a8[,6]>=a8[1000,6],6]))/1000 # ser/com
result[7,4]=length(sort(a8[a8[,7]>=a8[1000,7],7]))/1000 # luz/ser/com

# Medida 9 (g.hi); fixo: a9

result[1,5]=length(sort(a9[a9[,1]>=a9[1000,1],1]))/1000 # luz
result[2,5]=length(sort(a9[a9[,2]>=a9[1000,2],2]))/1000 # ser
result[3,5]=length(sort(a9[a9[,3]>=a9[1000,3],3]))/1000 # com
result[4,5]=length(sort(a9[a9[,4]>=a9[1000,4],4]))/1000 # luz/ser
result[5,5]=length(sort(a9[a9[,5]>=a9[1000,5],5]))/1000 # luz/com
result[6,5]=length(sort(a9[a9[,6]>=a9[1000,6],6]))/1000 # ser/com
result[7,5]=length(sort(a9[a9[,7]>=a9[1000,7],7]))/1000 # luz/ser/com

# Medida 10 (g.tm); fixo: a10

result[1,6]=length(sort(a10[a10[,1]>=a10[1000,1],1]))/1000 # luz
result[2,6]=length(sort(a10[a10[,2]>=a10[1000,2],2]))/1000 # ser
result[3,6]=length(sort(a10[a10[,3]>=a10[1000,3],3]))/1000 # com
result[4,6]=length(sort(a10[a10[,4]>=a10[1000,4],4]))/1000 # luz/ser
result[5,6]=length(sort(a10[a10[,5]>=a10[1000,5],5]))/1000 # luz/com
result[6,6]=length(sort(a10[a10[,6]>=a10[1000,6],6]))/1000 # ser/com
result[7,6]=length(sort(a10[a10[,7]>=a10[1000,7],7]))/1000 # luz/ser/com

# Medida 11 (r.dsh); fixo: a11

result[1,7]=length(sort(a11[a11[,1]>=a11[1000,1],1]))/1000 # luz
result[2,7]=length(sort(a11[a11[,2]>=a11[1000,2],2]))/1000 # ser
result[3,7]=length(sort(a11[a11[,3]>=a11[1000,3],3]))/1000 # com
result[4,7]=length(sort(a11[a11[,4]>=a11[1000,4],4]))/1000 # luz/ser
result[5,7]=length(sort(a11[a11[,5]>=a11[1000,5],5]))/1000 # luz/com
result[6,7]=length(sort(a11[a11[,6]>=a11[1000,6],6]))/1000 # ser/com
result[7,7]=length(sort(a11[a11[,7]>=a11[1000,7],7]))/1000 # luz/ser/com

# Medida 12 (r.h); fixo: a12

result[1,8]=length(sort(a12[a12[,1]>=a12[1000,1],1]))/1000 # luz
result[2,8]=length(sort(a12[a12[,2]>=a12[1000,2],2]))/1000 # ser
result[3,8]=length(sort(a12[a12[,3]>=a12[1000,3],3]))/1000 # com
result[4,8]=length(sort(a12[a12[,4]>=a12[1000,4],4]))/1000 # luz/ser
result[5,8]=length(sort(a12[a12[,5]>=a12[1000,5],5]))/1000 # luz/com
result[6,8]=length(sort(a12[a12[,6]>=a12[1000,6],6]))/1000 # ser/com
result[7,8]=length(sort(a12[a12[,7]>=a12[1000,7],7]))/1000 # luz/ser/com

# Medida 13 (r.nl); fixo: a13

result[1,9]=length(sort(a13[a13[,1]>=a13[1000,1],1]))/1000 # luz
result[2,9]=length(sort(a13[a13[,2]>=a13[1000,2],2]))/1000 # ser
result[3,9]=length(sort(a13[a13[,3]>=a13[1000,3],3]))/1000 # com
result[4,9]=length(sort(a13[a13[,4]>=a13[1000,4],4]))/1000 # luz/ser
result[5,9]=length(sort(a13[a13[,5]>=a13[1000,5],5]))/1000 # luz/com
result[6,9]=length(sort(a13[a13[,6]>=a13[1000,6],6]))/1000 # ser/com
result[7,9]=length(sort(a13[a13[,7]>=a13[1000,7],7]))/1000 # luz/ser/com

# Medida 14 (r.hi); fixo: a14

result[1,10]=length(sort(a14[a14[,1]>=a14[1000,1],1]))/1000 # luz
result[2,10]=length(sort(a14[a14[,2]>=a14[1000,2],2]))/1000 # ser
result[3,10]=length(sort(a14[a14[,3]>=a14[1000,3],3]))/1000 # com
result[4,10]=length(sort(a14[a14[,4]>=a14[1000,4],4]))/1000 # luz/ser
result[5,10]=length(sort(a14[a14[,5]>=a14[1000,5],5]))/1000 # luz/com
result[6,10]=length(sort(a14[a14[,6]>=a14[1000,6],6]))/1000 # ser/com
result[7,10]=length(sort(a14[a14[,7]>=a14[1000,7],7]))/1000 # luz/ser/com

# Medida 15 (r.tm); fixo: a15

result[1,11]=length(sort(a15[a15[,1]>=a15[1000,1],1]))/1000 # luz
result[2,11]=length(sort(a15[a15[,2]>=a15[1000,2],2]))/1000 # ser
result[3,11]=length(sort(a15[a15[,3]>=a15[1000,3],3]))/1000 # com
result[4,11]=length(sort(a15[a15[,4]>=a15[1000,4],4]))/1000 # luz/ser
result[5,11]=length(sort(a15[a15[,5]>=a15[1000,5],5]))/1000 # luz/com
result[6,11]=length(sort(a15[a15[,6]>=a15[1000,6],6]))/1000 # ser/com
result[7,11]=length(sort(a15[a15[,7]>=a15[1000,7],7]))/1000 # luz/ser/com

return(result)
}