Ir ao conteúdo
  • Cadastre-se
felipecastroo

Simulando o efeito TCL

Recommended Posts

Fala galera, boa noite !

Olha só, estou precisando de uma ajuda enorme de vocês...

Estou com um código no R e não estou conseguindo trocar a amostra que o exemplo utilizou pela minha..

O código é o seguinte :

tcl.exp=function(n, N=200, titulo=" ", yl=c(0, .4)) 

{ ## início da função – tcl.exp
medias=numeric(N)
for (i in 1:N) medias[i]= mean(rexp(n,1/3))
hist(medias, xlim=c(-1,10), ylim=yl, freq=F, main=titulo)
x=seq(-1,10, .02)
points(x, dnorm(x, 3, 3*sqrt(1/n) ), type="l", lwd=3)
} ## fim da função

graphics.off()
par(mfrow=c(3,3), mai=c(.3,.4,.1,.1))
tcl.exp(1,titulo="n=1")
tcl.exp(2,titulo="n=2")
tcl.exp(3,titulo="n=3")
tcl.exp(4,titulo="n=4")
tcl.exp(5,titulo="n=5")
tcl.exp(6,titulo="n=6")
tcl.exp(10,titulo="n=10",yl=c(0,.6))
tcl.exp(15,titulo="n=15",yl=c(0,.6))
tcl.exp(20,titulo="n=20",yl=c(0,.6))

A primeira parte corre de 1 a 200 ( número de repetições pra cada valor de n ( 1,2,3,4,5,10,15,20) e a segunda mostra os histogramas..
O meu problema é que não estou trocar a amostra que o exemplo utilizou pela minha..
Estou com a amostra aqui no pc, e já passei pro R da seguinte forma:

banco <- read.csv("E:/Desktop/3 PERIODO/dados.csv",sep="\t",dec=",",header = TRUE)
set.seed(112041461)
amostra <- banco[runif(200, 1, nrow(banco)),]

É uma amostra aleatória de um bando de dados passado pelo professor !

Mas é isso, agradeço qualquer ajuda e fico no aguardo !

Abraços !

Compartilhar este post


Link para o post
Compartilhar em outros sites

Crie uma conta ou entre para comentar

Você precisar ser um membro para fazer um comentário

Criar uma conta

Crie uma nova conta em nossa comunidade. É fácil!

Crie uma nova conta

Entrar

Já tem uma conta? Faça o login.

Entrar agora





Sobre o Clube do Hardware

No ar desde 1996, o Clube do Hardware é uma das maiores, mais antigas e mais respeitadas publicações sobre tecnologia do Brasil. Leia mais

Direitos autorais

Não permitimos a cópia ou reprodução do conteúdo do nosso site, fórum, newsletters e redes sociais, mesmo citando-se a fonte. Leia mais

×