reamostragem automatizada

Boa tarde caros colegas da r-br, Gostaria da opinião dos senhores em um processo de automatização de geração de amostras (reamostragem) que iniciei. A ideia é gerar n amostras de tamanhos diferentes e verificar se as médias são estatísticamente iguais. Para fins de exemplificação, é como se eu fosse no campo e medisse inicialmente altura de 5 árvores, voltasse nessa área e medisse 3, posteriormente 2 e assim por diante para testar minha suficiência amostral. ## Segue um breve RCMR #### Reamostragem #### h <- c(18.79, 12.34, 14.11, 15.24, 13.82, 14.54, 13.77, 13.99, 10.00, 16.80, 10.49, 20.00, 14.19, 14.09, 14.65, 12.22, 16.21, 13.90, 16.29, 15.39); parc <- rep(1:4, each=5); arv <- rep(1:5, 4); dados <- data.frame(parc, arv, h); dados[1:6,] rm(arv, parc, h); dados_sorteio <- dados[unlist( with(dados, tapply(seq(nrow(dados)), parc, sample, size=3) ) ),] dados_sorteio # neste passo eu gero um novo dataframe reduzindo o tamanho da minha amostra, porém fico restrito a uma única comparação e a ideia é gerar esse procedimento n vezes e comparar as médias e também a estabilização do coeficiente de variação (curva de CV x n). Agradeço pela gentileza e contribuição de todos. Abraços Samuel -- Prof. Samuel P C Carvalho Manejo Florestal e Métodos Quantitativos Faculdade de Engenharia Florestal - UFMT (065)3615-8679 / Ramal 218 www.ufmt.br/fenf

Samuel, boa tarde! Segue uma sugestão, ### <code r> dados <- data.frame(h=c(18.79, 12.34, 14.11, 15.24, 13.82, 14.54, 13.77, 13.99, 10.00, 16.80, 10.49, 20.00, 14.19, 14.09, 14.65, 12.22, 16.21, 13.90, 16.29, 15.39), parc=rep(1:4, each=5), arv=rep(1:5, 4)) {set.seed(765); sorteio <- lapply(10:1, function(x) unlist(sample(seq(dados[,1]),x)))} sorteio datasets <- lapply(sorteio, function(x) dados[unlist(x),]); datasets medias <- sapply(sorteio, function(x) c(colMeans(dados[unlist(x),]), n=length(x))); medias soma <- sapply(sorteio, function(x) c(colSums(dados[unlist(x),]), n=length(x))); soma ### </code> Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
participantes (2)
-
Samuel Carvalho
-
Éder Comunello