<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>Olá a todos<div><br></div><div>Estou escrevendo um código que calcula o lance de uma moeda jogada e a probabilidade dela acontecer em respostas de 1 para cara e 0 para coroa, num bloco do meu código fiquei com uma dúvida relacionada a vetorização. Abaixo passarei todo o código e irei destacar o que me traz a dúvida:<br><br><div>####Laboratorio 1 de R</div><div><br></div><div>###cara eh 1</div><div>###coroa eh 0</div><div>vetor_moedas <- NULL</div><div>vetor <- NULL</div><div>a <- NULL</div><div><br></div><div>moeda <- function(p){</div><div><br></div><div>if(runif(1) < p){</div><div><br></div><div> y <- 1 </div><div><br></div><div> }else{</div><div><br></div><div> y <- 0</div><div> </div><div> }</div><div> return(y)</div><div> </div><div> }</div><div> </div><div>####separando uma funcao da outra####</div><div><br></div><div>moedas <- function(numeros,probabilidade){</div><div><br></div><div> for (i in 1:numeros){</div><div><br></div><div> vetor[i] <- moeda(probabilidade)</div><div><br></div><div> }</div><div><br></div><div> return(vetor)</div><div><br></div><div> } </div><div><br></div><div><br></div><div>####separando uma funcao da outra####</div><div><br></div><div>proporcao <- function(moedas){</div><div><br></div><div> numeros_1 <- sum(vetor_moedas)</div><div><br></div><div> numeros_totais <- length(vetor_moedas) </div><div><br></div><div> proporcao_final <- numeros_1/numeros_totais</div><div><br></div><div> return(proporcao_final)</div><div><br></div><div> }</div><div><br></div><div>####separando uma funcao da outra####</div><div><br></div><div><span style="font-size: 12pt;">simulacao <- function(proporcao, moedas){</span></div><div><br></div><div>###variando o "n" e o "p"</div><div><br></div><div>for (i in 1:500){</div><div><br></div><div><span style="color: rgb(38, 114, 236); background-color: rgb(255, 255, 255);">\\Aqui começa a minha, pois como fazer com que o vetor a1 aloque como data.frame todos os dados feitos pela função moedas, como fazer com que a2 também aloque um outro data.frame que aloque informações da função moedas e assim por diante... </span></div><div><br></div><div> n <- c(10, 20, 30, 40, 50)</div><div> p <- c(0.2, 0.3, 0.4, 0.5) </div><div> </div><div>a[i] <- data.frame(moedas(n[i],p[i])) </div><div><span style="background-color: rgb(0, 0, 255);"><br></span></div><div><font color="#5133ab">\\ As vezes ele só aloca o que é produzido em moedas(n[1],p[1]) em a1 e recicla o restante. Se alguém quiser testar fiquem a vontade para ajudar!</font></div><div><br></div><div> }</div><div><br></div><div><br></div><div>return(a)</div><div><br></div><div> }</div><div><br></div><div>##############################</div><div>##############################</div><div>###Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda</div><div>###lancamento de 10 moedas - funcao moedas</div><div><br></div><div>###comeco das chamas</div><div><br></div><div>#########################################</div><div>prob <- 0.4</div><div><br></div><div>moeda(prob)</div><div><br></div><div>#########################################</div><div>n_1 <- 100</div><div><br></div><div>moedas(n_1, prob)</div><div><br></div><div>#########################################</div><div><br></div><div>vetor_moedas <- moedas(n_1, prob)</div><div><br></div><div>#########################################</div><div><br></div><div>proporcao(vetor_moedas)</div><div><br></div><div>#########################################</div><div><br></div><div>simulacao(proporcao,moedas)</div><div><br></div><div>###</div><div>###</div><div>##############################</div> </div><div><br></div><div>Um grande abraço a todos!!<br><br>Paulo Henrique de A. S. Pimenta.<div><br></div><div>Graduando em Meteorologia (Bacharelado) - IAG/USP.<br>Fone: +5511981318435.</div></div> </div></body>
</html>