[R-br] Dúvida em looping e data.frame
Paulo Henrique Pimenta
paulopimenta6 em hotmail.com
Quarta Janeiro 27 04:33:10 BRST 2016
Olá a todos
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:
####Laboratorio 1 de R
###cara eh 1###coroa eh 0vetor_moedas <- NULLvetor <- NULLa <- NULL
moeda <- function(p){
if(runif(1) < p){
y <- 1
}else{
y <- 0 } return(y) } ####separando uma funcao da outra####
moedas <- function(numeros,probabilidade){
for (i in 1:numeros){
vetor[i] <- moeda(probabilidade)
}
return(vetor)
}
####separando uma funcao da outra####
proporcao <- function(moedas){
numeros_1 <- sum(vetor_moedas)
numeros_totais <- length(vetor_moedas)
proporcao_final <- numeros_1/numeros_totais
return(proporcao_final)
}
####separando uma funcao da outra####
simulacao <- function(proporcao, moedas){
###variando o "n" e o "p"
for (i in 1:500){
\\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...
n <- c(10, 20, 30, 40, 50) p <- c(0.2, 0.3, 0.4, 0.5) a[i] <- data.frame(moedas(n[i],p[i]))
\\ 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!
}
return(a)
}
###############################################################Probabilidade de ser cara ou 1 eh 0.5 - funcao moeda###lancamento de 10 moedas - funcao moedas
###comeco das chamas
#########################################prob <- 0.4
moeda(prob)
#########################################n_1 <- 100
moedas(n_1, prob)
#########################################
vetor_moedas <- moedas(n_1, prob)
#########################################
proporcao(vetor_moedas)
#########################################
simulacao(proporcao,moedas)
####################################
Um grande abraço a todos!!
Paulo Henrique de A. S. Pimenta.
Graduando em Meteorologia (Bacharelado) - IAG/USP.
Fone: +5511981318435.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160127/683e3d06/attachment.html>
Mais detalhes sobre a lista de discussão R-br