vegas.t=function(dados,fator, nsim=1000, teste="bi") { ## funcao feita na disciplina de Introduçao ao R - 20 março 2009 ## calcula o p e retorna a matriz de permutação para o similar do teste t ## deve entrar primeiro o vetor de valores e depois o vetor de fatores ##cria o vetor de resultados resulta<-rep(NA,nsim) ## coloca o valor da diferença observada na primeira posiçao resulta[1]<-diff(tapply(dados,fator,mean)[]) ## ciclos de permutaçao for(i in 2:nsim) { resulta[i]=diff(tapply(sample(dados),fator,mean)[]) } ## cria o grafico hist(resulta) abline(v = resulta[1], col="red") abline(v = resulta[1]*-1, col="red") ###cálculo do nosso p de PATO!!! ### É signficativa a diferença?? ## define se o resultado é uni ou bicaudal if(teste=="bi") { bicaudal=sum(abs(resulta)>=abs(resulta[1])) proba=bicaudal/nsim } if(teste=="uni") { if(resulta[1]<0) { proba=sum(resulta<=resulta[1])/nsim } if(resulta[1]>=0) { proba=sum(resulta>=resulta[1])/nsim } } cat("\n\t probabilidade erro tipo I =" , proba) return(resulta) ##return(list(proba,resulta)) }