Agosto 2012 - R-br

dúvida - modelo
by Sérgio Henrique almeida da silva ju 17 Ago '12

17 Ago '12

Olá Estou modelando um conjunto de dados e a variável resposta parece ter distribuição geométrica, qual é a melhor forma de modelar esses dados? Abraços dput(dados) structure(list(AFtotal = structure(c(103L, 76L, 58L, 105L, 125L, 163L, 106L, 146L, 227L, 187L, 5L, 117L, 147L, 39L, 232L, 2L, 67L, 150L, 46L, 118L, 18L, 171L, 157L, 226L, 80L, 96L, 142L, 162L, 128L, 152L, 109L, 10L, 108L, 102L, 47L, 208L, 122L, 78L, 176L, 3L, 196L, 40L, 209L, 119L, 179L, 29L, 1L, 81L, 190L, 24L, 204L, 28L, 92L, 136L, 104L, 85L, 91L, 15L, 141L, 13L, 11L, 121L, 112L, 210L, 31L, 6L, 70L, 156L, 114L, 68L, 107L, 112L, 211L, 167L, 1L, 215L, 88L, 213L, 194L, 136L, 220L, 199L, 136L, 9L, 136L, 118L, 38L, 192L, 42L, 62L, 155L, 19L, 79L, 48L, 138L, 172L, 134L, 231L, 231L, 216L, 136L, 195L, 130L, 101L, 217L, 205L, 174L, 36L, 111L, 140L, 106L, 99L, 22L, 218L, 170L, 7L, 8L, 136L, 19L, 62L, 25L, 123L, 1L, 137L, 220L, 97L, 52L, 169L, 190L, 216L, 1L, 53L, 172L, 100L, 79L, 86L, 201L, 61L, 198L, 17L, 218L, 67L, 106L, 193L, 1L, 138L, 164L, 200L, 190L, 133L, 129L, 74L, 186L, 72L, 41L, 177L, 48L, 113L, 50L, 229L, 178L, 106L, 124L, 75L, 14L, 148L, 218L, 191L, 19L, 219L, 30L, 1L, 71L, 112L, 172L, 21L, 221L, 16L, 203L, 27L, 69L, 143L, 82L, 207L, 223L, 32L, 54L, 60L, 59L, 1L, 97L, 149L, 63L, 132L, 212L, 214L, 1L, 1L, 165L, 77L, 20L, 73L, 26L, 168L, 160L, 233L, 139L, 1L, 188L, 95L, 175L, 206L, 93L, 144L, 135L, 56L, 6L, 189L, 197L, 45L, 30L, 35L, 1L, 4L, 136L, 154L, 89L, 65L, 184L, 225L, 230L, 151L, 1L, 173L, 43L, 44L, 33L, 182L, 34L, 153L, 1L, 1L, 90L, 106L, 110L, 98L, 1L, 219L, 1L, 1L, 116L, 49L, 64L, 57L, 131L, 185L, 145L, 17L, 94L, 1L, 180L, 1L, 120L, 159L, 51L, 172L, 181L, 183L, 7L, 87L, 106L, 228L, 55L, 126L, 157L, 218L, 158L, 115L, 133L, 222L, 84L, 166L, 1L, 37L, 1L, 1L, 202L, 66L, 12L, 1L, 83L, 23L, 110L, 233L, 231L, 161L, 106L, 224L, 127L), .Label = c("0,00", "1.031,00", "1.035,00", "1.042,50", "1.064,00", "1.080,00", "1.120,00", "1.155,00", "1.171,50", "1.177,50", "1.179,00", "1.200,00", "1.204,50", "1.210,00", "1.218,00", "1.220,00", "1.230,00", "1.233,00", "1.260,00", "1.275,00", "1.287,00", "1.290,00", "1.314,00", "1.320,00", "1.356,00", "1.359,50", "1.400,00", "1.410,00", "1.418,00", "1.440,00", "1.466,00", "1.480,00", "1.500,00", "1.566,00", "1.572,00", "1.590,00", "1.611,00", "1.652,00", "1.680,00", "1.722,00", "1.734,00", "1.740,00", "1.752,00", "1.800,00", "1.818,00", "1.821,00", "1.833,00", "1.836,00", "1.854,00", "1.872,00", "1.900,00", "1.920,00", "1.938,00", "10.080,00", "10.230,00", "10.557,00", "10.680,00", "11.238,00", "12.000,00", "12.270,00", "12.420,00", "120,00", "13.644,00", "13.860,00", "14.590,00", "150,00", "165,00", "180,00", "2.002,50", "2.060,00", "2.078,00", "2.100,00", "2.120,00", "2.157,00", "2.320,00", "2.332,50", "2.376,00", "2.382,00", "2.430,00", "2.457,00", "2.475,00", "2.504,00", "2.520,00", "2.526,00", "2.532,00", "2.556,00", "2.572,00", "2.580,00", "2.616,00", "2.640,00", "2.645,00", "2.655,00", "2.706,00", "2.739,00", "2.745,00", "2.767,50", "2.784,00", "2.798,00", "2.820,00", "2.841,00", "2.895,00", "2.937,00", "2.940,00", "2.965,00", "21.515,00", "240,00", "247,50", "259,00", "26.862,00", "264,00", "273,00", "297,00", "3.019,00", "3.035,00", "3.036,00", "3.282,00", "3.306,00", "3.360,00", "3.497,00", "3.504,00", "3.510,00", "3.540,00", "3.600,00", "3.615,00", "3.660,00", "3.710,00", "3.840,00", "3.892,50", "3.900,00", "3.960,00", "3.977,00", "3.996,00", "300,00", "318,00", "330,00", "360,00", "372,00", "396,00", "4.020,00", "4.086,00", "4.092,00", "4.100,00", "4.140,00", "4.158,00", "4.200,00", "4.219,50", "4.331,00", "4.359,00", "4.554,00", "4.618,00", "4.680,00", "4.770,00", "407,50", "410,00", "420,00", "424,00", "480,00", "495,00", "5.052,00", "5.160,00", "5.340,00", "5.622,00", "5.688,00", "5.715,00", "5.745,00", "5.832,00", "5.940,00", "5.958,00", "510,00", "516,00", "537,00", "540,00", "565,50", "570,00", "575,00", "576,00", "585,00", "6.097,50", "6.285,00", "6.495,00", "6.556,00", "6.576,00", "6.720,00", "6.746,00", "6.828,00", "639,00", "645,00", "655,00", "66,00", "660,00", "678,00", "690,00", "693,00", "696,00", "7.056,00", "7.074,00", "7.200,00", "7.236,00", "7.245,00", "7.500,00", "7.600,00", "7.680,00", "7.860,00", "7.920,00", "720,00", "738,00", "754,00", "754,50", "780,00", "795,00", "8.268,00", "8.467,50", "8.592,00", "8.755,00", "8.820,00", "80,00", "801,00", "810,00", "812,00", "840,00", "876,00", "882,00", "895,50", "9.120,00", "9.150,00", "9.310,00", "9.508,00", "924,00", "930,00", "933,00", "960,00", "99,00", "990,00"), class = "factor")), .Names = "AFtotal", class = "data.frame", row.names = c(NA, -299L)) -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

4 3

Anova em faixas
by Juliano Farias 17 Ago '12

17 Ago '12

Srs, Preciso de uma ajuda. Preciso fazer uma ANOVA de um experimento com três fatores, porém o problema está no delineamento que foi em faixas. Tenho o quadro da análise da variança para esse tipo de delineamento, porém não consegui montar a programação e não achei um pacote para esse tipo de delineamento. O experimento foi realizado em blocos (B) ao acaso com quatro repetições em esquema fatorial em faixas 2 x 2 x 2. Os dois níveis do fator A, com os dois níveis do fator C foram casualizados nas parcelas principais, e os dois níveis do fator D foram casualizados em faixas perpendiculares às parcelas principais. Grato pela atenção. -- Atenciosamente, Juliano Ricardo Farias Eng. Agrônomo Doutorando em Entomologia Lab. de Resistência de Artrópodes a Pesticidas ESALQ/USP

3 4

Somar vendas por semana e por mês
by Hugo Sena 17 Ago '12

17 Ago '12

Boas! Vou fazer análise de séries temporais e preciso de colocar as vendas de alguns produtos tendo: em cada linha as vendas totais por semana e também tendo em cada linha, as vendas totais por mês. Tenho os dados em Excel, alguma ajuda para isto ser rápido e automático?. Cumps.

2 1

Regressão Gama
by Luiz Henrique da Conceicao Leal 15 Ago '12

15 Ago '12

Bom dia pessoal Construi um regressão gama e gostaria de construir um intervalo de confiança (no gráfico) para reta estimada. Alguém tem alguma sugestão? modelo=glm(y~x,family=Gamma(link = "identity")) plot(x,y,pch=20) abline(modelo,col="blue",lwd=2) Desde já agradeço Luiz Henrique Leal

1 0

Duvida no TAPPLY
by Edson Lira 15 Ago '12

15 Ago '12

Estou usando esta rotina e obtenho tapply(acj$Quantidade,list(Laboratório=acj$lab,Mês=acj$Mês),sum,na.rm=T)#Laboratórios E obtenho esta saída abaixo, que me dá o quantitativo de exames mensal, para obter a soma nas linhas e a soma nas colunas,o que devo acrescdentar a rotina? Mês Laboratório Jan Fev Mar Abr Mai Jun Jul Ago Set Out Bioquímica 24927 25529 24589 23736 29551 30175 28760 29220 24291 22189 Estudo de Hemoglobinopatias 300 226 232 201 293 360 209 243 233 208 H. Pardini 546 479 436 360 433 403 405 455 357 310 HEMATOLOGIA 8764 8877 8398 7969 9389 9360 9264 9585 8145 7548 HORMÔNIOS 1137 1276 1307 1271 1579 1522 1479 1519 1257 1286 IMUNOLOGIA ERITROCITÁRIA 1961 1942 1890 1820 2132 2206 1996 1992 1736 1564 PARASITOLOGIA/URINÀLISE 4365 4265 3477 3526 4183 4581 4581 4582 3898 3667 HEMOSTASIA 1174 987 1025 928 1266 1167 1257 1134 792 732 IMUNOLOGIA 1839 1570 1428 2220 1715 1522 1266 1490 1349 1253 SOROLOGIA 3327 3116 2863 3431 4437 4426 4294 4134 2834 1993 [ ]'s. Edson Lira Estatístico Manaus-Amazonas

3 4

Utilizando DEA com R
by Cleyton Fábio 15 Ago '12

15 Ago '12

Ola gente, tudo bem? Estou trabalhando em um projeto onde quero analisar a eficiência de municípios brasileiros utilizando Analise Envoltória de Dados. A ideia é implementar uma solução(100% baseada em software livre) que facilite este trabalho, conheci o R e descobri que existe alguns pacotes como o DEA e o Benchmarking, porem ainda estou tendo um pouco de dificuldades em utilizar a aplicação. Gostaria de saber se alguem tem alguma documentação sobre a aplicação de DEA ou Benchmarking, poderia ser artigos, monografias, tutoriais e a possibilidade de disponibilizar. O trabalho é a minha monografia! grato, Cleyton Fábio

5 4

Duvidas com spDists
by Julianna Trindade 15 Ago '12

15 Ago '12

Pessoal, tenho alguns pontos num shapefile e preciso saber as distâncias entre todos os pontos. Consegui fazer isso utilizando spDists mas a matriz que resulta não identifica os pontos pela coluna de atributo que eu preciso, apresenta apenas uma matriz quadrada utilizando 1,1; 1,2...2,1;2,2...e etc pra definir as linhas e as colunas e eu preciso que ao invés de 1,1 apareçam os atributos de uma coluna específica. Isso seria possível?? Além disso, qual unidade é utilizada nas distancias?Km?metros? O codigo eh o seguinte: rm(list=ls(all=TRUE)) library(base) library(maps) library(sp) library(maptools) library(spdep) library(classInt) library(RColorBrewer) library(shape) library(SDMTools) library(latticeExtra) pontos = readShapePoints("c:/users/Julianna/Documents/SecSaude/hospitais_produtividade.shp") proj4string(pontos) <- CRS("+init=epsg:4291") distancias = spDists(pontos) Resultando: pontosSpatialPointsDataFrame[1] distancias23x23 double matrix -- Julianna Trindade "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." Salmos 21; 13

2 2

[OFF] Biblioteca C para matrizes.
by Pedro Rafael 15 Ago '12

15 Ago '12

Pessoal gostaria de saber qual uma boa biblioteca livre de C para trabalhar com operações matriciais como inversas e produtos de matrizes. Qualquer dica será interessante. [ ] Pedro Rafael Diniz Marinho.

3 28

Vaga para Estágio - Desenvolvimento em R - Porto Alegre
by Fernando Henrique Ferraz Pereira da Rosa 15 Ago '12

15 Ago '12

PROGRAMA DE ESTÁGIO O Programa de estágio da Souza Cruz desenvolve jovens talentos, agregando valor à sua formação acadêmico-profissional. Ao final do Programa, os estagiários têm a oportunidade de participar do Programa de Trainee ou serem indicados para outros desafios na Empresa. Vaga para os cursos de Estatística / Economia / Engenharia de Produção / Processamento de Dados / Ciência da Computação e correlatos Bolsa auxílio no valor de R$1.053 (6h) e R$702 (4h) , auxílio transporte, restaurante, plano de saúde (opcional) e seguro de vida. Área de lotação: Product Centre – Suporte Estatístico. A área provê suporte estatístico para as diversas áreas da empresa (laboratórios, fábricas, desenvolvimento de produtos, entre outros). Resumo de atividades: atender às necessidades de manutenção e desenvolvimento de softwares e ferramentas estatísticas da área e dos seus clientes. A demanda por uso e aplicação de conhecimentos estatísticos será balanceada de acordo com a área de formação do candidato. A posição poderá ser assumida por alguém sem conhecimento prévio em estatística, caso no qual será dado mais enfoque ao desenvolvimento de software e programação, ou caso seja preenchida por alguém com formação em estatística, haverá maior enfoque nos métodos estatísticos. Áreas de atuação: Planos amostrais Planejamento de Experimentos Programas intra e interlaboratoriais Relatórios Gerenciais (MIS) Otimização de Processos Controle Estatístico de Processo / 6-sigma PRÉ-REQUISITOS Universitários cursando o penúltimo ou último ano de graduação. Bons conhecimentos em informática (preferencialmente deve conhecer ao menos uma linguagem de programação). DESEJÁVEL Conhecimento em R e/ou SAS, Excel VBA, lógica de programação. LOCAL DE TRABALHO Unidade Fabril de Cachoeirinha (Região Metropolitana de Porto Alegre) INSCRIÇÕES Enviar CV para fernando.henrique.rosa(a)souzacruz.com.br -- "Though this be randomness, yet there is structure in't." Fernando H Rosa - Statistician http://www.fernandohrosa.com.br / http://www.feferraz.net - Estatística, Matemática e Computação AprendaAlemao.net - Seu ponto de partida para melhorar seu Alemão! @fhrosa

1 0

Re: [R-br] microdados censo IBGE 2010
by Roney Fraga 14 Ago '12

14 Ago '12

Manoel, O Marcos F Silva criou algunas funções maravilhosas para ajudar a trablhar com microdados no R. No endereço: https://sites.google.com/site/marcosfs2006/textos baixar o arquivo Usando R para Importar Microdados.pdf e as funções no link .zip. Com base no censo de 2000, onde apenas a tabela para o estado de São Paulo tem aproximadamente 3,5 GB quando em formato .csv, acho que será um pouco complicado trabalhar com esse volume de dados diretamente no R. Para lidar melhor com isso optei por colocar os arquivos para cada estado como uma tabela do banco de dados PostgreSQL, isso resulta em um bom ganho de tempo, até porque as consultas básicas podem ser feitas no próprio PostgreSQL ou via algum pacote como RpgSQL ou RPostgreSQL. Atenciosamente Roney

3 2