----- Mensagem encaminhada ----
De: Paulo Justiniano Ribeiro Jr <[email protected]> Para: [email protected] Enviadas: Sexta-feira, 31 de Outubro de 2008 8:21:23 Assunto: Re: [R_STAT] Reamostragem Andre 1. obtenha as reamostras: reamostras <- array(sample( e, length(e)*1000, rep=T), dim=c(length( e)/2,2,1000) ) 2. faca os testes T nas 1000 reamostras TT <- apply(reamostras, 3, function(x) t.test(x[,1] , x[,2])) 3. Veja o resultado em uma das amostras: TT[[1]] 4. extraia uma quantidade que quiser como por exemplo os p-valores pvals <- sapply(TT, function(x) x$p.value) OBS: note que estao sendo usadas as opcoes default do teste t para comparacao de duas amostras dadas em t.test. (bilateral, nao pareado, var. diferentes, etc) Se quiser mudar acrescente argumentos desejados na chamada do apply(): TT <- apply(reamostras, 3, function(x) t.test(x[,1] , x[,2]), var.equal=TRUE) Paulo Justiniano Ribeiro Jr LEG (Laboratorio de Estatistica e Geoinformacao) Universidade Federal do Parana Caixa Postal 19.081 CEP 81.531-990 Curitiba, PR - Brasil Tel: (+55) 41 3361 3573 Fax: (+55) 41 3361 3141 e-mail: paulojus AT ufpr br http://www.leg. ufpr.br/~ paulojus On Fri, 31 Oct 2008, André Oliveira Souza wrote: > ola; queria fazer uma amostragem e tirar as informações. > > Olha tenho duas mostras. > > a<-c(1,2,4,5, 3,4,5,6,6, 7,2) > b<-c(5,3,4,5, 3,4,5,3,4, 3,5) > > Aplico o teste t e tenho um P-valor Po. > > Eu gostaria de juntar estes dois vetores num mesmo e fazer 1000 > reamostragem neste vetor de tamanho do vetor a e com reposição. > > e<- c(5,3,4,5,3, 4,5,3,4,3, 5,1,2,4,5, 3,4,5,6,6, 7,2) > > Depois eu queria ver aplicar o teste t(nas amostras) e ver como se > comporta estes p valores. > > att > obrigado. > > __._,_.___ Mensagens neste tópico (4) Responder (através da web) | Adicionar um novo tópico Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes Alterar configurações via web (Requer Yahoo! ID) Alterar configurações via e-mail: Alterar recebimento para lista diária de mensagens | Alterar formato para o tradicional Visite seu Grupo | Termos de uso do Yahoo! Grupos | Sair do grupo Atividade nos últimos dias * 10 Novos usuáriosVisite seu Grupo Yahoo! Mail Conecte-se ao mundo Proteção anti-spam Muito mais espaço Yahoo! Barra Instale grátis Buscar sites na web Checar seus e-mails . Yahoo! Grupos Crie seu próprio grupo A melhor forma de comunicação . __,_._,___
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
