[R-br] Amostragem em um vetor

[Fernando Souza]

Acredito que isso deva te ajudar. Não entendi muito bem o que quis dizer
com  " para cada  tam repetir  nam (vezes)" , por isso o código abaixo não
faz nenhuma repetição dos valores amostrados , porém acredito que é fácil
de ajustar.

populacao <- rchisq(100000, df = 10)
tamanho <- c(6,14,16,18,20) # diferentes tamanhos para amostrar

amostra <- function(x,n){
    a <- list()

    for(i in 1:length(n)){

        a[[i]] <- sample(x,size = n[i])


        }
   return(a)

}
amostrados <- amostra(populacao,tamanho)

#install.packages("plyr")
library(plyr)

resultado <- ldply (amostrados,function(x){
   media <- mean(x)
   desvio <- sd(x)

   resposta <- data.frame(media,desvio)
   return(resposta)

})

##você pode converter o dataframe para matrix
as.matrix(resultado)


Em 26 de setembro de 2017 22:42, Andre Oliveira via R-br <
r-br em listas.c3sl.ufpr.br> escreveu:

> Boa noite,
> preciso de amostrar em um vetor de nome populacao amostras de tamanho n,
> com n variando de tamanhos  [2; tam] e para cada  tam repetir  nam
> (vezes) e assim pegar a estimativas de media e dos desvios -padrão das
> amostras.
>
>
> Para cada amostra obtida extrair desvio padrão, média e guardar em uma
> matriz com vazia com NAs!  Ao final plotar a média dos desvios e a media
> das medias em função de cada n.
>
> Tentei fazer,  mas sem sucesso! Caso alguém pude ajudar ficarei grato!
>
>
> ############################################################
> #########################
>
> n=10
> populacao <- rchisq(100000, df = n)
> mu=n
> s=sqrt(2*n)
> mu
> s
>
> par(mfrow = c(1, 3))
> hist(populacao, prob=T, main="Y ~ Quiquadrado(GL=10)",xlab=" ",
> yla="Densidade",col="limegreen")
>
> amostragem<-function(populacao, tam, nam)
>   {
>    SD<-NULL
>    media<-NULL
>     X <-matrix(NA, tam, 2)
>     {
>     for(i in 2:nam)
>     amostra<-sample(populacao,tam)
>     SD[i]=sd(amostra)
>     media[i]<-mean(amostra)
>    X[i,] <- c(mean(SD), mean(media))
>     }
> print(cbind(mean(SD),mean(media))
> print(X)
> plot.ts(SD, pch=19, type = "p", xlab="Número de folhas amostrada",
> ylab="Estimativa da Variabiliade", lty=2,xaxp=c(0,150,10))
> plot.ts(media, pch=19, type = "p", xlab="Número de folhas amostrada",
> ylab="Estimativa da média", lty=2,xaxp=c(0,150,10))
> }
>
>
> amostragem(populacao, 150, 50000)
>  ###########################################################
> ################################
>
> obg
>
> _______________________________________________
> 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.
>



-- 
=========================================
Fernando Souza
Zootecnista, DSc. Produção e Alimentação Animal
Celular: (31)99796-8781 (Vivo)
E-mail:nandodesouza em gmail.com <e-mail%3Anandodesouza em gmail.com>
Lattes: http://lattes.cnpq.br/6519538815038307
Blog: https://producaoanimalcomr.wordpress.com/
==========================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170927/d9654d16/attachment.html>