[R-br] bootstrap

[Pedro Emmanuel Alvarenga Americano do Brasil]

Maria,

Um pedaço de script que talvez te ajude.

data <- as.data.frame(cbind(c(20:60),c(30:70),c(40:80)))
data
B = 500
medias <- matrix(NA,ncol=3,nrow=B)
for(i in 1:B){
  medias[i,] <- colMeans(data[sample(1:nrow(data),nrow(data)-1,T),]) #
bootstrap com todas os dados -1 uma linha
  }
medias

Dr. Pedro Emmanuel A. A. do Brasil
Curriculum Lattes:  http://lattes.cnpq.br/6597654894290806
Instituto de Pesquisa Clínica Evandro Chagas
Fundação Oswaldo Cruz
Rio de Janeiro - Brasil
Av. Brasil 4365,
CEP 21040-360,
Tel 55 21 3865-9648
email: pedro.brasil em ipec.fiocruz.br
email: emmanuel.brasil em gmail.com

---Apoio aos softwares livres
www.zotero.org - gerenciamento de referências bibliográficas.
www.broffice.org ou www.libreoffice.org - textos, planilhas ou
apresentações.
www.epidata.dk - entrada de dados.
www.r-project.org - análise de dados.
www.ubuntu.com - sistema operacional



Em 22 de abril de 2012 23:20, Maria Papa <mceliamat em yahoo.com.br> escreveu:

> Olá Pessoal, boa noite
>
> Estou tentando, sem sucesso, realizar uma reamostragem bootstrap para
> médias de 3 amostras estão organizados em uma matriz na qual, que sendo
> cada coluna representa uma amostra. Porém, a cada reamostragem, eu
> preciso que seja feita para as três amostras, ou seja, quando um dado for
> selecionado para formar a amostra bootstrap, deve ser retirada a linha
> toda.  O resultado final deve ser 3  médias reamostradas, que foram
> obtidas em um único processo de reamostragem. Alguém sabe como fazer isso
> no R? Eu consegui somente quando eu realizo 3 procedimentos de
> reamostragem distintos.
>
> Desde já agradeço a todos
>
> célia
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120423/17c020fe/attachment.html>