- R-br - DINF UFPR

[OFF] Scientific Linux
by Pedro Rafael 02 Set '12

02 Set '12

Pessoal alguém já ouviu falar no Scientific Linux? (http://en.wikipedia.org/wiki/Scientific_Linux) Vocês sabem informar se essa distrbuição é uma boa distribuição Linux? [ ], Pedro Rafael Diniz Marinho.

2 2

Plotar pontos no mapa de Monte Santo - Bahia com as bibliotecas
by Mauricio Cardeal 02 Set '12

02 Set '12

Desculpem, não havia escrito o carregamento das bibliotecas Olá caros colegas ! Estou tentando plotar alguns pontos no mapa de Monte Santo, Bahia através do Rgoogle e rgdal, mas o mapa que aparece ao final não é o de Monte Santo. Por favor, o que estou fazendo errrado. Segue o CMR: #bibliotecas library(RgoogleMaps) library(rgdal) #dados de latitude bdlt <- c(-43.565,-43.564,-43.563,-43. 565,-43.564) #dados de longitude bdlg <- c(-10.518,-10.466,-10.412,-10.516,-10.474) dados <- as.data.frame(cbind(bdlt,bdlg)) names(dados) <- c("lat","lon") attach(dados,warn.conflicts=F) dados dados$size <- "tiny" dados$col <- "red" dados$char <- "" #get the bounding box: bb <- qbbox(lat = dados[,"lat"], lon = dados[,"lon"]) #download the map: mapa <- GetMap.bbox(bb$lonR, bb$latR, destfile = "dados.png", GRAYSCALE =F, markers = dados); #determine the max zoom, so that all points fit on the plot zoom <- min(MaxZoom(latrange=bb$latR,lonrange=bb$lonR)); plot.new() PlotOnStaticMap(mapa,lat = dados[,"lat"], lon = dados[,"lon"], cex=1.5,pch=20,col=c('blue', 'green', 'red'), add=F); Agradeço a ajuda

1 0

Plotar pontos no mapa de Monte Santo - Bahia
by Mauricio Cardeal 02 Set '12

02 Set '12

Olá caros colegas ! Estou tentando plotar alguns pontos no mapa de Monte Santo, Bahia através do Rgoogle e rgdal, mas o mapa que aparece ao final não é o de Monte Santo. Por favor, o que estou fazendo errrado. Segue o CMR: #dados de latitude bdlt <- c(-43.565,-43.564,-43.563,-43.565,-43.564) #dados de longitude bdlg <- c(-10.518,-10.466,-10.412,-10.516,-10.474) dados <- as.data.frame(cbind(bdlt,bdlg)) names(dados) <- c("lat","lon") attach(dados,warn.conflicts=F) dados dados$size <- "tiny" dados$col <- "red" dados$char <- "" #get the bounding box: bb <- qbbox(lat = dados[,"lat"], lon = dados[,"lon"]) #download the map: mapa <- GetMap.bbox(bb$lonR, bb$latR, destfile = "dados.png", GRAYSCALE =F, markers = dados); #determine the max zoom, so that all points fit on the plot zoom <- min(MaxZoom(latrange=bb$latR,lonrange=bb$lonR)); plot.new() PlotOnStaticMap(mapa,lat = dados[,"lat"], lon = dados[,"lon"], cex=1.5,pch=20,col=c('blue', 'green', 'red'), add=F); Agradeço a ajuda Mauricio Cardeal UFBA

1 0

Tutorial do pacote Hmisc
by Carlos Andrade 01 Set '12

01 Set '12

Olá, pessoal da lista. Alguém referência de links de tutorias para o pacote Hmisc? -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

2 1

Duvida de estatistica - GLM binomial
by Augusto Ribas 01 Set '12

01 Set '12

Ola pessoal. Eu estava olhando aqui uns dados com um amigo e me surgiu uma duvida. O problema é o seguinte: Os dados são referentes a ocorrência de uma doença em passarinhos. Pegamos muitos passarinhos em redes em um lugar e testamos se ele esta doente ou não. Durante a captura, registramos varias caracteristicas dos passarinhos, entre elas o sexo e a idade. Ai gostaríamos de perguntar se essas variáveis influenciam na chance do passarinho estar doente ou não. O problema começa no seguinte, a idade é fácil de determinar, são 3 classes, jovem, sub-adulto e adulto e isso é feito sem muito erro, Agora o sexo, era pra ser macho ou femea, mas quando o passarinho é muito novo, é impossível saber o sexo dele, então a gente tava colocando indeterminado, pq visualmente não da pra saber o sexo do passarinho, e não tem outro jeito. Então eu tenho 2 variáveis, de 3 níveis cada uma e fiz um glm binomial. So que qd eu faço a analise o sexo da significativo (a classe indeterminado), sendo que os indeterminados tem menor chance de estar doente. Mas na verdade eu acho que o que interessa é a idade (Que faz sentido biológico ser significativo). Se vc olha o gráfico, vc ve o padrão é o mesmo no sexo e na idade. E a inabilidade de identificar o sexo nos passarinhos vem exatamente dele ser jovem demais, não desenvolveu ainda. Ai uma coisa confunde com a outra. Será que alguem tem aluma sugestão de o que fazer? Tem como representar o sexo como "Não sei", é um dos 2 níveis mas pra esse passarinho eu não sei? Eu não sei como fazer isso. É possível fazer isso pro sexo no glm no meu caso? Segue um exemplo para ilustrar mais ou menos o que eu to vendo: #Gerando dados de exemplo set.seed(666) n<-90 #eu acredito que a idade tem um efeito significativo #ja que por exemplo, o cara muito jovem vai ter tempo de exposição menor #a indivíduos doentes idade<-sample(c("Adulto","Sub-Adulto","Jovem"),n,replace=T) idade efeito<-ifelse(idade=="Jovem",1,0) efeito diagnostico<-rbinom(n,1,0.75+efeito*(-0.25)) #E o sexo do passarinho não tem efeito, mas os jovens sempre são inderterminados sexo<-NA for(i in 1:length(idade)) { if(idade[i]=="Jovem") { sexo[i]<-c("Indeterminado") } else{ sexo[i]<-sample(c("Macho","Femea"),1) } } dados<-data.frame(sexo,idade,diagnostico) #so completando os casos dados<-rbind(dados,c("Femea","Jovem",0)) dados<-rbind(dados,c("Indeterminado","Adulto",0)) dados<-rbind(dados,c("Macho","Jovem",0)) dados<-rbind(dados,c("Indeterminado","Sub-Adulto",0)) #amostras table(dados$sexo,dados$idade) head(dados) #Grafico #no final o que eu vejo é mais ou menos isso: par(mfrow=c(1,2),cex=0.6) barplot(table(diagnostico,sexo),ylim=c(0,40),legend.text =c("Parasitado","Não Parasitado")) barplot(table(diagnostico,idade),ylim=c(0,40),legend.text =c("Parasitado","Não Parasitado")) #O indeterminado e Jovem é mais meio a meio no numero de parasitados que as outras classes #mas eu vejo um resultado significativo pro sexo, ou jovem se mudar a ordem #isso esta errado modelo01<-glm(diagnostico~factor(idade)+factor(sexo),family="binomial") summary(modelo01) plogis(coef(modelo01)[1:2]) #eu acredito que cai naquele problema do aov, de experimento não balanceado e tal #se olhar individualmente tudo ok, mas não consigo olhar as 2 variaveis juntas. modelo02<-glm(diagnostico~factor(idade),family="binomial") summary(modelo02) plogis(coef(modelo02)[1:2]) modelo03<-glm(diagnostico~factor(sexo),family="binomial") summary(modelo03) plogis(coef(modelo03)[1:2]) Se alguem puder dar uma olhada e dar alguma sugestão, fiz um exemplo, mas os dados deixam um grafico exatamente como aquele, se alguem quiser olhar os dados originais eu posso mandar. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

6 16

Mais códigos HTML e função readlines
by Paulo Nogueira 01 Set '12

01 Set '12

Pessoal, não sei se a lista se presta a esclarecer dúvidas como essa minha, mas pelo menos estou usando R para resolver os problemas. rs Ultimamente tenho trabalhado bastante com determinadas buscas na web para coletar dados e estruturá-los. Porém, como meus conceitos de web se restringem ao conhecimento de usuário, e não de desenvolvedor, acabo tendo alguns problemas. Vou exemplificar uma das minhas dificuldades com um problema prático: Utilizar readlines() e os pacores RCurl, XML, etc, para tratar o código fonte é algo que já faço até que bem. Mas problemas surgem quando o conteúdo que eu preciso extrair de informação não vem no código fonte? Vejam esse exemplo: http://www.skyscanner.com.br/passagens-aereas/saoa/slz/121002/130409/tarifa… Após carregar toda a busca, o site SkyScanner exibe os preços encontrados para o trecho São Paulo - São Luis. [1] No navegador Firefox ou Chrome, ao clicar em "Exibir Código Fonte", não aparecem as informações do vôos, apenas o código que estrutura a página. [2] Ao clicar em "Inspecionar Elemento" eu consigo analisar o código fonte do meu interesse. [3] Ao utilizar readlines passando o a url acima o R importa o código fonte [1] e não o [2], com as informações de interesse. Imagino que isso tenha a ver com as funções de javascript e coisas desse tipo... Alguém sabe contornar essa situação com o R? Se não, o que devo estudar para agilizar a busca dessa solução? Estou meio perdido nessa parte e agradeço ajudas de qualquer tipo. rs Peço desculpas se por acaso isso for muito off-topic. Abraços, Paulo Nogueira

4 5

Ajuda com o Gedit
by Luciano F. Sgarbi 01 Set '12

01 Set '12

Olá pessoal, Alguém ai usa o gedit no R, com o sistema operacional Windows 7? Estou tentando faze-lo funcionar mas não tem jeito... Se alguém der um help fico muito agradecido. Abraços Luciano -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732

2 2

Bug do Rcmdr?
by Emerson Cotta Bodevan 31 Ago '12

31 Ago '12

Prezados, recorro aos senhores para mais uma dúvida. No momento estou utilizando o R para windows... Instalei o R e depois o Rcmdr. Após a instalação do R criei um atalho na barra de ferramentas. Estou com o seguinte problema: - ao iniciar o R deste atalho e carregar o Rcmdr pelo comando (require(Rcmdr)) aparece no menu *Dados *a opção de *importar arquivos *e lá a opção de *arquivos excel, access, etc*. - ao iniciar o R por um atalho em um subdiretório de um projeto esta opção *arquivos excel, access, etc*. não aparece. Alguém sabe como resolver? Agradeço qualquer ajuda. Abs, Emerson. -- "A educação não muda o mundo. A educação muda as pessoas. As pessoas mudam o mundo." * Paulo Freire.*

2 2

[Off-topic] Corretor ortográfico pt_br Kile 2.1
by Andre Oliveira 31 Ago '12

31 Ago '12

Pessoal como posso colocar corretor ortográfico para o Kile 2.1 no Ubuntu 12.04? obrigado André Oliveira Souza

1 0

Manual sobre o pacote tcltk
by Ary Casado 31 Ago '12

31 Ago '12

Pessoal Estou desenvolvendo um pequeno sistema com interfaces gráficas para o meu trabalho de dissertação. Escolhi o pacote tcltk para esse trabalho porque achei muito fácil de lidar com ele, porém quanto mais escrevo e crio forms mais aparecem funções, eventos, propriedades e métodos. Alguém conhece algum manual bem didático sobre o tcltk???? O site www.sciviews.org/rgui já me ajudou bastante com vários exemplos, porém não explica tudo. Arykerne Nascimento Mestrando Proten/UFPE

1 0