Maio 2014 - R-br - DINF UFPR

Formatação com sprintf
by Rubem Kaipper Ceratti 07 Mai '14

07 Mai '14

Caros listeiros, Tenho um problema que consiste em transformar um certo número para formato 'character' com comprimento 14. Para isso tentei usar a função base::sprintf, mas para números com mais de 10 dígitos, recebo a mensagem de erro abaixo: > sprintf('%.14d', 1202342234) [1] "00001202342234" > sprintf('%.14d', 12023422344) Error in sprintf("%.14d", 12023422344) : invalid format '%.14d'; use format %f, %e, %g or %a for numeric objects Alguém sabe uma forma de contornar isso? Grato, Rubem

3 5

Emprego para estatísticos com conhecimentos de R
by Wenceslau Teixeira 06 Mai '14

06 Mai '14

With the hiring of an IT developer - Statistics, CropDesign wants to strengthen its team of statisticians that is responsible for experimental design and data analysis in the field of agricultural crops. The IT developer - Statistics will develop and improve local and global IT infrastructure, and maintain database systems to enhance the team's efficiency to perform data analysis. He/she will develop web-based systems for reporting analysis results to customers, and improve procedures through automation and standardization. In addition, he/she will be active as a statistician, and contribute to data analysis, interpretation of results and reporting to biologists. *What we expect:* You have a master degree or higher in a relevant field (IT, statistics, or related) or at least 3 years of relevant work experience. You have ample experience in computer programming (*particularly in R*) and databasing (SQL). *We are:* BASF Plant Science develops plant biotechnology solutions to increase the yield and quality of crop plants. CropDesign, the Belgian division of BASF Plant Science, focuses on rice as model crop. The company applies genomics, informatics and robotics technologies and employs over 140 employees in state-of-the-art laboratories and an automated greenhouse complex. https://www.researchgate.net/job/813404_IT_Developer-Statistics?source=prfs… -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

1 0

[off] Recomendações para programação
by FHRB Toledo 06 Mai '14

06 Mai '14

Senhores, Postado no R-Bloggers http://www.econometricsbysimulation.com/2014/05/how-to-code-something-new-i… att, FH

1 0

Selecionar e copiar conteúdo de página usando apenas o terminal linux
by walmes . 06 Mai '14

06 Mai '14

Saudações, Preciso ler os dados de uma página mas da forma tradicional, usando XML para o exame e retenção dos dados eu não consigo. Os dados são de altura, idade, peso, enfim, de jogadores de futebol do campeonato brasileiro. Estão disponíveis nessa página, por exemplo, para o Fluminense *http://www.whoscored.com/Teams/1232 <http://www.whoscored.com/Teams/1232>* O que me interessa é a tabela com o nome dos jogadores, dados biométricos e esportivos. Um selecionar e copiar dessa porção que me interessa traz esse conteúdo. R Name Pos Age cm kg Apps Goal A Yel Red SpG PS% AW MoM Rt 1 Rafael Sóbis AM(LR),FW 28 172 67 2 2 - - - 3.5 78.8 1.5 - 7.82 2 Carlinhos D(L) 27 176 72 2 - - - - 1 81.1 1 1 7.57 3 Gum D(C) 28 189 86 2 - - - - 0.5 75.2 2.5 1 7.51 4 Wagner AM(C) 29 172 71 2 - - 1 - 3 90 - - 7.43 5 Bruno Vieira D(R) 28 178 69 2 - - - - - 85.9 1.5 - 7.23 6 Darío Conca M 30 167 58 2 - - 1 - 3 85.9 - - 7.21 7 Fred FW 30 185 75 2 1 2 - - 3 75.6 1.5 - 7.19 8 Diguinho DM(C) 31 171 70 2 - - - - 0.5 84.7 1.5 - 7.14 9 Elivelton D 21 179 81 2 - - - - - 74.6 0.5 - 6.86 10 Diego Cavalieri GK 31 191 86 2 - - - - - 46.2 - - 6.69 11 Jean DM(C) 27 170 70 2 - - - - 1 88.6 - - 6.61 12 Biro Biro AM(L) 19 164 55 0(1) - - - - 1 60 - - 6.47 13 Edwin Valencia M 29 181 85 0(2) - - - - - 81.8 1 - 6.31 14 Rafinha M(R) 20 177 74 0(1) - - - - - 50 - - 6 15 Walter FW 24 178 88 0(1) - - - - - 60 - - 5.95 16 Felipe GK 26 193 90 - - - - - - - - - - 17 Chiquinho AM(C) 24 170 69 - - - - - - - - - - 18 Kenedy FW 18 182 77 - - - - - - - - - - 19 Wellington Carvalho D 21 185 84 - - - - - - - - - - 20 Wellington Silva D 26 177 64 - - - - - - - - - - Quando eu peço o código fonte da página (acesso pelo botão direito do mouse no navegador firefox), não aparecem esses valores lá. Se aparecessem dava para usar o XML. Pelo jeito os dados são trazidos por uma consulta à uma base de dados on the fly para colocar na página. A única forma que eu sei de pegar é copiando para área de transferência. Eu tenho feito assim: 1. Abro a página de cada time no navegador e dou ctrol+a para selecionar todo o conteúdo; 2. Dou ctrol+c para copiar o conteúdo para área de transferência; 3. Dou control+v dentro de um arquivo texto; 4. Leio o arquivo texto no R, elimino a sujeira, organizo em colunas até chegar um data.frame. O que eu quero é substituir passos gráficos de 1 à 3 por algum comando do teminal linux que pudesse pegar a URL que eu passar e salvar o conteúdo em um aquivo que eu der o nome. Com isso elimino intervenção humana e rapidamente baixo todos os times, de todos os torneios para ter excelentes dados para trabalhar em aula. Se alguém se interessar, os dados da liga européia (56 times) que peguei no braço estão em *http://www.leg.ufpr.br/~walmes/data/euro_football_players.txt <http://www.leg.ufpr.br/~walmes/data/euro_football_players.txt>* Dei umas googladas e vi alguns comandos que, infelizmente, ainda não trazem o que eu preciso, como *sudo apt-get install lynx* *lynx -dump http://www.whoscored.com/Teams/1232 <http://www.whoscored.com/Teams/1232>* O recorte do resultado do comando está abaixo. View: [41]Overall [42]Home [43]Away R Name Pos Age cm kg Apps Goal A Yel Red SpG PS% AW MoM Rt <------ a tabela que eu queria era pra estar *Players shaded are players who are not currently active in team. <------ entre essas duas linhas, mas não. (Loaned, sold, etc..) Tournaments: [44]Brasileirão View: [45]Overall [46]Home [47]Away Caso alguém tenha algum caminho a apontar, fico 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 ==========================================================================

6 8

Problemas no ajuste de semivariogramas teóricos
by Helder Gramacho 05 Mai '14

05 Mai '14

Boa noite pessoal, Gostaria de uma ajuda, estou tentando fazer o ajuste de semivariogramas teóricos pelo método dos mínimos quadrados, entretanto, aparentemente não está sendo possível realizar o ajuste, ficando como na figura abaixo. Será que se trata de algum problema com o código que estou executando, ou os dados é que não permitem mesmo o ajuste. Reproduzi o código abaixo da figura e o arquivo está no Dropbox. Desde já agradeço qualquer ajuda, [image: Imagem inline 1] # Download de um arquivo do Dropbox links <- c("https://www.dropbox.com/s/8blubgumcsss834/cobre.csv") tokens <- gsub("^.*/s/","",dirname(links)) fileNames <- basename(links) newLinks <- file.path("http://dl.dropbox.com/s", tokens, fileNames); newLinks for (a in newLinks) { tryCatch(download.file(a, dest=basename(a), mode='wb'), error=function(...) print("Falha no download!"))} cobre<-read.table(file="cobre.csv", sep=",", header=T, dec=".") attach(cobre) cobre<-as.geodata(cbind(cobre$V1,cobre$V2,cobre$V3)) # Semivariograma Experimental par(mfrow=c(1,1),xpd=F) var.cu<-variog(cobre,uvec=seq(0.1,3,l=10),pairs.min=30,estimador.type="classical", direction="omnidirectional",tolerance=pi/8) plot(var.cu, xlab="distância", ylab="semivariância",main='Cobre',font.main = 3) #-------------------------------------------------------------------------------------------------------- # Ajustando Semivariograma Teórico pelo Método dos Mínimos Quadrados #-------------------------------------------------------------------------------------------------------- # Modelo exponencial exp.ols.cu<-variofit(var.cu,ini=c(235, 0.41),weights="equal", cov.model="exp") exp.ols.cu plot(var.cu, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Exponencial') lines(exp.ols.cu, col="blue") summary(exp.ols.cu) ## Modelo esférico sph.ols.cu<-variofit(var.cu,ini=c(235, 0.41),weights="equal", cov.model="sph") sph.ols.cu plot(var.cu, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Esférico') lines(sph.ols.cu, col="blue") summary(sph.ols.cu) ## Modelo gaussiano gaus.ols.cu<-variofit(var.cu,ini=c(235, 0.41),weights="equal", cov.model="gaus") gaus.ols.cu plot(var.cu, xlab='Alcance', ylab='Semivariância', main='Semivariograma OLS - Gaussiano') lines(gaus.ols.cu, col="blue") summary(gaus.ols.cu) *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

4 8

Mover coordenadas de um raster (Georreferenciamento)
by ASANTOS 05 Mai '14

05 Mai '14

Boa noite pessoal, Tenho um novo problema, tenho um raster (dem.crop.final.tif) que representa a elevação do terreno (parte de um objeto Topodata do INPE) e preciso georreferenciar o raster, sendo que tenho as coordenadas atuais (x.orig e y.orig) que estão erradas e preciso corrigir com as coordenadas certas (x.new e y.new) para todos os pixels do raster baseado nestes 6 pontos, alguém sabe se é possível realizar essa transformação sem utilizar o GRASS ou QGIS, segue CRM dos dados: ### <code r> lapply(c('rgdal', 'maptools', 'shapefiles'), require, character.only=TRUE) #Download do controno da área em shapefile links <- c( "https://www.dropbox.com/s/3ph2630xicpexo0/dem.crop.final.tif") if (sum(sapply(basename(links), file.exists))<1) { lapply(links, function(a) tryCatch(download.file(a, dest=basename(a), mode='wb'), error=function(...) message("ERRO!")))} else message('OK!') dem <- readGDAL('dem.crop.final.tif') x.orig<-c(-42.33960, -42.33937, 42.34546, -42.34738, -42.21652, -42.40492) y.orig<-c(-19.23349, -19.23514,-19.23217, -19.22818,-19.05137,-19.50393) x.new<-c(-42.33929, -42.33921, -42.34515, -42.34714, -42.21620, -42.40471) y.new<-c(-19.24486, -19.24655,-19.24350,-19.23955,-19.06276,-19.51528) ### </code> Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

2 1

Análise de Cluster
by Roberto Guevara 05 Mai '14

05 Mai '14

Prezados colegas usuários de R, boa tarde. Estou tentando realizar uma análise de cluster e minha dúvida é se há como utilizar valores de similaridade, ao invés de dissimilaridade, entre as unidades a serem agrupadas. Sei que posso obter os valores de similaridade subtraindo de 1 os valores de dissimilaridade. No entanto, não achei nenhum argumento na função hclust para que a análise seja feita com base nas similaridades. Quem puder me ajudar, agradeço. Peço desculpas se for um tópico repetido ou caso se trate de uma dúvida por falta de conhecimento da análise e não da linguagem R em si. Atenciosamente, -- Roberto Guevara Ferreira Lima Biólogo - CRBio 6ª Região: 73146/06-D Mestre em Zoologia Laboratório de Ecologia e Zoologia de Vertebrados Instituto de Ciências Biológicas Universidade Federal do Pará Rua Augusto Corrêa, nº 01 - Bairro: Guamá. CEP: 66075-110 Currículo Lattes: http://lattes.cnpq.br/7363382159007600

1 0

Re: [R-br] Análise de dados categorizados pela metodologia GSK
by Ivan Bezerra Allaman 05 Mai '14

05 Mai '14

Uma breve googlada http://r.789695.n4.nabble.com/Repeated-measures-with-categorical-data-td468… (s,f,p) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Análise de dados categorizados pela metodologia GSK
by Wirton Macedo Coutinho 05 Mai '14

05 Mai '14

Bom dia pessoal, Alguém conhece algum pacote do R que faz análise de dados categorizados por modelos lineares com resposta dicotômica, utilizando a metodologia GSK? Att., Wirton Macedo Coutinho Fitopatologia - Embrapa Algodão Campina Grande, PB

1 0

Fatoração LU no R
by Andre Oliveira 05 Mai '14

05 Mai '14

Pessoal alguém sabe como obter a fatoração LU de uma matriz no R? obrigado André Oliveira Souza

2 2