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

•  10
  Novos usuários

Visite seu Grupo