Junho 2014 - R-br - DINF UFPR

Plotar Mapa
by Edson Lira 11 Jun '14

11 Jun '14

Estou tentando usar a rotina abaixo para plotar o mapa library(maptools) setwd(choose.dir());getwd() ### defina o diretório de trabalho mapmanac <- readShapePoly("13MIE250GC_SIR.shp") mapmanac@data Está me dando o erro: Shapefile type: PolygonM (25), # of Shapes: 13 Erro em read.shape(filen = fn, verbose = verbose, repair = repair) : Shapefile type not (yet) handled by this function Abaixo o link com o arquivo zipado. http://www.datafilehost.com/d/d5c3e0f2 [ ]'s. Edson Lira Estatístico Manaus-Amazonas

3 2

Teste de Hipótese de Coeficientes de uma regressão
by Adriano Borges Costa 11 Jun '14

11 Jun '14

Boa Tarde, Sou iniciante no uso do R e estou com uma dívida bastante simples. Estou fazendo a análise de uma política pública e quero medir a relação entre o grau de execução da política e características das gestões municipais. Para tanto, estou usando como variáveis de controle informações socioeconômicas das cidades (população, renda, IDH, etc). Para tanto, estou rodando uma regressão linerar simples mais ou menos desta forma: Modelo<-lm(ExecuçãodoPrograma~VariáveisdeGestãoPública+Controles) O Modelo deu certo e estou conseguindo fazer a análise dos parâmetros. No entanto, preciso fazer um teste de hipótese para verificar se alguns parâmetros são iguais a zero e se uns são iguais aos outros. Como faço este teste no R? Eu devo usar as funções test.t() e var.test()? Obrigado -- Adriano Borges Costa

7 10

Sugestões sobre a melhor abordagens estatística
by Alisson Lucrecio 11 Jun '14

11 Jun '14

Caro Colegas, Boa tarde. Estou fazendo uma análise para um colega, mas estou com duvida sobre a melhor abordagem usar. Ele me enviou uma planilha com dados de matéria seca, numero nódulos (rizóbio), matéria seca dos nódulos, numero de vargens por planta, rendimendo de grão de feijão em 10 regioẽs diferente submetido a 8 tratamento diferente como com e sem nitrogênio e inoculado com 6 estirpes de rizóbio diferentes. O objetivo dele é comparar a eficiência das estirpes de rizóbio em relação a com e sem nitrogênio e selecionar a melhor estiper de rizóbio nos diferente locais. Algum colega poderia me ajudar com sugestões sobre a melhor abordagem estatística? Obrigado -- Alisson Lucrecio da Costa

2 5

Plotar categorias em clusplot
by Rodrigo Muniz 11 Jun '14

11 Jun '14

Caros membros da lista, gostaria de solicitar a ajuda de vocês mais uma vez, para plotar categorias em gráficos de cluster. Segue um exemplo. No exemplo abaixo temos um data frame com com as variáveis x,y e z, e tais variáveis estão dentro de categorias elencadas de A-F. A partir da função "kmeans" as variáveis foram agrupadas em 5 grupos, e plotados utilizando a função "clusplot". Eu gostaria de plotar no mesmo gráfico algum identificador das categorias, como cores, ou mesmo símbolos. Desculpem pelo exemplo simplista. Desde já agradeço. ex=data.frame(cat=rep(c("A", "B","C","D","E","F"), c(10, 20, 5, 15, 20,30)) , x=rnorm(100), y=rnorm(100), z=rnorm(100)) clus=scale(ex[,2:3]) fit=kmeans(clus, 5) clusplot(clus, fit$cluster, color=T, shade=T, label=2, lines=0) Att... -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) E-mail - muniz.ra(a)usp.br Cel (19) 98300-4333 (Pessoal)

1 0

Evolução de elementos em grupos
by Roney Fraga Souza 11 Jun '14

11 Jun '14

Caros, Estou tentando resolver um problema que nem sei por onde começar. Tenho dois data.frame() > data1 name grupo 1 Banker, 1986, V32, P1613 1 2 Banker, 1986, V32, P30 1 3 Macmillan, 1986, V60, P41 1 4 Charnes, 1985, V30, P91 1 5 Sherman, 1985, V9, P297 3 6 Nunamaker, 1985, V6, P50 2 7 Sherman, 1984, V4, P35 3 8 Banker, 1984, V17, P35 1 9 Banker, 1984, V30, P1078 1 10 Nunamaker, 1983, V18, P183 2 11 Parks, 1983, V12, P337 2 12 Charnes, 1981, V27, P668 3 13 Bessent, 1980, V16, P57 2 14 Charnes, 1978, V2, P429 1 > data2 name grupo 1 Sengupta, 1987, V18, P2279 2 2 Sengupta, 1987, V8, P93 2 3 Thanassoulis, 1987, V38, P397 1 4 Sengupta, 1987, V14, P117 2 5 Banker, 1986, V32, P1613 1 6 Banker, 1986, V32, P30 1 7 Macmillan, 1986, V60, P41 1 8 Charnes, 1985, V30, P91 1 9 Sherman, 1985, V9, P297 3 10 Nunamaker, 1985, V6, P50 2 11 Sherman, 1984, V4, P35 3 12 Banker, 1984, V17, P35 1 13 Banker, 1984, V30, P1078 1 14 Nunamaker, 1983, V18, P183 2 15 Parks, 1983, V12, P337 3 16 Charnes, 1981, V27, P668 3 17 Bessent, 1980, V16, P57 2 18 Charnes, 1978, V2, P429 1 e desejo comparar quantos elementos de cada grupo de data2 fazem parte de cada grupo de data1. Tentando ser mais claro: #quantos elementos de subset(data2, grupo==1) # fazem parte de subset(data1, grupo==1) subset(data1, grupo==2) subset(data1, grupo==3) #depois o raciocínio se repete para o grupo 2 de data2 #quantos elementos de subset(data2, grupo==2) # fazem parte de subset(data1, grupo==1) subset(data1, grupo==2) subset(data1, grupo==3) Por fim, desejo acompanhar quanto o grupo 1 do objeto data1 contribuiu para constituir o grupo 1 do objeto data2, e assim por diante. Para baixar esses dados, segue link: http://www.datafilehost.com/d/4428e6f8 Atenciosamente Roney

2 10

Encoding e função postForm
by Paulo Nogueira Starzynski 09 Jun '14

09 Jun '14

Senhores, acho que a questão é simples mas eu me confundo com encodings no R. No exemplo abaixo estou utilizando a função postForm para informar um número de rastreador no site dos correios. Em seguida, utilizo o pacote XML para limpar o código fonte e ficar com o que me interessa, que são os dados da carta em trânsito, desde Foz do Iguacu até minha casa. Está funcionando, mas eu gostaria que os acentos estivessem corretos. O resultado que obtenho no R é esse: [1] "\r04/06/2014 \r17:06 FOZ DO IGUACUÂ /Â PR\r\r\r\r\rObjeto postado \r \r" Marquei o enconding como "ISO-8859-1" porque é o que aparece no começo do código fonte da página: <meta http-equiv="Content-Type" content="text/html; charset=ISO-8859-1"> Como posso melhorar? Segue o código: require(RCurl) require(XML) add<-"http://www2.correios.com.br/sistemas/rastreamento/resultado.cfm" src<-postForm(add, .params = c(objetos="JG310431462BR", btnPesq="Buscar", .enconding="ISO-8859-1")) tree<-htmlTreeParse(src, useInternalNodes=T) xpathSApply(tree, '//table', xmlValue) Grato, Paulo Nogueira Starzynski

2 5

Dúvidas read.geodata
by Alex 09 Jun '14

09 Jun '14

Bom dia aos colegas da lista! Após usar a função "read.geodata", percebi que no resultado da leitura não são exibidas todas as casas decimais presentes no arquivo de origem. Na verdade o mesmo é arrendondado. Pergunta: Este seria um erro da função ou apenas uma forma de exibição truncada/arredondada? O trecho do código, a exibição parcial dos resultados e os dados presentes no arquivo de origem estão abaixo. *# Entrada de dados* (dadosP<- read.geodata("arquivo.csv", head=T, coords=4:5, data.col=9, sep="\t", dec=",")) *Resultado:* $coords E N [1,] 693067.3 6997731 [2,] 693568.8 7000175 [3,] 695651.3 6999927 [4,] 701238.3 7001212 [5,] 703768.7 7000372 [6,] 712850.4 6997434 [7,] 712848.8 7001789 [8,] 715729.2 6999106 [9,] 722978.1 6999740 [10,] 725936.6 6999752 [11,] 729037.8 7001935 [12,] 737063.3 7001920 ................ *Nota:* observem que a coordenada E apresenta 1 (uma) casa decimal, enquanto a coordenada N nenhuma. As coordenadas do arquivo de origem estão abaixo. *Dados aquivo de origem.* E N 693067,3 6997730,91 693568,8 7000174,95 695651,27 6999927,37 701238,26 7001211,82 703768,71 7000371,5 712850,37 6997434,3 712848,78 7001788,7 715729,22 6999105,77 722978,07 6999739,94 725936,55 6999752,02 729037,81 7001934,93 737063,33 7001920,08 Desde já agradeço a todos, Abraços, Alex Santos

2 2

[OFF-Topic] Lista de discussão R de Portugal - R-pt.
by walmes . 06 Jun '14

06 Jun '14

Senhores, Venho divulgar que os usuários de Portugal, talvez seguindo nosso exemplo, criaram uma lista de discussão para eles. A lista está no google groups, é fechada para não membros e teve a primeira mensagem hoje (começou a algumas horas). A administradora, Marta Rufino, me procurou e eu sugeri que ela redirecionasse mensagens para o Nabble tal como a nossa lista faz. Com isso a lista se torna pública para consulta ao histórico. Para assinar a lista portuguesa entre em https://groups.google.com/forum/#!forum/r-pt Para assinar a lista espanhola https://groups.google.com/forum/#!forum/r-es-socios Se ela acatar a sugestão de redirecionar para o Nabble e enviar o link eu divulgo junto na página R-br no LEG e na home do Nabble. Participar da lista deles (e eles da nossa) é uma forma muito positiva de estreitar relações científicas/acadêmicas além de promover/difundir o R. Participem. Grato. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

4 4

reamostragem automatizada
by Samuel Carvalho 06 Jun '14

06 Jun '14

Boa tarde caros colegas da r-br, Gostaria da opinião dos senhores em um processo de automatização de geração de amostras (reamostragem) que iniciei. A ideia é gerar n amostras de tamanhos diferentes e verificar se as médias são estatísticamente iguais. Para fins de exemplificação, é como se eu fosse no campo e medisse inicialmente altura de 5 árvores, voltasse nessa área e medisse 3, posteriormente 2 e assim por diante para testar minha suficiência amostral. ## Segue um breve RCMR #### Reamostragem #### h <- c(18.79, 12.34, 14.11, 15.24, 13.82, 14.54, 13.77, 13.99, 10.00, 16.80, 10.49, 20.00, 14.19, 14.09, 14.65, 12.22, 16.21, 13.90, 16.29, 15.39); parc <- rep(1:4, each=5); arv <- rep(1:5, 4); dados <- data.frame(parc, arv, h); dados[1:6,] rm(arv, parc, h); dados_sorteio <- dados[unlist( with(dados, tapply(seq(nrow(dados)), parc, sample, size=3) ) ),] dados_sorteio # neste passo eu gero um novo dataframe reduzindo o tamanho da minha amostra, porém fico restrito a uma única comparação e a ideia é gerar esse procedimento n vezes e comparar as médias e também a estabilização do coeficiente de variação (curva de CV x n). Agradeço pela gentileza e contribuição de todos. Abraços Samuel -- Prof. Samuel P C Carvalho Manejo Florestal e Métodos Quantitativos Faculdade de Engenharia Florestal - UFMT (065)3615-8679 / Ramal 218 www.ufmt.br/fenf

2 1

parâmetro com variavel não definida
by Waldelene Moura 05 Jun '14

05 Jun '14

Olá pessoal, se alguem puder me ajudar nesse problema agradeço. É o seguinte Preciso passar como parâmetro de uma função uma matriz que possui nos seus elementos uma variável p, sendo que essa variável p só vai assumir valores dentro da função chamada, exemplo, A = cbind(c(p+5, 6), c(10, p+8)) calcular_matrizes (A) dentro da função calcular_matrizes tenho que achar valores para essa matriz para varios valores de p A pergunta é, como passar esse parâmetro, no caso a matriz A na função, sem dar erro de que a variável p não foi definida? (object 'p' not found). Waldelene Moura Analista de Sistemas Universidade Federal do Maranhão - PROEN

2 4