Estatística no Programa R
Gostaria de reagir a esta mensagem? Crie uma conta em poucos cliques ou inicie sessão para continuar.

Simulando o efeito TCL

Ir para baixo

Simulando o efeito TCL Empty Simulando o efeito TCL

Mensagem  felipecastro Qui Jul 04, 2013 10:15 pm

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 eh 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 !

felipecastro

Mensagens : 1
Data de inscrição : 04/07/2013

Ir para o topo Ir para baixo

Ir para o topo


 
Permissões neste sub-fórum
Não podes responder a tópicos