[R-br] Amostragem em um vetor
Andre Oliveira
andreolsouza em yahoo.com.br
Terça Setembro 26 22:42:56 -03 2017
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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170927/d8fbe6a4/attachment.html>
Mais detalhes sobre a lista de discussão R-br