Maio 2011 - R-br - DINF UFPR

Re: [R-br] Microdados PNAD
by Elias T. Krainski 09 Mai '11

09 Mai '11

Régis, A função read.fwf() é idicada para ler microdados se vc tem MUITA memoria RAM ou se seu arquivo é pequeno (o que não é o seu caso). Com essa função vc lê todas as colunas dos dados, o que pode "estourar" a memória RAM. Eu acho apropriado ler apenas as colunas e linhas que interessam. Por exemplo, se quero estimar a proporção de mulheres que trabalham em segurança, basta eu ler a coluna sexo (e peso amostral) e apenas as linhas referentes aos indivíduos que trabalham com segurança. Eu uso uma função disponível em http://leg.ufpr.br/~elias/ensino/ce223/aula18.R No exemplo, eu considero dados da pnad 2006. Att. Elias T. Krainski

2 1

Redes neurais
by Marcelo Claro de Souza 09 Mai '11

09 Mai '11

Alguém teria um tutorial bem básico sobre redes neurais no R? Muito obrigado.

4 3

Rstudio Muitos dados
by eder＠leg.ufpr.br 09 Mai '11

09 Mai '11

Pessoal gostaria de comportilhar com vocês, Primeiramente gostaria de dizer que o Rstudio tem uma interface muito amigavel, porem estou usando o RStudio Desktop 0.93.89 - Windows XP/Vista/7 em um computador I7, 4gb Ram, 3Tb hd, Windows 7 64 bit e R 2.12.0, Quando gero por exemplo uma matrix de 500 x 500 e mando executa-la ele trava, travando toda a maquina, quando um procedimento destes é feito no RGui tudo roda sem problemas algum, desta forma gostaria de saber se alguem de vocês já teve problema similar com o Rstudio no Windows. Outra coisa que percebi é que quando se gera um objeto por exemplo uma matrix 500x500, avaliando o Rgui com o R executado no Rstudio o Rstudio aloca muito mais memória para o processo, acho isso muito estranho. Att

1 0

Microdados PNAD
by Régis B. de Oliveira 09 Mai '11

09 Mai '11

Boa Tarde a todos. Sou novato na lista e no uso do R. Trabalho com microdados do IBGE, mais especificamente os dados da Pesquisa Nacional por Amostra de Domicílios (PNAD). Alguem da lista tem experiência na manipulação destes dados com o R. Achei muito confusa a explicação que contém no site do IBGE (leitura em R). Preciso de um "empurrãozinho" inicial. Trabalho como o SAS e pretendo migrar. A dificuldade inicial é justamente ler os microdados. No SAS eu costumo ler todo o arquivo txt e salvar na extesão .sas7bdat. A partir daí, faço minhas análises (basicamente regressões e tabelas de frequencias). É isso. Desde já agradeço, -- Régis

5 5

Erros no output do Vim!
by Ivan Bezerra Allaman 08 Mai '11

08 Mai '11

Bom dia senhores!! Ao utilizar o Vim, estou encontrando alguns erros no output. Primeiro, é necessário eu executar várias vezes o comando ";l " para aparecer o output desejado. Segundo, o output está aparecendo pela metade. Por exemplo, quando rodo uma base de dados com 400 linhas por exemplo, na primeira execução só aparece as 50 primeiras, na segunda execução aparece os 150 e assim por diante até aparecer toda a base. Alguém já se deparou com este problema? Qual a solução? Só para complementar, já instalei a última atualização do "screen.vba". Abraços. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Reenvidando o email multcompBoxplot
by Edson Lira 08 Mai '11

08 Mai '11

Boa noite a todos, A rotina abaixo: require(multcompView) difat<-c(3.078444,3.567257,4.023913,3.7525714) names(difat) <- c("A-B", "A-C", "B-C", "A-D") (mcT4 <- multcompTs(difat)) mcT4.1 <- plot(mcT4, label.groups=0.05) (mcL4 <- multcompLetters(difat, Letters=LETTERS)) mcL4.1 <- plot(mcL4, label.groups=0.05) Mostra a seguinte mensagem: Erro em plot.window(...) : valores finitos são necessários para 'xlim' Além disso: Warning messages: 1: In min(x) : nenhum argumento não faltante para min; retornando Inf 2: In max(x) : nenhum argumento não faltante para max; retornando -Inf 3: In min(x) : nenhum argumento não faltante para min; retornando Inf 4: In max(x) : nenhum argumento não faltante para max; retornando -Inf 5: In plot.window(...) : "label.groups" não é um parâmetro gráfico []'s Edson Lira Estatístico Manaus-Amazonas

1 0

Erro no multcompBoxplot
by Edson Lira 08 Mai '11

08 Mai '11

Boa noite, a rotina abaixo: require(multcompView) difat<-c(3.078444,3.567257,4.023913,3.7525714) names(difat) <- c("A-B", "A-C", "B-C", "A-D") (mcT4 <- multcompTs(difat)) mcT4.1 <- plot(mcT4, label.groups=0.05) (mcL4 <- multcompLetters(difat, Letters=LETTERS)) mcL4.1 <- plot(mcL4, label.groups=0.05) Está mostrando a mensagem abaixo: Erro em plot.window(...) : valores finitos são necessários para 'xlim' Além disso: Warning messages: 1: In min(x) : nenhum argumento não faltante para min; retornando Inf 2: In max(x) : nenhum argumento não faltante para max; retornando -Inf 3: In min(x) : nenhum argumento não faltante para min; retornando Inf 4: In max(x) : nenhum argumento não faltante para max; retornando -Inf 5: In plot.window(...) : "label.groups" não é um parâmetro gráfico O erro parece ser nos limites, não consegui identificar aonde inseri-los. Alguém teria alguma idéia? Edson Lira Estatístico Manaus-Amazonas

1 0

Re: [R-br] Microdados PNAD
by Régis B. de Oliveira 07 Mai '11

07 Mai '11

Edson, obrigado pela atenção. Os arquivos que utilizo estão disponíveis em (dados.zip): ftp://ftp.ibge.gov.br/Trabalho_e_Rendimento/Pesquisa_Nacional_por_Amostra_de_Domicilios_anual/microdados/2009/ Eu não preciso especificar para o R onde começa e onde termina cada variável do arquivo? Aproveito envio outra questão/dúvida. Você sugere o uso de softwares tipo o Tinn-R? Obrigado, Régis.

3 3

Erro no multcompBoxplot
by Edson Lira 07 Mai '11

07 Mai '11

Boa noite a todos, estou rodando um gráfico do pacote multcompView, tenho a variável altura (h), e a variável at que é a combinação dos níveis de dois fatores (ap,am, nap,nam). Quando rodo, aparece na tela gráfica, mas some rapidamente e me mostra a mensagem de erro. Erro em plot.window(...) : valres finitos são necesários para 'xlim' Alem disso: Warning messages: 1: In min(x) : nenhum argumento não faltante para min; retrnando Inf 2: In min(x) : nenhum argumento não faltante para max; retrnando -Inf Alguém poderia e dar uma dica? Edson Lira Estatístico Manaus-Amazonas

3 2

Modelo de efeito misto com correlacao AR1
by Bernardo Rangel Tura 06 Mai '11

06 Mai '11

Caros colegas, Estou precisando de um ajuda para um modelo de efeito misto com correlação AR1 A auto correlação se justifica pois os pacientes devem apreender e obter melhores resultado de um bloco para outro então usei os seguinte comando ------------------------------------------------------- base <- read.csv('jogov2.csv') require(nlme) modelo <- lme(valor~Gr+Idd+Escol,random=~1|Bloco,correlation=corAR1(, form=~1|Bloco) ,data=base) summary(modelo) ----------------------------------------------------- Que resultou em ----------------------------------------------------- Linear mixed-effects model fit by REML Data: base AIC BIC logLik 3928.017 3958.571 -1957.009 Random effects: Formula: ~1 | Bloco (Intercept) Residual StdDev: 2.480127 6.846336 Correlation Structure: AR(1) Formula: ~1 | Bloco Parameter estimate(s): Phi 0.1407996 Fixed effects: valor ~ Gr + Idd + Escol Value Std.Error DF t-value p-value (Intercept) -1.7057456 1.8249072 577 -0.934703 0.3503 GrB 1.7242603 0.6606029 577 2.610131 0.0093 Idd -0.0479666 0.0178547 577 -2.686502 0.0074 Escol 0.3441736 0.0926290 577 3.715615 0.0002 Correlation: (Intr) GrB Idd GrB -0.068 Idd -0.415 -0.101 Escol -0.602 -0.227 0.001 Standardized Within-Group Residuals: Min Q1 Med Q3 Max -3.16689925 -0.56838036 -0.09641692 0.55523377 3.66357685 Number of Observations: 585 Number of Groups: 5 ------------------------------------------------------- Agora vem o problema eu preciso saber se o efeito aleatório tem um média diferente diferente entre o grupo B e o grupo A e além disso se o valor de Phi é diferente para os dois grupos Alguém pode me ajudar ? Ps. Base em anexo -- []s Tura

4 4