Junho 2013 - R-br - DINF UFPR

Converter topodata tifcom resolução espacial de 90m para 30m
by ASANTOS 21 Jun '13

21 Jun '13

Bom dia Membros, Baixei uma imagem TOPODATA no sitio do INPE referente a elevação horizontal de uma região em Minas Gerais, mas a resolução espacial é de 90m, para eu criar um slope para comparar com uma imagem de landsat preciso padronizar a resolução espacia para 30m, no meu caso obtive a imagem, sendo: # temp <- tempfile()## Objeto temporario download.file("http://www.webmapit.com.br/inpe/topodata/download/?ms=137174069031256&fbn=1…",temp)### Obtendo a imagem de elevação horizontal r1<- raster(c("19S435HN.tif"))## Coletado dentro da pasta temporária, pois não foi possivel fazer unz() para tif image(r1) E agora tenho que usar algum interpolador e pergunto para alguém que já trabalhou com isso, qual técnica utilizar, inverso do quadrado da distancia, krigagem ou outra. Pergunto também se alguém teria algum exemplo para postar, pois imagino que vou ter o mesmo problema que em alguns posts recentes em extrair a coordenadas dos pixels e criar uma coluna Data com o valor dos pixels. Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia 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 ======================================================================

3 3

erro com teste de wilcoxon
by alanarocha＠sapo.pt 21 Jun '13

21 Jun '13

bom dia, alguém sabe ajudar-me a padronizar estes dados em coluna? Casos (nº de doente) DoseTH9 DoseTH6 DoseTH3 DoseTC3 DoseTC6 DoseTC9 54020104 930 900 900 390 1.170,00 2.700,00 54020178 1.860,00 930 900 420 840 630 54020303 1.860,00 2.790,00 2.700,00 1.620,00 1.440,00 1.530,00 54020492 930 930 900 420 840 720 54020637 1.080,00 480 870 2.970,00 2.700,00 2.520,00 54020698 2.790,00 1.590,00 930 390 720 2.340,00 54020714 930 930 900 390 390 390 54030032 1.860,00 750 780 720 780 720 54030048 3.720,00 3.720,00 3.600,00 840 840 840 54030135 780 420 60 390 390 420 54030164 2.700,00 2.790,00 2.790,00 720 780 360 54030208 930 930 900 420 420 360 54030297 2.700,00 3.720,00 3.720,00 1.170,00 1.170,00 540 54030300 1.860,00 930 900 420 420 720 54030303 930 930 1.800,00 420 750 780 54030338 930 900 930 390 360 390 54030344 1.680,00 1.860,00 1.860,00 780 780 840 54030349 1.860,00 1.800,00 1.800,00 840 780 780 54030358 360 390 270 150 150 60 e como corrigir o erro: wilcox.test(dados[,-c(1,2,5,6,7)]$DoseTH3, dados[,-c(1,2,5,6,7)]$DoseTC3, paired=TRUE) Cumprimentos Ana

2 1

Re: [R-br] ggplot2
by Thiago V. dos Santos 20 Jun '13

20 Jun '13

Por favor mantenha a discussão na lista. Isso pode ajudar futuramente pessoas que tenham a mesma dúvida. Se eu entendi bem, os seus dados compreendem idades que variam de 34 a 75 anos. Porém, você parece querer mostrar, digamos, de 0 a 75. Se for realmente isso, eu usaria expand.grid. Veja esse exemplo *auto-reprodutível*: require(ggplot2) # Aqui eu produzo um conjunto de dados sintéticos # Repare que as idades variam de 34 a 75 dat <- data.frame(age=c(34:75), y=rnorm(42)) # Agora, eu expando a grade para o intervalo desejado, hipoteticamente 0 a 75 nesse caso all <- expand.grid(age = 0:75) # E em seguida junto com os originais dat <- merge(dat, all, all.y = TRUE) # Agora é fácil criar o gráfico p <-ggplot(dat) p + geom_line(aes(age,y)) Caso não seja isso, por favor formule melhor a pergunta e mostre um conjunto hipotético de dados que imitem o seu conjunto real. Saudações, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/ind… Phone: (612) 323 9898 ________________________________ From: geovane barbosa <geovanecb(a)yahoo.com.br> To: "thi_veloso(a)yahoo.com.br" <thi_veloso(a)yahoo.com.br> Sent: Thursday, June 20, 2013 8:16 PM Subject: Re: [R-br] ggplot2 Eu amigo minha idade varia de 34 a 75 anos. Quero todas as idades. Enviado do Yahoo! Mail no Android ________________________________ From: Thiago V. dos Santos <thi_veloso(a)yahoo.com.br>; To: r-br(a)listas.c3sl.ufpr.br <r-br(a)listas.c3sl.ufpr.br>; Benilton Carvalho <beniltoncarvalho(a)gmail.com>; geovane barbosa <geovanecb(a)yahoo.com.br>; Subject: Re: [R-br] ggplot2 Sent: Thu, Jun 20, 2013 11:45:42 PM Geovane, Como está a escala atual: idade em anos? Anos desde alguma data? Está em formato de data? E que tipo de modificação você quer fazer: mostrar cada 5 ou 10 anos? Em formato de inteiro ou em formato dd-mm-aaa? O excerto do código (ou a mensagem anterior) não ajuda a entender o objetivo. Saudações, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/ind… Phone: (612) 323 9898 ________________________________ From: geovane barbosa <geovanecb(a)yahoo.com.br> To: Benilton Carvalho <beniltoncarvalho(a)gmail.com> Cc: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Sent: Thursday, June 20, 2013 6:13 PM Subject: Re: [R-br] ggplot2 Olá desculpe c <- ggplot(dados, aes(y=dados$SA, x=dados$idade, colour=factor(dados$raca))) c + stat_smooth(size = 1)+xlab("Idade")+ylab("")+scale_colour_hue(name="") #################################### #################################### Prof. Geovane Carlos Barbosa UCL - Faculdade do Centro Leste ################################### ################################### ________________________________ De: Benilton Carvalho <beniltoncarvalho(a)gmail.com> Para: geovane barbosa <geovanecb(a)yahoo.com.br> Cc: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Enviadas: Quinta-feira, 20 de Junho de 2013 20:03 Assunto: Re: ggplot2 Sim... como nao ha' um exemplo minimo reproduzivel (explicando mais claramente seu objetivo), tbm nao ha' codigo para demonstrar como resolver o problema. Em 20 de junho de 2013 19:43, geovane barbosa <geovanecb(a)yahoo.com.br> escreveu: > > Olá pessoa boa noite, estou com dúvida no ggplot2 > > Tenho no eixo x a variável idade, porém como posso modificar a sua escala? > > abraços > > > > > > > > > > #################################### > #################################### > Prof. Geovane Carlos Barbosa > UCL - Faculdade do Centro Leste > ################################### > ################################### _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

Correlações
by Adriano Carvalho Costa 20 Jun '13

20 Jun '13

Boa noite, Estou estudando sobre correlações, alguém ai tem material sobre correlação: - Variáveis qualitativo e variáveis qualitativa - Variáveis discretas e variáveis discretas - Variáveis discretas e variáveis qualitativas - Variáveis qualitativas e variáveis continuas - Variáveis discretas e variáveis contínuas Obrigado!

1 0

sobrescrito
by Peterson Gardin 20 Jun '13

20 Jun '13

Olá pessoal 1) Como colocar subscrito para fórmula de nitrato NO3, colocar o 3 em subescrito. 2) Como colocar sobrescrito no g-1 e h-1. Comandos gráficos abaixo. bar.group(comparison$groups,horiz=FALSE,col=rainbow(5),border="blue", xlim=c(0,6), ylim=c(0,2.8),las=1) title(cex.lab=1.2, xlab="Porcentagem de N", ylab="Micromol NO3 g-1 h-1") João Peterson Pereira Gardin Doutor em Fisiologia Vegetal, Eng. Agr., Pesquisador Fones: (49) 3533 5626 e 9933-4290 MSN: petersongardin(a)yahoo.com.br - SKYPE: petersongardin Endereço: Rua João Zardo, 1660, Bairro Campo Experimental, Videira - SC - CEP 89560-000

4 7

Comparar dois data frames
by Fátima Lima Paula 20 Jun '13

20 Jun '13

Olá pessoal. Existe alguma forma de comparar dois data frames para saber se há alguma diferença entre eles? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

3 6

Fábio, please add me to your LinkedIn network
by MÁRCIO SOUZA 20 Jun '13

20 Jun '13

LinkedIn ------------ MÁRCIO SOUZA requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from MÁRCIO SOUZA http://www.linkedin.com/e/-33g73l-hi6et0wh-1d/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View profile of MÁRCIO SOUZA http://www.linkedin.com/e/-33g73l-hi6et0wh-1d/rso/265728365/Zfs1/name/82777… ------------------------------------------ You are receiving Invitation emails. This email was intended for Fábio Corrêa. Learn why this is included: http://www.linkedin.com/e/-33g73l-hi6et0wh-1d/plh/http%3A%2F%2Fhelp%2Elinke… (c) 2012, LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 20 Jun '13

20 Jun '13

Aos membros do R-br por favor, voltando ao meu propósito inicial eu preciso apenas de uma orientação em relação a uma função que retorne o índice (i,j) e o valor do pixel de todo raster no R, nada mais do que isto. O R possui bastante recursos para executar esta tarefa, mas infelizmente não encontrei solução. Depois disto exportar/converter para shapefile será mais simples. Grato, Alex Santos

3 3

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 20 Jun '13

20 Jun '13

Elias, informo que minhas amostas não são pixels. Minha metodologia está consolidada, preciso apenas de uma diretriz para encontrar está função no R (pois acredito que o mesmo seja capaz para realizar o que desejo. Abaixo, está em linha gerais minha metodologia, mas por favor, não se atente a mesma pois meu foco é encontra/contrair a função que extraia o índice (i,j) de cada pixel. 1 - Eu tenho uma articulação do mapeamento sistemático que recobre minha área de interesse (no formato shapefile - poligono). Sobre estas áreas que houve a amostragem sistemática. (obs: Numa solução paralela utilizando o FME (software comercial) já agreguei as informações de (i,j) em cada quadrícula, podendo utilizar esta informação como X, Y e representa minhas quadrículas de forma contínua. Mas quero substituir o FME (comercial) pelo R (livre), pois acredito no potencial do R e acho ele capaz de executar esta tarefa. 2 - converti a articulação para raster, já com a proposta de encontrar as informações do índice (i,j) de cada pixel. Sendo que cada pixel possui o mesmo tamanho de minhas quadrículas. Mas espero que compreenda, pois o que desejo realmente é uma função que leia um arquivo raster - extraia o índice (i,j) de cada pixel (compondo uma matriz com os índices mais o valor do atributo) e exporte esta matriz para um shapefile. obs.: Adianto que o mesmo não será georreferenciado e nem há necessidade neste momento, pois haverá o georreferenciamento após a krigagem através dos extremos da articulação original. Att, Alex Santos

2 1

sortear 100x e correlacionar
by Hélio Gallo Rocha 20 Jun '13

20 Jun '13

Bom dia a todos, Peço ajuda na seguinte questão destes dados que comparam dois modelos de previsão de safra em café, onde : colheita = dados da parcela colhida modelo 1 de previsão modelo 2 de previsão CRM coleta=c(rep(1,25),rep(2,25)) colheita=c(10.2,9.22,8.6,15.57,11.93,12.28,12.43,16.47,18.72,12.77,16.03,13.4,16.95, 15.4,18.64,14.38,15.02,16.98,14.38,16.35,18.28,16.06,17.92,12.85,13.83,4.48,3.35,7.24, 7.73,8.45,4.83,3.8,2.2,5.99,4.24,4.67,8.31,4.25,7.32,5.45,4.24,5.72,5.09,4.84,5.62,5.5, 4.54,4.72,7.69,3.45) modelo1=c(17.08,12.24,13.44,18.52,11.82,11.37,16.17,16.92,16.66,13.32,18.95,14.78, 15.54,16.12,25.62,18.2,15.86,20.5,20.47,20.45,17.15,22.15,20.29,13.82,12.86,3.13, 1.84,6.9,6.35,8.9,5.64,4.79,2,4.64,5.06,4.85,5.61,4.82,6.59,4.15,4.46,5.1,4.78, 4.05,5.02,5.65,3.04,3.61,9.48,4.13) modelo2=c(12.43,10.99,10.32,12.46,9.85,10.15,12.68,12.46,12.16,9.92,12.58,10.63, 11.18,10.6,15.71,13.22,10.82,12.85,13.48,14.13,11.41,14.09,12.23,10.28,8.75, 5.2,3.56,9.89,9.03,12,8.16,7.15,3.2,6.58,6.84,7.37,8.5,7.87,8.38,5.58,6.62, 6.36,6.84,6.48,7.58,7.35,4.99,5.3,11.99,6.21) dados=cbind(coleta,colheita,modelo1,modelo2);dados x=cbind(colheita,modelo1,modelo2) cor(x) #sorteio de 20 medidas, r=do.call(rbind, by(dados, dados[,1], function(x, ...) x[sample(1:nrow(x), ...), ], size=10)) y=cbind(r$colheita,r$modelo1,r$modelo2);y cor(y) Gostaria de sortear 100 vezes com 20 amostras, sendo o resultado a média da correlação e a média dos dados Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 4