[R-br] Reamostragem

[Fabio Mathias Corrêa]

----- Mensagem encaminhada ----

De: Paulo Justiniano Ribeiro Jr <paulojus em c3sl.ufpr.br>
Para: R_STAT em yahoogrupos.com.br
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
. 
__,_._,___    


      
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110317/d8685f3d/attachment.html>