Dificuldade na utilização de um loop duplo para a criação de várias amostras
Boa tarde! Senhores, peço sua ajuda para criar um código que gere um número 'n' de amostras de tamanho 'n'. Atualmente, consigo fazer apenas uma amostra de tamanho 'n', mas gostaria de fazer de um modo prático a geração de outras 'n-1' amostras. No final gostaria de armazenar tudo numa matriz 'nxn'. Segue o código que estou usando para a geração de uma amostra: ## Tamanho da Amostra n=100 ## Parâmetros 1 a <- c(4.5,1.5) b <- c(3.0,9.0) p <- c(0.5,0.5) ## Função Acumulada Inversa fmenos1_kuma = function(x,a,b){ (1-(1-x)^(1/b))^(1/a) } ## Mistura mkuma = function(x){ (a[1]*b[1]*(x^(a[1]-1))*(1-x^a[1])^(b[1]-1))*p[1] + (a[2]*b[2]*(x^(a[2]-1))*(1-x^a[2])^(b[2]-1))*p[2] } ## Simulação da amostra u1 <- runif(n) u2 <- runif(n) x <- numeric(n) for (i in 1:n){ if(u1[i] < p[1]) x[i]= fmenos1_kuma(u2[i],a[1],b[1]) if(u1[i] >= p[1]) x[i]= fmenos1_kuma(u2[i],a[2],b[2]) } O que eu pensei em fazer foi adicionar um outro 'For' no passo da simulação da amostra, mas não consegui alimentar um objeto 'y' com os valores gerados para cada amostra. Obrigado. -- Att, Thiago Morais de Carvalho (61 - 8169-5094)
Opa, ali na sua parte de Simulação da amostra, veja que eu crio uma matriz de tamanho n*n, dai vou realizando sua simulação e guardando nela. Da pra otimizar isso usando apply, se colocar o seu codigo da simulação dentro de uma função, mas colocar outro for é relativamente mais simples. Basta a cada iteração, guardar sua simulação, uma linha de tamanho n da matriz, na matriz que você criou. Veja que agora cada linha da matriz saida é uma simulação. ## Tamanho da Amostra n=10 ## Parâmetros 1 a <- c(4.5,1.5) b <- c(3.0,9.0) p <- c(0.5,0.5) ## Função Acumulada Inversa fmenos1_kuma = function(x,a,b){ (1-(1-x)^(1/b))^(1/a) } ## Mistura mkuma = function(x){ (a[1]*b[1]*(x^(a[1]-1))*(1-x^a[1])^(b[1]-1))*p[1] + (a[2]*b[2]*(x^(a[2]-1))*(1-x^a[2])^(b[2]-1))*p[2] } ## Simulação da amostra #Crie uma matriz para receber os dados saida<-matrix(NA,n,n) saida for(j in 1:n) { u1 <- runif(n) u2 <- runif(n) x <- numeric(n) for (i in 1:n){ if(u1[i] < p[1]) x[i]= fmenos1_kuma(u2[i],a[1],b[1]) if(u1[i] >= p[1]) x[i]= fmenos1_kuma(u2[i],a[2],b[2]) } saida[,j]<-x } saida Espero ter ajudado. Se não for fazer muitas simulações isso deve resolver :) Em 19 de novembro de 2014 15:23, Thiago Touya <thiagotouya@gmail.com> escreveu:
Boa tarde!
Senhores, peço sua ajuda para criar um código que gere um número 'n' de amostras de tamanho 'n'.
Atualmente, consigo fazer apenas uma amostra de tamanho 'n', mas gostaria de fazer de um modo prático a geração de outras 'n-1' amostras. No final gostaria de armazenar tudo numa matriz 'nxn'.
Segue o código que estou usando para a geração de uma amostra:
## Tamanho da Amostra n=100
## Parâmetros 1 a <- c(4.5,1.5) b <- c(3.0,9.0) p <- c(0.5,0.5)
## Função Acumulada Inversa fmenos1_kuma = function(x,a,b){ (1-(1-x)^(1/b))^(1/a) }
## Mistura mkuma = function(x){ (a[1]*b[1]*(x^(a[1]-1))*(1-x^a[1])^(b[1]-1))*p[1] + (a[2]*b[2]*(x^(a[2]-1))*(1-x^a[2])^(b[2]-1))*p[2] }
## Simulação da amostra u1 <- runif(n) u2 <- runif(n) x <- numeric(n) for (i in 1:n){ if(u1[i] < p[1]) x[i]= fmenos1_kuma(u2[i],a[1],b[1]) if(u1[i] >= p[1]) x[i]= fmenos1_kuma(u2[i],a[2],b[2]) }
O que eu pensei em fazer foi adicionar um outro 'For' no passo da simulação da amostra, mas não consegui alimentar um objeto 'y' com os valores gerados para cada amostra.
Obrigado.
-- Att, Thiago Morais de Carvalho (61 - 8169-5094)
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056
Olá, Necessito de auxílio para abrir no RKward 3.02 (Ubuntu versão 14.04 LTS) um arquivo da base de dados da Câmara Federal, extensão csv. Parte do conteúdo segue abaixo com sep "; "... Utilizei o seguinte comando para tentar abri-lo propleg2014 <- read.csv("proposicoes.csv", sep=";", dec=",", header=TRUE) e me retornou a seguinte mensagem Warning message: In scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : EOF within quoted string Podem me ajudar a compor o comando válido que eu consiga abrir o arquivo proposicoes.csv que se encontra armazenado em minha pasta? Obrigado. DATAPRESENTACAOPROPOSICAO;DATTRANSFPROPOSICAOLEI;COD_TIPO_PARLAMENTAR;NOM_CIVIL_PARLAMENTAR;NOM_PARLAMENTAR;COD_SEXO;DES_TIPO_PARLAMENTAR;COD_PARLAMENTAR_ORIGEM;NUM_LEGISLATURA_APRESEN;NUM_SESSAO_LEGISLATIVA_APRESEN;COD_TIPO_SESSAO_LEGISLATIVA_APRESEN;DAT_INICIO_LEGISLATURA_APRESEN;DAT_FIM_LEGISLATURA_APRESEN;DAT_INICIO_SESSAO_LEGISLATIVA_APRESEN;DAT_FIM_SESSAO_LEGISLATIVA_APRESEN;DES_TIPO_SESSAO_LEGISLATIVA_APRESEN;DAT_INICIAL_RECESSO_APRESEN;DAT_FINAL_RECESSO_APRESEN;NUM_LEGISLATURA_LEI;NUM_SESSAO_LEGISLATIVA_LEI;COD_TIPO_SESSAO_LEGISLATIVA_LEI;DAT_INICIO_LEGISLATURA_LEI;DAT_FIM_LEGISLATURA_LEI;DAT_INICIO_SESSAO_LEGISLATIVA_LEI;DAT_FIM_SESSAO_LEGISLATIVA_LEI;DES_TIPO_SESSAO_LEGISLATIVA_LEI;DAT_INICIAL_RECESSO_LEI;DAT_FINAL_RECESSO_LEI;COD_PROPOSICAO_ORIGEM;ANO_PROPOSICAO;NUM_PROPOSICAO;NOM_PROPOSICAO;AREAS_TEMATICAS_APRESENTACAO;SIG_NORMA_JURIDICA;DES_TIPO_NORMA_JURIDICA;COD_TIPO_NORMA_JURIDICA_ORIGEM;SIG_UF;TEX_REGIAO_GEOGRAFICA_AUTOR;NOM_UF;NOM_PARTIDO_POLITICO;COD_PARTIDO_POLITICO;SIG_PARTIDO_POLITICO;DES_SITUACAO_PROPOSICAO;SIG_TIPO_PROPOSICAO;DES_TIPO_PROPOSICAO;DES_COMPLETA;IND_ATIVO_TIPO_PROPOSICAO;DES_ATIVO_TIPO_PROPOSICAO;DAT_VETO_TOTAL;NOM_ORGAO;IND_ATIVO_ORGAO;DES_ATIVO_ORGAO;COD_ORGAO_ORIGEM;DES_APELIDO_ORGAO;COD_TIPO_AUTOR_ORIGEM;DES_TIPO_AUTOR;NUM_LEGISLATURA_VET_TOTAL;NUM_SESSAO_LEGISLATIVA_VET_TOTAL;COD_TIPO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_INICIO_LEGISLATURA_VET_TOTAL;DAT_FIM_LEGISLATURA_VET_TOTAL;DAT_INICIO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_FIM_SESSAO_LEGISLATIVA_VET_TOTAL;DES_TIPO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_INICIAL_RECESSO_VET_TOTAL;DAT_FINAL_RECESSO_VET_TOTAL 11-03-2002;;10000;FRANCISCO DO NASCIMENTO GOMES;DR. GOMES;M;DEPUTADO;179;51;4;1;01-02-1999;31-01-2003;15-02-2002;16-12-2002;Ordinária;02-07-2002;01-08-2002;-1;-1;-1;;;;;N/A;;;46244;2002;6245;PL 6245/2002;VIAÇÃO E TRANSPORTES;N/A;N/A;-1;AM;NORTE;AMAZONAS;N/A;-1;N/A;N/A;PL;PROJETO DE LEI;PROJETO DE LEI;1;ATIVADO;;N/A;-1;N/A;-1;N/A;TipoParlamentar_10000;DEPUTADO;-1;-1;-1;;;;;N/A;; Em 19 de novembro de 2014 15:23, Thiago Touya <thiagotouya@gmail.com> escreveu:
Boa tarde!
Senhores, peço sua ajuda para criar um código que gere um número 'n' de amostras de tamanho 'n'.
Atualmente, consigo fazer apenas uma amostra de tamanho 'n', mas gostaria de fazer de um modo prático a geração de outras 'n-1' amostras. No final gostaria de armazenar tudo numa matriz 'nxn'.
Segue o código que estou usando para a geração de uma amostra:
## Tamanho da Amostra n=100
## Parâmetros 1 a <- c(4.5,1.5) b <- c(3.0,9.0) p <- c(0.5,0.5)
## Função Acumulada Inversa fmenos1_kuma = function(x,a,b){ (1-(1-x)^(1/b))^(1/a) }
## Mistura mkuma = function(x){ (a[1]*b[1]*(x^(a[1]-1))*(1-x^a[1])^(b[1]-1))*p[1] + (a[2]*b[2]*(x^(a[2]-1))*(1-x^a[2])^(b[2]-1))*p[2] }
## Simulação da amostra u1 <- runif(n) u2 <- runif(n) x <- numeric(n) for (i in 1:n){ if(u1[i] < p[1]) x[i]= fmenos1_kuma(u2[i],a[1],b[1]) if(u1[i] >= p[1]) x[i]= fmenos1_kuma(u2[i],a[2],b[2]) }
O que eu pensei em fazer foi adicionar um outro 'For' no passo da simulação da amostra, mas não consegui alimentar um objeto 'y' com os valores gerados para cada amostra.
Obrigado.
-- Att, Thiago Morais de Carvalho (61 - 8169-5094)
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Você teria o link de acesso ao arquivo? Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("") Em 20 de novembro de 2014 09:58, Antonio Armando Ulian do Lago Albuquerque < ulagoalbuquerque@gmail.com> escreveu:
Olá, Necessito de auxílio para abrir no RKward 3.02 (Ubuntu versão 14.04 LTS) um arquivo da base de dados da Câmara Federal, extensão csv. Parte do conteúdo segue abaixo com sep "; "... Utilizei o seguinte comando para tentar abri-lo propleg2014 <- read.csv("proposicoes.csv", sep=";", dec=",", header=TRUE) e me retornou a seguinte mensagem Warning message: In scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : EOF within quoted string Podem me ajudar a compor o comando válido que eu consiga abrir o arquivo proposicoes.csv que se encontra armazenado em minha pasta? Obrigado.
DATAPRESENTACAOPROPOSICAO;DATTRANSFPROPOSICAOLEI;COD_TIPO_PARLAMENTAR;NOM_CIVIL_PARLAMENTAR;NOM_PARLAMENTAR;COD_SEXO;DES_TIPO_PARLAMENTAR;COD_PARLAMENTAR_ORIGEM;NUM_LEGISLATURA_APRESEN;NUM_SESSAO_LEGISLATIVA_APRESEN;COD_TIPO_SESSAO_LEGISLATIVA_APRESEN;DAT_INICIO_LEGISLATURA_APRESEN;DAT_FIM_LEGISLATURA_APRESEN;DAT_INICIO_SESSAO_LEGISLATIVA_APRESEN;DAT_FIM_SESSAO_LEGISLATIVA_APRESEN;DES_TIPO_SESSAO_LEGISLATIVA_APRESEN;DAT_INICIAL_RECESSO_APRESEN;DAT_FINAL_RECESSO_APRESEN;NUM_LEGISLATURA_LEI;NUM_SESSAO_LEGISLATIVA_LEI;COD_TIPO_SESSAO_LEGISLATIVA_LEI;DAT_INICIO_LEGISLATURA_LEI;DAT_FIM_LEGISLATURA_LEI;DAT_INICIO_SESSAO_LEGISLATIVA_LEI;DAT_FIM_SESSAO_LEGISLATIVA_LEI;DES_TIPO_SESSAO_LEGISLATIVA_LEI;DAT_INICIAL_RECESSO_LEI;DAT_FINAL_RECESSO_LEI;COD_PROPOSICAO_ORIGEM;ANO_PROPOSICAO;NUM_PROPOSICAO;NOM_PROPOSICAO;AREAS_TEMATICAS_APRESENTACAO;SIG_NORMA_JURIDICA;DES_TIPO_NORMA_JURIDICA;COD_TIPO_NORMA_JURIDICA_ORIGEM;SIG_UF;TEX_REGIAO_GEOGRAFICA_AUTOR;NOM_UF;NOM_PARTIDO_POLITICO;COD_PARTIDO_POLITICO;SIG_PARTIDO_POLITICO;DES_SITUACAO_PROPOSICAO;SIG_TIPO_PROPOSICAO;DES_TIPO_PROPOSICAO;DES_COMPLETA;IND_ATIVO_TIPO_PROPOSICAO;DES_ATIVO_TIPO_PROPOSICAO;DAT_VETO_TOTAL;NOM_ORGAO;IND_ATIVO_ORGAO;DES_ATIVO_ORGAO;COD_ORGAO_ORIGEM;DES_APELIDO_ORGAO;COD_TIPO_AUTOR_ORIGEM;DES_TIPO_AUTOR;NUM_LEGISLATURA_VET_TOTAL;NUM_SESSAO_LEGISLATIVA_VET_TOTAL;COD_TIPO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_INICIO_LEGISLATURA_VET_TOTAL;DAT_FIM_LEGISLATURA_VET_TOTAL;DAT_INICIO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_FIM_SESSAO_LEGISLATIVA_VET_TOTAL;DES_TIPO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_INICIAL_RECESSO_VET_TOTAL;DAT_FINAL_RECESSO_VET_TOTAL 11-03-2002;;10000;FRANCISCO DO NASCIMENTO GOMES;DR. GOMES;M;DEPUTADO;179;51;4;1;01-02-1999;31-01-2003;15-02-2002;16-12-2002;Ordinária;02-07-2002;01-08-2002;-1;-1;-1;;;;;N/A;;;46244;2002;6245;PL 6245/2002;VIAÇÃO E TRANSPORTES;N/A;N/A;-1;AM;NORTE;AMAZONAS;N/A;-1;N/A;N/A;PL;PROJETO DE LEI;PROJETO DE LEI;1;ATIVADO;;N/A;-1;N/A;-1;N/A;TipoParlamentar_10000;DEPUTADO;-1;-1;-1;;;;;N/A;;
Em 19 de novembro de 2014 15:23, Thiago Touya <thiagotouya@gmail.com> escreveu:
Boa tarde!
Senhores, peço sua ajuda para criar um código que gere um número 'n' de amostras de tamanho 'n'.
Atualmente, consigo fazer apenas uma amostra de tamanho 'n', mas gostaria de fazer de um modo prático a geração de outras 'n-1' amostras. No final gostaria de armazenar tudo numa matriz 'nxn'.
Segue o código que estou usando para a geração de uma amostra:
## Tamanho da Amostra n=100
## Parâmetros 1 a <- c(4.5,1.5) b <- c(3.0,9.0) p <- c(0.5,0.5)
## Função Acumulada Inversa fmenos1_kuma = function(x,a,b){ (1-(1-x)^(1/b))^(1/a) }
## Mistura mkuma = function(x){ (a[1]*b[1]*(x^(a[1]-1))*(1-x^a[1])^(b[1]-1))*p[1] + (a[2]*b[2]*(x^(a[2]-1))*(1-x^a[2])^(b[2]-1))*p[2] }
## Simulação da amostra u1 <- runif(n) u2 <- runif(n) x <- numeric(n) for (i in 1:n){ if(u1[i] < p[1]) x[i]= fmenos1_kuma(u2[i],a[1],b[1]) if(u1[i] >= p[1]) x[i]= fmenos1_kuma(u2[i],a[2],b[2]) }
O que eu pensei em fazer foi adicionar um outro 'For' no passo da simulação da amostra, mas não consegui alimentar um objeto 'y' com os valores gerados para cada amostra.
Obrigado.
-- Att, Thiago Morais de Carvalho (61 - 8169-5094)
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
vou compartilhar em dropbox com vc ok! Em 20 de novembro de 2014 10:05, David Feitosa <contato@davidfeitosa.com> escreveu:
Você teria o link de acesso ao arquivo?
Atenciosamente,
David Feitosa
(\_(\ (=°;°) (("")("")
Em 20 de novembro de 2014 09:58, Antonio Armando Ulian do Lago Albuquerque <ulagoalbuquerque@gmail.com> escreveu:
Olá,
Necessito de auxílio para abrir no RKward 3.02 (Ubuntu versão 14.04 LTS) um arquivo da base de dados da Câmara Federal, extensão csv. Parte do conteúdo segue abaixo com sep "; "... Utilizei o seguinte comando para tentar abri-lo propleg2014 <- read.csv("proposicoes.csv", sep=";", dec=",", header=TRUE) e me retornou a seguinte mensagem Warning message: In scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : EOF within quoted string Podem me ajudar a compor o comando válido que eu consiga abrir o arquivo proposicoes.csv que se encontra armazenado em minha pasta? Obrigado.
DATAPRESENTACAOPROPOSICAO;DATTRANSFPROPOSICAOLEI;COD_TIPO_PARLAMENTAR;NOM_CIVIL_PARLAMENTAR;NOM_PARLAMENTAR;COD_SEXO;DES_TIPO_PARLAMENTAR;COD_PARLAMENTAR_ORIGEM;NUM_LEGISLATURA_APRESEN;NUM_SESSAO_LEGISLATIVA_APRESEN;COD_TIPO_SESSAO_LEGISLATIVA_APRESEN;DAT_INICIO_LEGISLATURA_APRESEN;DAT_FIM_LEGISLATURA_APRESEN;DAT_INICIO_SESSAO_LEGISLATIVA_APRESEN;DAT_FIM_SESSAO_LEGISLATIVA_APRESEN;DES_TIPO_SESSAO_LEGISLATIVA_APRESEN;DAT_INICIAL_RECESSO_APRESEN;DAT_FINAL_RECESSO_APRESEN;NUM_LEGISLATURA_LEI;NUM_SESSAO_LEGISLATIVA_LEI;COD_TIPO_SESSAO_LEGISLATIVA_LEI;DAT_INICIO_LEGISLATURA_LEI;DAT_FIM_LEGISLATURA_LEI;DAT_INICIO_SESSAO_LEGISLATIVA_LEI;DAT_FIM_SESSAO_LEGISLATIVA_LEI;DES_TIPO_SESSAO_LEGISLATIVA_LEI;DAT_INICIAL_RECESSO_LEI;DAT_FINAL_RECESSO_LEI;COD_PROPOSICAO_ORIGEM;ANO_PROPOSICAO;NUM_PROPOSICAO;NOM_PROPOSICAO;AREAS_TEMATICAS_APRESENTACAO;SIG_NORMA_JURIDICA;DES_TIPO_NORMA_JURIDICA;COD_TIPO_NORMA_JURIDICA_ORIGEM;SIG_UF;TEX_REGIAO_GEOGRAFICA_AUTOR;NOM_UF;NOM_PARTIDO_POLITICO;COD_PARTIDO_POLITICO;SIG_PARTIDO_POLITICO;DES_SITUACAO_PROPOSICAO;SIG_TIPO_PROPOSICAO;DES_TIPO_PROPOSICAO;DES_COMPLETA;IND_ATIVO_TIPO_PROPOSICAO;DES_ATIVO_TIPO_PROPOSICAO;DAT_VETO_TOTAL;NOM_ORGAO;IND_ATIVO_ORGAO;DES_ATIVO_ORGAO;COD_ORGAO_ORIGEM;DES_APELIDO_ORGAO;COD_TIPO_AUTOR_ORIGEM;DES_TIPO_AUTOR;NUM_LEGISLATURA_VET_TOTAL;NUM_SESSAO_LEGISLATIVA_VET_TOTAL;COD_TIPO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_INICIO_LEGISLATURA_VET_TOTAL;DAT_FIM_LEGISLATURA_VET_TOTAL;DAT_INICIO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_FIM_SESSAO_LEGISLATIVA_VET_TOTAL;DES_TIPO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_INICIAL_RECESSO_VET_TOTAL;DAT_FINAL_RECESSO_VET_TOTAL 11-03-2002;;10000;FRANCISCO DO NASCIMENTO GOMES;DR. GOMES;M;DEPUTADO;179;51;4;1;01-02-1999;31-01-2003;15-02-2002;16-12-2002;Ordinária;02-07-2002;01-08-2002;-1;-1;-1;;;;;N/A;;;46244;2002;6245;PL 6245/2002;VIAÇÃO E TRANSPORTES;N/A;N/A;-1;AM;NORTE;AMAZONAS;N/A;-1;N/A;N/A;PL;PROJETO DE LEI;PROJETO DE LEI;1;ATIVADO;;N/A;-1;N/A;-1;N/A;TipoParlamentar_10000;DEPUTADO;-1;-1;-1;;;;;N/A;;
Em 19 de novembro de 2014 15:23, Thiago Touya <thiagotouya@gmail.com> escreveu:
Boa tarde!
Senhores, peço sua ajuda para criar um código que gere um número 'n' de amostras de tamanho 'n'.
Atualmente, consigo fazer apenas uma amostra de tamanho 'n', mas gostaria de fazer de um modo prático a geração de outras 'n-1' amostras. No final gostaria de armazenar tudo numa matriz 'nxn'.
Segue o código que estou usando para a geração de uma amostra:
## Tamanho da Amostra n=100
## Parâmetros 1 a <- c(4.5,1.5) b <- c(3.0,9.0) p <- c(0.5,0.5)
## Função Acumulada Inversa fmenos1_kuma = function(x,a,b){ (1-(1-x)^(1/b))^(1/a) }
## Mistura mkuma = function(x){ (a[1]*b[1]*(x^(a[1]-1))*(1-x^a[1])^(b[1]-1))*p[1] + (a[2]*b[2]*(x^(a[2]-1))*(1-x^a[2])^(b[2]-1))*p[2] }
## Simulação da amostra u1 <- runif(n) u2 <- runif(n) x <- numeric(n) for (i in 1:n){ if(u1[i] < p[1]) x[i]= fmenos1_kuma(u2[i],a[1],b[1]) if(u1[i] >= p[1]) x[i]= fmenos1_kuma(u2[i],a[2],b[2]) }
O que eu pensei em fazer foi adicionar um outro 'For' no passo da simulação da amostra, mas não consegui alimentar um objeto 'y' com os valores gerados para cada amostra.
Obrigado.
-- Att, Thiago Morais de Carvalho (61 - 8169-5094)
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Não tenho como acessar o Dropbox de onde estou agora. Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("") Em 20 de novembro de 2014 10:17, Antonio Armando Ulian do Lago Albuquerque < ulagoalbuquerque@gmail.com> escreveu:
vou compartilhar em dropbox com vc ok!
Em 20 de novembro de 2014 10:05, David Feitosa <contato@davidfeitosa.com> escreveu:
Você teria o link de acesso ao arquivo?
Atenciosamente,
David Feitosa
(\_(\ (=°;°) (("")("")
Em 20 de novembro de 2014 09:58, Antonio Armando Ulian do Lago Albuquerque <ulagoalbuquerque@gmail.com> escreveu:
Olá,
Necessito de auxílio para abrir no RKward 3.02 (Ubuntu versão 14.04 LTS) um arquivo da base de dados da Câmara Federal, extensão csv. Parte do conteúdo segue abaixo com sep "; "... Utilizei o seguinte comando para tentar abri-lo propleg2014 <- read.csv("proposicoes.csv", sep=";", dec=",", header=TRUE) e me retornou a seguinte mensagem Warning message: In scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : EOF within quoted string Podem me ajudar a compor o comando válido que eu consiga abrir o arquivo proposicoes.csv que se encontra armazenado em minha pasta? Obrigado.
DATAPRESENTACAOPROPOSICAO;DATTRANSFPROPOSICAOLEI;COD_TIPO_PARLAMENTAR;NOM_CIVIL_PARLAMENTAR;NOM_PARLAMENTAR;COD_SEXO;DES_TIPO_PARLAMENTAR;COD_PARLAMENTAR_ORIGEM;NUM_LEGISLATURA_APRESEN;NUM_SESSAO_LEGISLATIVA_APRESEN;COD_TIPO_SESSAO_LEGISLATIVA_APRESEN;DAT_INICIO_LEGISLATURA_APRESEN;DAT_FIM_LEGISLATURA_APRESEN;DAT_INICIO_SESSAO_LEGISLATIVA_APRESEN;DAT_FIM_SESSAO_LEGISLATIVA_APRESEN;DES_TIPO_SESSAO_LEGISLATIVA_APRESEN;DAT_INICIAL_RECESSO_APRESEN;DAT_FINAL_RECESSO_APRESEN;NUM_LEGISLATURA_LEI;NUM_SESSAO_LEGISLATIVA_LEI;COD_TIPO_SESSAO_LEGISLATIVA_LEI;DAT_INICIO_LEGISLATURA_LEI;DAT_FIM_LEGISLATURA_LEI;DAT_INICIO_SESSAO_LEGISLATIVA_LEI;DAT_FIM_SESSAO_LEGISLATIVA_LEI;DES_TIPO_SESSAO_LEGISLATIVA_LEI;DAT_INICIAL_RECESSO_LEI;DAT_FINAL_RECESSO_LEI;COD_PROPOSICAO_ORIGEM;ANO_PROPOSICAO;NUM_PROPOSICAO;NOM_PROPOSICAO;AREAS_TEMATICAS_APRESENTACAO;SIG_NORMA_JURIDICA;DES_TIPO_NORMA_JURIDICA;COD_TIPO_NORMA_JURIDICA_ORIGEM;SIG_UF;TEX_REGIAO_GEOGRAFICA_AUTOR;NOM_UF;NOM_PARTIDO_POLITICO;COD_PARTIDO_POLITICO;SIG_PARTIDO_POLITICO;DES_SITUACAO_PROPOSICAO;SIG_TIPO_PROPOSICAO;DES_TIPO_PROPOSICAO;DES_COMPLETA;IND_ATIVO_TIPO_PROPOSICAO;DES_ATIVO_TIPO_PROPOSICAO;DAT_VETO_TOTAL;NOM_ORGAO;IND_ATIVO_ORGAO;DES_ATIVO_ORGAO;COD_ORGAO_ORIGEM;DES_APELIDO_ORGAO;COD_TIPO_AUTOR_ORIGEM;DES_TIPO_AUTOR;NUM_LEGISLATURA_VET_TOTAL;NUM_SESSAO_LEGISLATIVA_VET_TOTAL;COD_TIPO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_INICIO_LEGISLATURA_VET_TOTAL;DAT_FIM_LEGISLATURA_VET_TOTAL;DAT_INICIO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_FIM_SESSAO_LEGISLATIVA_VET_TOTAL;DES_TIPO_SESSAO_LEGISLATIVA_VET_TOTAL;DAT_INICIAL_RECESSO_VET_TOTAL;DAT_FINAL_RECESSO_VET_TOTAL 11-03-2002;;10000;FRANCISCO DO NASCIMENTO GOMES;DR. GOMES;M;DEPUTADO;179;51;4;1;01-02-1999;31-01-2003;15-02-2002;16-12-2002;Ordinária;02-07-2002;01-08-2002;-1;-1;-1;;;;;N/A;;;46244;2002;6245;PL 6245/2002;VIAÇÃO E TRANSPORTES;N/A;N/A;-1;AM;NORTE;AMAZONAS;N/A;-1;N/A;N/A;PL;PROJETO DE LEI;PROJETO DE LEI;1;ATIVADO;;N/A;-1;N/A;-1;N/A;TipoParlamentar_10000;DEPUTADO;-1;-1;-1;;;;;N/A;;
Em 19 de novembro de 2014 15:23, Thiago Touya <thiagotouya@gmail.com> escreveu:
Boa tarde!
Senhores, peço sua ajuda para criar um código que gere um número 'n' de amostras de tamanho 'n'.
Atualmente, consigo fazer apenas uma amostra de tamanho 'n', mas gostaria de fazer de um modo prático a geração de outras 'n-1' amostras. No final gostaria de armazenar tudo numa matriz 'nxn'.
Segue o código que estou usando para a geração de uma amostra:
## Tamanho da Amostra n=100
## Parâmetros 1 a <- c(4.5,1.5) b <- c(3.0,9.0) p <- c(0.5,0.5)
## Função Acumulada Inversa fmenos1_kuma = function(x,a,b){ (1-(1-x)^(1/b))^(1/a) }
## Mistura mkuma = function(x){ (a[1]*b[1]*(x^(a[1]-1))*(1-x^a[1])^(b[1]-1))*p[1] + (a[2]*b[2]*(x^(a[2]-1))*(1-x^a[2])^(b[2]-1))*p[2] }
## Simulação da amostra u1 <- runif(n) u2 <- runif(n) x <- numeric(n) for (i in 1:n){ if(u1[i] < p[1]) x[i]= fmenos1_kuma(u2[i],a[1],b[1]) if(u1[i] >= p[1]) x[i]= fmenos1_kuma(u2[i],a[2],b[2]) }
O que eu pensei em fazer foi adicionar um outro 'For' no passo da simulação da amostra, mas não consegui alimentar um objeto 'y' com os valores gerados para cada amostra.
Obrigado.
-- Att, Thiago Morais de Carvalho (61 - 8169-5094)
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
participantes (4)
-
Antonio Armando Ulian do Lago Albuquerque -
Augusto Ribas -
David Feitosa -
Thiago Touya