Boa tarde pessoal,

estou trabalhando com a função disponível no link <http://morotalab.org/Mrode2005/relmat/quass.txt>

e gostaria de aumentar a velocidade de execução da funcão.

Qualquer opinião é bem vinda.

Segue abaixo um CMR:

sire <- rep(1:50,each=50)

dam <- rep(1:50,50)

s <- c(rep(0,50),sire)

d = c(rep(0,50),dam) 

system.time(quass(p,m)) 

Desde já obrigado