Maio 2012 - R-br - DINF UFPR

Criando Arquivo Shape com Linha
by Ricardo Dobrovolski 07 Mai '12

07 Mai '12

Caros R Users, Eu preciso criar um conjunto de arquivos shape, sendo que cada arquivo deve conter como objeto apenas uma linha que conecta dois pontos conhecidos. Desde já, agradeço. Atenciosamente, Ricardo Dobrovolski

1 0

windows.options() - Linux
by Luís Gustavo 07 Mai '12

07 Mai '12

Boa tarde pessoal, Estou com a seguinte dúvida: qual o comando que é similar ao windows.options() no Linux, pois este não existe para o Linux. O objetivo final é criar uma função que me gera a figura da legenda, porém esta função tem rodar tanto no Windows quanto no Linux. Pois gostaria que a figura da legenda se ajuste para os mais variados números de faixas. Abaixo segue um exemplo. (Acredito que função não a mais eficiente, qualquer sugestão é bem-vinda) Obrigado. brks <- c(0,5,10,15,20) r <- 3 cols <- brewer.pal(4,"Blues") gera_figura_legenda=function(brks, cols,r) gera_figura_legenda=function(brks, cols,r) { brks = round(brks,r) dest.fig = tempfile(pattern="legenda",fileext=".png",tmpdir=getwd()) windows.options(width=2.5, height=1.7, pointsize=12, reset = TRUE) CairoPNG(filename = dest.fig, width = 290, height = 290, bg="#FFFFFF00", pointsize = 10) par(bg="transparent", mar=c(0, 0, 0, 0)) brks = round(brks,r) plot(1,1,col="white",axes=F,col.axis="white",xlab="",ylab="") aux1 = 2:(length(brks)-1) aux2 = expand.grid(aux1,aux1) d = aux2[,1]-aux2[,2] aux3 = aux2[d==-1,] legend( "center", title = "Legenda", c(paste("<",brks[2]), paste(brks[aux3[,1]],brks[aux3[,2]],sep="-"), paste(">",brks[(length(brks)-1)])), fill=c(cols), cex=1.3, bg="#FFFFFF") dev.off() return (dest.fig) } -- Luís Gustavo Silva e Silva

5 5

Favor desconsidera o email anterior com inscrição RBRAS.
by Carlos Pereira Consultoria Estatística 07 Mai '12

07 Mai '12

Favor desconsidera o email anterior com inscrição RBRAS. Carlos Pereira da Silva Bacharel em Estatística – UEPB Mestrando em Estatística e Experimentação Agropecuária - UFLA Lavras - MG Cel: 035 - 9119 - 2245 Email: ccpsilva81(a)hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8

1 0

Dúvida no uso de 2 versões do R
by Marco Pollo 06 Mai '12

06 Mai '12

Caros, Utilizo atualmente o R versão 2.14.1 na plataforma windows 64bits. No entanto, preciso utilizar um pacote que só roda em versões anteriores à versão 2.10("Erro: package ‘esamp2’ was built before R 2.10.0: please re-install it"). Posso utilizar duas versões diferentes sem prejuízo do uso que faço da versão mais atual? atencipo meus agradecimentos. Marco Pollo Almeida

2 6

Dúvida básica regressão
by Peterson Gardin 04 Mai '12

04 Mai '12

Pessoal Estou orientando um TCC de uma aluna que não é da minha área. Eu quando trabalho com regressão só uso a polinomial, sou agrônomo. Ela está fazendo um TCC que usa doses (0,1, 1, 10 ,100, 1000) de fungicidas (5 marcas de fungicidas) e 5 repetições. Devo extrair o log10 de doses e analisar por regressão linear?? as doses ficariam -1, 0,1, 2, 3 ou devo usar um modelo com log?? Não conheço esses modelos. Obrigado desde já. João Peterson Pereira Gardin Doutor em Fisiologia Vegetal, Eng. Agr., Pesquisador Fones: (49) 3566-0054 e 9933-4290 MSN: petersongardin(a)yahoo.com.br - SKYPE: petersongardin Endereço: Rua João Zardo, 1660, Bairro Campo Experimental, Videira - SC - CEP 89560000

3 3

Loop varrendo quantidade enorme de arquivos
by Thiago Veloso 04 Mai '12

04 Mai '12

Boa noite pessoal, Preciso de ajuda para fazer um loop que percorra uma lista enorme de arquivos. Organizei uma série (+50GB) de imagens de satélite que está organizada em anos (2001 a 2011) e, por ano, a cada oito dias (julianos). lai200101.tif lai200109.tif lai200117.tif lai200125.tif . . . lai2011361.tif onde a nomenclatura segue o padrão: lai{ano}{diajuliano}.tif. Percebam que o incremento entre os dias é sempre oito. Preciso fazer duas operações nesse banco de dados: 1) um loop para abrir cada um dos arquivos, nomeando a variável com o ano e mês, e multiplicar a mesma por 0.001. Algo assim: library(raster) lai{ano}{diajuliano}<-raster('lai{ano}{diajuliano}.tif') #carregando imagem lai{ano}{diajuliano}<-lai{ano}{diajuliano}*0.001 #calibrando imagem 2) Ao longo de todos os arquivos, somar cada quatro para produzir um mês de dados, por exemplo: lai_jan_2001 <- lai200101.tif + lai200109.tif + lai200117.tif + lai200125.tif lai_feb_2001 <- lai200133.tif + lai200141.tif + lai200149.tif + lai200157.tif . . . lai_dec_2011 <- lai2011337.tif + lai2011345.tif + lai2011353.tif + lai2011361.tif Um verdadeiro desafio de programação, concordam? Alguém pode me dar uma ideia de como começar? Agradece desde já, Thiago.

4 5

Predição de dados - Modelagem
by Sérgio Henrique almeida da silva ju 04 Mai '12

04 Mai '12

Olá Rodei um modelo: reg1<-lm(est ~ e1 + f1 + g1 + i1 + k1 + m1 + o1 + r1 + s1 + t1, data = dados) summary (reg1)$coefficients Estimate Std. Error t value Pr(>|t|) (Intercept) 33.811287 0.7983445 42.351753 1.557462e-202 e1 -3.300681 1.0234877 -3.224934 1.313683e-03 f1 -1.602262 0.6537985 -2.450697 1.448143e-02 g1 -1.337623 0.6721036 -1.990203 4.692499e-02 i1 -6.224414 0.9810556 -6.344608 3.820429e-10 k1 -1.279617 0.5544328 -2.307975 2.126718e-02 m1 -1.569079 0.5133349 -3.056638 2.316421e-03 o1 2.984603 1.6371984 1.822994 6.869551e-02 r1 1.340099 0.6338281 2.114295 3.481391e-02 s1 1.598583 0.4717877 3.388352 7.392205e-04 t1 6.286149 0.9227351 6.812518 1.946057e-11 E gostaria de criar uma variável em um outro banco (dados2) com a estimação da minha variável gerada pelo meu modelo, estou usando o seguinte comando, porém dá um erro. dados2$estim=predict(idfu.reg1, dados2, se.fit = F) Error in `$<-.data.frame`(`*tmp*`, "esqui", value = list(fit = c(37.4586673781906, : replacement has 4 rows, data has 5565 Como posso fazer isso? -- 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 4

Sobre distância euclidiana
by Marcio Leles 03 Mai '12

03 Mai '12

Calculo de distância euclidiana Bom dia. Tenho um arquivo que ajusto um modelo de regressão gerando três betas para cada árvore. Tenho casos de ter a medição de mil árvores. Após fazer os ajustes passo a ter beta 0, beta 1 e beta 2 para cada árvore. Gostaria de calcular a distância euclidiana do beta 0 com o beta 0, beta 1 com o beta 1 e beta 2 com o beta 2, entre todas as árvores. O meu objetivo é ter uma matriz de distância euclidiana entre todos os betas de todas as árvores para saber por exemplo qual árvore é mais próxima da árvore de número 21. Abaixo um pedaço do arquivo. Obrigado pela ajuda. ARV y x1 x2 coefb_0 coefb_1 coefb_2 1 2,631605 0,011628 0,000135 2,443489 -7,54249 6,537685 1 1 0,151163 0,02285 2,443489 -7,54249 6,537685 1 0,830123 0,267442 0,071525 2,443489 -7,54249 6,537685 2 1,531406 0,007576 5,74E-05 1,341076 -2,17855 0,89864 2 1 0,098485 0,009699 1,341076 -2,17855 0,89864 2 0,832656 0,174242 0,03036 1,341076 -2,17855 0,89864 3 1,645648 0,007194 5,18E-05 1,397167 -2,65584 1,382927 3 1 0,093525 0,008747 1,397167 -2,65584 1,382927 3 0,808183 0,165468 0,02738 1,397167 -2,65584 1,382927 4 2,08642 0,007042 4,96E-05 1,652506 -3,76968 2,395018 4 1 0,091549 0,008381 1,652506 -3,76968 2,395018 4 0,840278 0,161972 0,026235 1,652506 -3,76968 2,395018 -- Atenciosamente; Marcio Leles

3 2

Write.table Posicional
by Vinicius Brito Rocha 03 Mai '12

03 Mai '12

Senhores, Estou precisando exportar uma tabela onde não há delimitador e as variáveis são determinadas através de posições específicas.. Ex: Var1 (ocupa a posição 1-10 ) Var2 (ocupa a posição 11-31) e assim por diante. A função write table permite fazer algo assim ou apenas exporta com delimitador? Abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

3 3

Problemas com R 2.15
by Tito Conte 03 Mai '12

03 Mai '12

Acabei de instalar o R para Ubuntu e não estou conseguindo instalar nenhum pacote novo, fiz isso para o pacote ncdf, e tentei pra mais alguns e nada todos dão o mesmo erro (achava que era problema com o NetCDF) a mensagem que aparece segue abaixo installing *source* package ‘ncdf’ ... ** package ‘ncdf’ successfully unpacked and MD5 sums checked checking for nc-config... /usr/bin/nc-config configure: creating ./config.status config.status: creating src/Makevars ** libs gcc -std=gnu99 -I/usr/share/R/include -DNDEBUG -I/usr/include -fpic -O3 -pipe -g -c ncdf.c -o ncdf.o ncdf.c:3:20: fatal error: netcdf.h: No such file or directory compilation terminated. make: *** [ncdf.o] Error 1 ERROR: compilation failed for package ‘ncdf’ * removing ‘/home/tito/R/i686-pc-linux-gnu-library/2.15/ncdf’ The downloaded source packages are in ‘/tmp/Rtmpn8hiQ4/downloaded_packages’ Warning message: In install.packages("ncdf") : installation of package ‘ncdf’ had non-zero exit status Obrigado Tito

4 4